Adding Constraint Equations
In practice, one may not have a geophone at shot B, or there may be geophones without observed first arrival times. In those cases, one must add constraint equations to the bottom of the matrix, G. For example, if we were unable to pick a time for the signal from Shot A to geophone 2, we might add a constraint that sets delay time at geophone 1 equal to that at geophone 2 (ie.
). That system would look like this:
|
(27) |
In this example, we set equal to with a weight of 9 fold in the least squares solution. This is the last row of matrix in equation 27. It is like the following equation 28 being repeated 9 times in the matrix equation 27:
|
(28) |
The effect is to strongly weight these two delay times equal. The OLS solution to Equation 27 is again given by
|
(29) |
where, in the case of Equation 27, the vector, , contains the delay times for the shots, receivers, and the refractor slowness.
|
(30) |
If we want to constrain the refractor velocity there is a way to do that too. Perhaps we lack reverse profiles or a second shot with at least a different angle between the shots and geophones which would yield more than one apparent refractor velocity. Modifying the above situation to include only shot A with a constrained refractor velocity, and a constrained delay time, , we could set up a system like this:
|
(31) |
Here, we set the velocity of the refractor at m/s and constrain the shot A delay time to seconds. Both of these values would need to come from some other source, like an intersecting survey that was absent the defects here.