In most practical situations, the traveltime curves for recorded waves can be described analytically as a family of functions that rely on several variables and parameters:

$T=F\left(x_1,x_2,\cdots ,c_1,c_2,\cdots ,c_i,\cdots \right),$ (1)

where $x_1,x_2,\cdots$ are variables and $c_1,c_2,\cdots ,c_i,\cdots$ are independent parameters.

The desired traveltime curves for transformed waves can also be described as a family of functions that are dependent on a several of variables and parameters:

${\tilde{T}}_0=P\left(y_1,y_2,\cdots ,d_1,d_2,\cdots ,d_j,\cdots \right),$ (2)

where $y_1,y_2,\cdots$ are variables and $d_1,d_2,\cdots ,d_j,\cdots$ are independent parameters.

If at least one of these parameters is the same for both families $d_j=c_i$, we can express $d_j$ using the remaining terms of the Equation (2)

$d_j=R\left({\tilde{T}}_0,y_1,y_2,\cdots ,d_1,d_2,\cdots ,d_{j-1},d_{j+1},\cdots \right)$ (3)

and substitute it into Equation (1) to obtain the new family of functions

$\tilde{T}=F\left(R\left({\tilde{T}}_0,y_1,y_2,\cdots ,d_1,d_2,\cdots ,d_{j-1},d_{j+1},\cdots \right),x_1,x_2,\cdots ,c_1,c_2,\cdots ,c_{i-1},c_{j+1},\cdots \right)\text{,}$ (4)

that connects families (1) and (2).

The envelope of the family (4) $\widetilde{\stackrel{˜}{T}}$ can be obtained by excluding any parameters of this family using the necessary condition for the existence of the envelope [32] [36] and so on:

${{\tilde{T}}^{\prime }}_{r_1}=0,\cdots ,{{\tilde{T}}^{\prime }}_{r_k}=0\text{,}$ (5)

where $r_k$ is an independent parameter of the family (4).

To perform stationary-phase transformations, seismic data must be integrated along the envelope of the family (4) and the outcome assigned to point ${\tilde{T}}_0$ of the some resultant trace located at point $\xi$. In this scenario, there exists a certain value of ${\tilde{T}}_0$ at which the condition $d_j=c_i$ is met. Consequently, one of the functions from family (1) will be coincide to one of the functions from family (4), and the envelope of the family (4) will be in touching with this function. During integration seismic data along $\widetilde{\stackrel{˜}{T}}$, the vicinity of the point of tangency is shifted to the desired location.

It is proposed that the difference, which is named by us “touch character”, between the second derivatives at the point of contact be kept constant during the integration. In this case, the wave distortions that unavoidably occur as a consequence of the transformations will be almost identical and will have little influence on the outcome of future correlation of the transformed waves. The concrete value of “touch character” must be determined experimentally.

If we denote “touch character” as $\alpha$ then:

$\alpha ={\widetilde{\stackrel{˜}{T}}}^{″}-T^{″}=const\text{.}$ (6)

If seismic conditions are almost similar, $\alpha$ is a very stable parameter that, once established, can be used for various transformations.

The “touch character” constancy can be established during the transformation process by uneven summing of seismic records. For certain simple stationary-phase transformations, the necessary summation order can be determined analytically using the difference of total differentials or a second-order Taylor series expansion around the point of tangency.

In most circumstances, however, an approximate technique is utilized, which consists of the stages listed below. Assuming that some point on the summation line (the envelope of the family of functions (4)) represents the point of tangency, will choose the next point so that the time on the summation (integration) curve varies from the predicted time on the traveltime curve by the value of “touch character”. In this case, the summation must be performed in sufficiently small steps to prevent missing the desired point.

The impact of the value of the “touch character” on the quality of the transformation is shown in Figures 1-4 .

Figure 1 depicts a fragment of the CMP timesection. Figures 2-4 depict the same fragment after prestack time migrations with “touch character” values of 5000, 500, and 10, respectively. According to the authors, the transformation resolution in Figure 2 is high, but noise immunity is low. Figure 4 depicts the opposite situation-low resolution and high noise immunity. Only in Figure 3 do we see an acceptable compromise between them.

The stationary-phase transformations are implemented using the same scheme for any transformation. First, equations for families of functions that describes the traveltime curves for waves recorded and sought after the transformation should be found. At least one parameter from each family should be shared. We express

this parameter using the parameters and variables from the second equation and then plug the resulting expression in the first equation.

The new family of functions connects both preceding families. Then should be find the envelope of this family and summarize the recorded seismic data along it unevenly. The order of summation is determined by the difference of total differentials or the decomposition of functions according to the Taylor formula, or approximately as shown above.

So, the sum of irregularly located seismic data is used to perform stationary-phase transformations. Because of this unevenness, the difference in the second derivatives at the point of tangency between the wave’s traveltime curve and the summation line may be preserved as a constant, allowing for a compromise between wave transformation resolution and noise immunity.

The technique of stationary-phase transformations will be illustrated in more detail below, using the example of an enough simple and extensively used procedure for transforming waves reflected from a flat boundary and recorded on time sections into time migrated sections.

The “touch character” value is determined experimentally and, to some extent, subjectively.

Time and depth migrations techniques are widely accepted. The works of [2] [3] [5] [15] [17] [18] and others played and are playing a significant role in the formation and promotion of this approach.

The corresponding procedure created by the stationary-phase transformation method is quite simple and serves as a good illustration of the implementation of this method.

The equation of the family of functions describing the reflection traveltime curves for the average velocity model with plane reflectors recorded on time sections has the form

$T=T_0+t*\eta \text{,}$ (7)

where $\eta$ is the location of the time section trace ( Figure 5 ).

Following the formulation of the problem, the transformation’s result must be located at the point $\left({\tilde{T}}_0,\xi \right)$, where ${\tilde{T}}_0$ is the double vertical depth in the time scale and $\xi$ is the location of the migrated section trace.

Because the cosine of the reflector dip angle is equal to $\sqrt{4-t^2*V^2}/2$ ( [32] [37] and so on) the equation of the family of traveltime curves for the reflected waves on migrated time sections is:

${\tilde{T}}_0=2*\left(T_0+t\xi \right)/\sqrt{4-t^2*V^2}\text{.}$ (8)

On Figure 5 for $T_0=2000$ ms, $t=0.2$ ms/m and $V=1.5$ m/ms we can see some functions of the family of the reflection traveltime curves on time section and them also on time migrated section.

Here, we express the parameter $T_0$ in terms of the rest of the family’s parameters and variables of the family (8):

$T_0={\tilde{T}}_0*\sqrt{4-t^2*V^2}/2-t\xi$ (9)

and substitute it into (7):

$\tilde{T}={\tilde{T}}_0*\sqrt{4-t^2*V^2}/2+t\left(\eta -\xi \right)\text{.}$ (10)

The acquired family connects the family of traveltime curves for reflected waves (7) with the desired family of traveltime curves for transformed waves (8) on time migrated sections.

It is easy to make sure that, if the value of ${\tilde{T}}_0$ is (8), then the equation of the function of the family $\tilde{T}$ equal $\tilde{T}=T_0+t*\eta$ regardless of the concrete value of $\xi$ and the envelope of the family (10) will be in touching with this function.

To find the envelope of the family (10) let’s use the condition ${T^{\prime }}_t=0$ for excluding the parameter $t$:

$t=4*\left(\eta -\xi \right)/\left(V*\sqrt{{\tilde{T}}_0^2*V^2+4*{\left(\eta -\xi \right)}^2}\right)\text{.}$ (11)

If $t$ is substituted into (10), the envelope of the family (10) will be obtained

$\widetilde{\stackrel{˜}{T}}=\sqrt{{\tilde{T}}_0^2+4*{\left(\eta -\xi \right)}^2/V^2}\text{.}$ (12)

For $t=0.2$ ms/m and $\xi =2500$ m on Figure 6 you can see some functions of the family (10) and coincidence of the function $T=2000+0.2*\eta$ from the family (7) with a function from the family (10) and touching this function with the envelope (12).

To satisfy condition (6) let’s expand (7) and (12) in a Taylor series at the point $\eta$ and yields:

$\begin{array}{l}{T^{\prime }}_{\eta }=t\text{,}\\ {T^{″}}_{{\eta }^2}=0\text{,}\end{array}$ (13)

${\widetilde{\stackrel{˜}{T^{\prime }}}}_{\eta }=\frac{4*\left(\eta -\xi \right)}{V^2*\sqrt{{\tilde{T}}_0^2+\frac{4*{\left(\eta -\xi \right)}^2}{V^2}}},$ (14)

${{{\widetilde{\stackrel{˜}{T}}}^{\prime }}^{\prime }}_{{\eta }^2}=\frac{4*{\tilde{T}}_0^2}{V^2*{\left(\sqrt{{\tilde{T}}_0^2+\frac{4*{\left(\eta -\xi \right)}^2}{V^2}}\right)}^3}\text{,}$ (15)

$T\left({\eta }_i+\Delta {\eta }_i\right)={\tilde{T}}_0*\sqrt{4-t^2*V^2}/2+t\left({\eta }_i-\xi \right)+\Delta {\eta }_i*t,$ (16)

$\begin{array}{c}\widetilde{\stackrel{˜}{T}}\left({\eta }_i+\text{Δ}{\eta }_i\right)=\sqrt{{\tilde{T}}_0{}^2+4*{\left({\eta }_i-\xi \right)}^2/V^2}+\frac{4*\Delta {\eta }_i*\left({\eta }_i-\xi \right)}{V^2*{\sqrt{{\tilde{T}}_0{}^2+4*{\left({\eta }_i-\xi \right)}^2/V^2}}^3}\\ +2*\frac{{\left(\text{Δ}{\eta }_i\right)}^2*{\left({\tilde{T}}_0\right)}^2}{V^2*{\sqrt{{\tilde{T}}_0{}^2+4*{\left({\eta }_i-\xi \right)}^2/V^2}}^3}\text{.}\end{array}$ (17)

If subtract (16) from (17), given that ${\eta }_i$ is the point of tangency, and equates this difference to $\alpha =const$, we will get:

$\text{Δ}{\eta }_i=\sqrt{\alpha *V^2*{\sqrt{{\tilde{T}}_0^2+4*{\left({\eta }_i-\xi \right)}^2/V^2}}^3}/\left(2*{\tilde{T}}_0\right).$ (18)

Therefore, to transform time section into time migrated section, we should sum the origin’s traces along (12) and put the results into the point $\left({\tilde{T}}_0,\xi \right)$ for all values, ${\tilde{T}}_0,{\eta }_i$ and $\xi$ ( Figure 7 ).

The Equation (12) is the equation of the diffracted wave with the diffraction center at point $\left({\tilde{T}}_0,\xi \right)$. Thus, the stationary-phase method for the average-velocity model of the medium in kinematic terms agrees with the Kirchhoff method and differs from it in dynamic features.

From Figure 7 we can see that the equation of the reflecting boundary is described as

$z=h_0+k*\xi ,$ (19)

where $h_0$ is vertical distance from the source to the reflector and $k=tg\phi =t*V/\sqrt{4-t^2*V^2}$.

Let’s use (19) to build the depth migration procedure.

If we take into account that

$h_0=T_0*V/\sqrt{4-t^2*V^2}\text{,}$ (20)

then the equation of the flat reflector is:

$z=\left(T_0+t*\xi \right)*V/\sqrt{4-t^2*V^2}\text{.}$ (21)

Let’s represent parameter $T_0$ by utilizing other parameters and variables from the family (21):

$T_0=z*\sqrt{4-t^2*V^2}/V-t*\xi$ (22)

and substitute it into (7):

$\tilde{T}=\left[z*\sqrt{4-t^2*V^2}+t*V*\left(\eta -\xi \right)\right]/V.$ (23)

For finding the envelope of this family excluding parameter $t$, using the next condition:

$t=2*\left(\eta -\xi \right)/\left[V*\sqrt{z^2+\left(\eta -\xi \right)}\right]$ (24)

${{\tilde{T}}^{\prime }}_t=-z*t*v/\sqrt{4-t^2*V^2}+\left(\eta -\xi \right)=0\text{,}$ (25)

Substituting (25) into (23) to obtain the equation of the envelope of the family (23):

$\widetilde{\stackrel{˜}{T}}=2*\sqrt{z^2+{\left(\eta -\xi \right)}^2}/V,$ (26)

which will be used for the receiving depth migration section.

As can be seen from Figure 7 , the envelopes used to realize time and depth migration coincide, but the results of summation along these envelopes refer to different points ${\tilde{T}}_0$ and $z$.

To keep the “touch character” during summation approximately constant, let us expand functions (7) and (26) by Taylor series at the point of tangency by analogy with (13)-(17).

In this case we will obtain

$\Delta {\eta }_i=\sqrt{\alpha *V*{\sqrt{z^2+{\left({\eta }_i-\xi \right)}^2}}^3}/2.$ (27)

The result of testing this procedure is shown in Figure 8 .

The traveltime curve for the wave reflected from a flat reflector in the CMP system is described by a hyperbola of the following shape [32] [37] :

$T=\sqrt{T_0^2+4*{\eta }^2/V_{stk}^2}\text{,}$ (28)

where $\eta$ is half of the source-receiver offset ( Figure 9 ).

Sometimes it is preferable to transform waves with hyperbolic traveltime curves into a wave with linear traveltime curves (for example, to increase the resolution of the CMP summation) of the following shape:

${\tilde{T}}_0=T_0+2*\xi /V_{stk}\text{,}$ (29)

where $\xi$ is the location of the result trace.

$T$ is the family of functions that depends on the parameters $\left(T_0,V_{stk}\right)$, and ${\tilde{T}}_0$ is the family of functions that also depends on parameters $\left(T_0,V_{stk}\right)$.

Let us express $T_0$ through other parameters and variables of the family (29):

$T_0={\tilde{T}}_0-2*\xi /V_{stk}$ (30)

and creates a new family of functions by substituting (30) into (28):

$\tilde{T}=\sqrt{{\left({\tilde{T}}_0-2*\xi /V_{stk}\right)}^2+4*{\eta }^2/V_{stk}^2}\text{.}$ (31)

The envelope of the family (31) is determined based on the Equation (5):

${{\tilde{T}}^{\prime }}_{(V_{stk})}=\left({\tilde{T}}_0-2*\xi /V_{stk}\right)*\xi -\left(2*{\xi }^2\right)/V_{stk}=0\text{;}{V}_{stk}=2*\left({\xi }^2+{\eta }^2\right)/\left({\tilde{T}}_0*\xi \right)\text{;}$ (32)

$\widetilde{\stackrel{˜}{T}}={\tilde{T}}_0*\eta /\sqrt{{\xi }^2+{\eta }^2}\text{.}$ (33)

Consequently, if we sum the recorded CMP traces along (33) for different values of the parameters and variables, we obtain a new CMP record with the traveltime curves for transformed waves which are linear functions ( Figure 10 ).

According to (33), values of $V_{stk}$ is not used directly to perform this procedure but their approximate values can be used to maintain the “touch character” during summing for constant. In this case we use the expansion of functions (28) and (33) in a Taylor series by analogy with (13)-(17) and obtain

$\text{Δ}{\eta }_i=\sqrt{\frac{\alpha }{\frac{2*{\left({\tilde{T}}_0-2*\xi /V_{stk}\right)}^2}{V_{stk}^2*{\left(\sqrt{{\left({\tilde{T}}_0-2*\xi /V_{stk}\right)}^2+4*{\eta }_i^2/V_{stk}^2}\right)}^3}-\frac{3*{\tilde{T}}_0*{\xi }^2*{\eta }_i}{2*{\left(\sqrt{{\xi }^2+{\eta }_i^2}\right)}^5}}}.$ (34)

If is no information about $V_{stk}$, the approximation technique mentioned in the introduction is employed to maintain the “touch character” as a constant.

The result of testing this procedure is shown in Figure 10 .

Authors didn’t find analogs for this procedure.

It is possible to transform time sections into time migrated sections using the values of $V_{stk}$ acquired in the process of realization of the CMP method ( Figure 11 ).

Due to the fact that the cosine of the reflector decline angle is equal $\sqrt{4-t^2*V^2}/2$

$V_{stk}=2*V/\sqrt{4-t^2*V^2}$ (35)

and considering that (8) is the equation of the family of traveltime curves for reflected waves on time migrated sections, we obtain the equation for family of reflection traveltime curves on time migrated sections when using the stack velocity:

${\tilde{T}}_0=\left(T_0+t\xi \right)*\sqrt{4-t^2*V_{stk}^2}/2\text{,}$ (36)

where $\xi$ is the location of the migrated section trace.

To obtain the equation of the family that connects the family for reflection traveltime curves on time sections with the family for reflection traveltime curves on time migrated sections the parameter $T_0$ from Equation (36) should be written in terms of the other parameters and variables of this equation and substituted into Equation (7):

$T_0=2*{\tilde{T}}_0/\sqrt{4+t^2*V_{stk}^2}-t*\xi$ (37)

and

$\tilde{T}=2*{\tilde{T}}_0/\sqrt{4+t^2*V_{stk}^2}+t*\left(\eta -\xi \right)\text{.}$ (38)

The necessary condition for the existence of an envelope of a family of functions yields the following equation:

$\left(\eta -\xi \right)=2*{\tilde{T}}_0*t*V_{stk}^2/{\sqrt{4+t^2*V_{stk}^2}}^3.$ (39)

To solve this equation, almost any numerical method may be utilised. One such strategy is Newton’s approach. In this case,

$F=\left(\eta -\xi \right)-2*{\tilde{T}}_0*V_{stk}^2/{\sqrt{4+t^2*V_{stk}^2}}^3=0,$ (40)

$t_{i+1}=t_i-F\left(t_i\right)/{F^{\prime }}_t\left(t_i\right)\text{,}$ (41)

${F^{\prime }}_t=-4*{\tilde{T}}_0*t*V_{stk}^2*\left(2-t^2*V_{stk}^2\right)/{\sqrt{4+t^2*V_{stk}^2}}^5=0,$ (42)

because the equation of the envelope of the family (38) cannot be written analytically, the approximation technique mentioned in the introduction is employed to maintain the “touch character” as a constant.

The result of testing this procedure is shown in Figure 12 .

Authors didn’t find analogs for this procedure.

Despite its undeniable benefits, the CMP method distorts the recorded wave field to some degree, similar to every other seismic data processing method. In many seismic-geological situations, it is preferable to acquire migrated sections immediately from an observed wave field.

Largely due to the works of [2] [3] [5] [6] [9] [15] [16] [18] - [20] [32] [38] - [40] and others, time and depth migrations have found wide application in the practice of seismic research.

To construct such a procedure using the stationary-phase method, we must first determine the equation of the family of traveltime curves for the waves reflected from the flat reflector for the average velocity medium model.

According to method of stationary-phase transformations of seismic data for implementing desirable transformations first of all we need to have equations of families of traveltime curves of reflected from a flat reflector in an average velocity medium registered on CMP seismograms and families of traveltime curves of reflected waves on time and depth migrated sections ( Figure 13 ).

The first of required equation can be obtained from [32] and [37] and so on, and the rest two we got earlier and they equal, respectively, (8) and (21). For the sake of completeness, we will provide all 3 equations.

$T=\sqrt{{\left(T_0+t*\eta \right)}^2+\left(4-t^2*V^2\right)*x^2/V^2}\text{,}$ (43)

${\tilde{T}}_0=2*\left(T_0+t\xi \right)/\sqrt{4-t^2*V^2}\text{,}$ (44)

$z=\left(T_0+t*\xi \right)*V/\sqrt{4-t^2*V^2}\text{.}$ (45)

When we express parameter $T_0$ in terms of the other parameters and variables of Equation (44)

$T_0={\tilde{T}}_0*\sqrt{4-t^2*V^2}/2-t\xi$ (46)

and substitute it into (43) we obtain the family of functions (47) that connects families (43) and (44).

$\tilde{T}=\sqrt{{\left[{\tilde{T}}_0*\sqrt{4-t^2*V^2}/2+t*\left(\eta -\xi \right)\right]}^2+\left(4-t^2*V^2\right)*x^2/V^2}\text{.}$ (47)

We’ll use condition ${{\tilde{T}}^{\prime }}_t=0$ to find the envelope of this family:

$\begin{array}{l}{{\tilde{T}}^{\prime }}_t=\left[{\tilde{T}}_0*\sqrt{4-t^2*V^2}/2+t*\left(\eta -\xi \right)\right]\\ *\left[\left(\eta -\xi \right)-{\tilde{T}}_0*t*V^2/\left(2*\sqrt{4-t^2*V^2}\right)\right]-t*x^2=0.\end{array}$ (48)

$t=\sqrt{2*\left(1-A/\sqrt{A^2+B}\right)}/V,$ (49)

where

$A=x^2-{\left(\eta -\xi \right)}^2+T_0^2*V^2/4\text{and}B=T_0^2*V^2*{\left(\eta -\xi \right)}^2.$

The envelope of the family (47) is equal

$\widetilde{\stackrel{˜}{T}}=\sqrt{{\tilde{T}}_0^2/4+{\left[x-\left(\eta +\xi \right)\right]}^2/V^2}+\sqrt{{\tilde{T}}_0^2/4+{\left[x-\left(\eta -\xi \right)\right]}^2/V^2}\text{.}$ (50)

Summing the seismic traces registered on the CMP seismograms along curve (50) results in a time migrated section ( Figure 14 ).

Because Equation (50) is the equation of a diffracted wave with the center of diffraction at point $\left({\tilde{T}}_0,\xi \right)$, the stationary-phase method for the average velocity model of the medium in kinematic terms is identical to the Kirchhoff method and differs from it in dynamic features.

The difference between complete differentials in tangent point is used to store “touch character” as a constant. These expressions are too extensive to fit in this study.

To obtain the depth migrated section, we need to use the equation of the family of functions (45).

By analogy (46)-(50) we express $T_0$ in terms of the other parameters and variables of the Equation (45) and substitute it into (43) and obtain the equation of the family of functions connecting these families:

$\tilde{T}=\sqrt{{\left[z*\sqrt{4-t^2*V^2}+t*V*\left(\eta -\xi \right)\right]}^2+\left(4-t^2*V^2\right)*x^2}/V.$ (51)

The envelope of this family is found using condition ${{\tilde{T}}^{\prime }}_t=0$ and by substituting found value of

$t=\sqrt{2*\left\{1-\frac{z^2+x^2-{\left(\eta -\xi \right)}^2}{\sqrt{4*z^2*{\left(\eta -\xi \right)}^2+{\left[z^2+x^2-{\left(\eta -\xi \right)}^2\right]}^2}}\right\}}/V$ (52)

into (51).

Representing this envelope explicitly is too cumbersome.

Figure 14 illustrates the outcomes of our test.

The approximation technique mentioned in the introduction is employed to maintain the “touch character” as a constant.

In the absence of information on the propagation average velocities of seismic waves, the use of this procedure to acquire a time migrated section appears sufficiently justifiable.

The equations necessary for the implementation of this procedure were obtained earlier (43) and (36), and they are presented below ( Figure 15 ):

The equation of families of traveltime curves of reflected from a flat reflector in an average velocity medium registered on CMP seismograms-

$T=\sqrt{{\left(T_0+t*\eta \right)}^2+\left(4-t^2*V^2\right)*x^2/V^2}\text{,}$ (53)

and the equation of family of reflection traveltime curves on time migrated sections when using the stack velocity:

${\tilde{T}}_0=\left(T_0+t\xi \right)*\sqrt{4-t^2*V_{stk}^2}/2\text{.}$ (54)

If we express parameter $T_0$ in terms of the other parameters and variables of Equation (54) and substitute it into (53) we obtain the equation of family of functions connecting (53) and (54) that has the following form:

$\tilde{T}=\sqrt{{\left[2*{\tilde{T}}_0/\sqrt{4+t^2*V_{stk}^2}+t*\left(\eta -\xi \right)\right]}^2+4*x^2/V_{stk}^2}.$ (55)

The envelope of this family will be found using condition ${{\tilde{T}}^{\prime }}_t=0$:

${{\tilde{T}}^{\prime }}_t=\left[2*{\tilde{T}}_0/\sqrt{4+t^2*V_{stk}^2}+t*\left(\eta -\xi \right)\right]*\left[\left(\eta -\xi \right)-\frac{2*{\tilde{T}}_0*t*V_{stk}^2}{\sqrt{4+t^2*V_{stk}^2}}\right]=0\text{.}$ (56)

This equation is being decided by using Newton’s method.

Because the equation of the envelope of the family (55) cannot be expressed analytically, the approximate method described in the introduction is utilized to keep the “touch character” as constant.

The result of our test is shown in Figure 16 .

Authors didn’t find analogs for this procedure.

The CMP approach for dipping reflectors, as is widely known [38] , includes records in the summing process that are unrelated to the common reflecting point. The DMO procedure [41] was created to address this problem to a certain extent.

The following approach is used to produce a time section that includes reflections from one point.

The equation of the family of traveltime curves of the reflected waves for a medium-velocity model with flat reflectors recorded on CMP seismograms equals (43), and (7) is equation of the family of traveltime curves for the same waves recorded on time sections ( Figure 17 ).

When we express parameter $T_0$ in terms of the other parameters and variables of Equation (7),

$T_0={\tilde{T}}_0-t*\xi ,$ (57)

where $\xi$ is the location of the resulting time section trace, and putting it into Equation (43) yields:

$\tilde{T}=\sqrt{{\left[{\tilde{T}}_0+t*\left(\eta -\xi \right)\right]}^2+\left(4-t^2*V^2\right)*x^2/V^2}.$ (58)

It is equation of the family of functions that connect families (43) and (7).

Condition

${{\tilde{T}}^{\prime }}_t={\tilde{T}}_0*\left(\eta -\xi \right)+t*{\left(\eta -\xi \right)}^2-t*x^2=0$ (59)

is used to find the parameter $t$:

$t={\tilde{T}}_0*\left(\eta -\xi \right)/\left[x^2-{\left(\eta -\xi \right)}^2\right]\text{.}$ (60)

The equation of the envelope used to convert reflections on CMP seismograms into reflections on time sections was obtained by substituting the value of $t$ into the equation of the family of the functions (58) and has the form:

$\widetilde{\stackrel{˜}{T}}=\sqrt{{\tilde{T}}_0^2*x^2/\left[x^2-{\left(\tau -\xi \right)}^2\right]+4*x^2/V^2}.$ (61)

Figure 18 displays the outcome of this procedure’s testing.

Authors didn’t find analogs for this procedure.

The following parametric equation describes the propagation of seismic waves in an inhomogeneous medium [42] - [44] :

$x=s-\underset{0}{\overset{z}{{\displaystyle \int }}}\frac{m*V\left(x,z\right)}{\sqrt{1-m^2*V^2\left(x,z\right)}}\text{d}z\text{,}t=\underset{0}{\overset{z}{{\displaystyle \int }}}\frac{\text{d}z}{V\left(x,z\right)*\sqrt{1-m^2*V^2\left(x,z\right)}}\text{,}$ (62)

where $s$ is the source point, $x$ is the abscissa of a wavefront point, $z$ is the ordinate of this point, and $t$ is the time of wave propagation from the point $\left(s,0\right)$ to the point $\left(x,z\right)$, $V\left(x,z\right)$ is the velocity of wave propagation in the point $\left(x,z\right)$, $i_z$ is the angle between the ray and vertical at a point $\left(x,z\right)$, and $m=\sin \left(i_z\right)/V\left(x,z\right)$ ( Figure 19 ).

If $V\left(x,z\right)$ is the linear function of $z$, then $V=V_0*\left(1+\beta *z\right)$ and equation (65) will have the following form [40] :

$\begin{array}{l}x=s-\frac{1}{m*V_0*\beta }*\left[\sqrt{1-m^2*V_0^2}-\sqrt{1-m^2*V_0^2*{\left(1+\beta *z\right)}^2}\right],\\ t=\frac{\frac{1}{V_0*\beta }*\ln \left(1+\beta *z\right)*\left(1+\sqrt{1-m^2*V_0^2}\right)}{1+\sqrt{1-m^2*V_0^2*{\left(1+\beta *z\right)}^2}}.\end{array}$ (63)

For a flat reflector $z=h+kx$ ( Figure 19 ) and a normal ray $m=\sin \left(\phi \right)/V\left(z\right)$, where $\phi$ is the dip angle of the reflector and $k=\text{tan}\left(\phi \right)$. If $t$ designated as $T$ and $s$ as $\eta$ , we have the following:

$x=\eta -\left(\sqrt{{\left(1+\beta *h+k*\beta *x\right)}^2*\left(1+k^2\right)-k^2}-\left(1+\beta *h+k*\beta *x\right)\right)/\left(\beta *k\right)\text{,}$(64)

$T=\frac{2}{\beta *V_0}*\text{ln}\frac{\left(1+\beta *h+k*\beta *x\right)*\sqrt{1+k^2}+\sqrt{{\left(1+\beta *h+k*\beta *x\right)}^2*\left(1+k^2\right)-k^2}}{1+\sqrt{1+k^2}}\text{.}$

(65)

Substituting variable x from Equation (63) into Equation (65) gives equation (66), which describes the traveltime curves for waves propagating at a velocity linearly dependent on depth and reflected from inclined boundaries:

$T=\frac{2}{\beta *V_0}*\text{ln}\frac{\left(1+\beta *h+k*\beta *\eta \right)+\sqrt{{\left(1+\beta *h+k*\beta *\eta \right)}^2*\left(1+k^2\right)-k^2}}{1+\sqrt{1+k^2}}\text{.}$ (66)

The equation of the family of reflection traveltime curves on depth sections has the form:

$z=h+k*\xi \text{,}$ (67)

where $\xi$ is the location of the result trace.

By substituting the parameter $h$ from Equation (67) into Equation (66), we obtain the family of functions that connect families (66) and (67), whose envelope is used to accomplish the required transformation. The equation of this family of functions is as follows:

$\tilde{T}=\frac{2}{\beta *V_0}*\text{ln}\frac{1+\beta *z+k*\beta *\left(\eta -\xi \right)+A}{1+\sqrt{1+k^2}}\text{,}$ (68)

where

$A=\sqrt{{\left(1+\beta *z+k*\beta *\left(\eta -\xi \right)\right)}^2+k^2}\text{.}$

$\begin{array}{l}{{\tilde{T}}^{\prime }}_k=\frac{\left(\left(1+\beta *z+k*\beta *\left(\eta -\xi \right)\right)*\beta *\left(\eta -\xi \right)+k+A\right)*\beta *\left(\eta -\xi \right)}{A}\\ -\frac{A+1+\beta *z+k*\beta *\left(\eta -\xi \right)}{\sqrt{1+k^2}*\left(1+\sqrt{1+k^2}\right)}=0.\end{array}$ (69)

This equation is solved by the numerical Newton method.

Figure 20 depicts the outcome of testing this procedure

The main advantage of VSP is the absence of surface waves, which makes it possible to obtain additional information in relation to ground-based observations. Therefore, improving the existing procedures for processing VSP data and developing new ones is important.

Galperin (1974) ( [45] ) proposed the VSP method. [32] [46] [47] and others have contributed to the development of this method.

First, it is necessary to find the equation of a family of traveltime curves for reflected from the inclined reflectors and recorded in the VSP data waves.

The equation of the family of traveltime curves for waves reflected from flat boundary and recorded on the VSP data were obtained from these works and has the next form

$T=\sqrt{{\left(2*T_0*V-h*\sqrt{4-t^2*V^2}+t*p*V\right)}^2+{\left(p*\sqrt{4-t^2*V^2}+h*t*V\right)}^2}/\left(2*V\right)\text{,}$

(70)

where $h$ is location of geophone in borehole and $p$ is the source location.

We assume that the velocity of the wave propagation is constant. As previously demonstrated, the equation for the family of traveltime curves for reflected from an inclined reflectors recorded on the time sections waves has the form (7).

First of all let’s transform registered on VSP data reflected waves into timesections ( Figure 21 ).

The equation for the family of traveltime curves for waves reflected from flat reflectors on times sections is as follows:

${\tilde{T}}_0=T_0+t*\xi \text{,}$ (71)

where $\xi$ is the location of the result time section trace. When we express the parameter $T_0$ through other terms of the expression (71)

$T_0={\tilde{T}}_0-t*\xi$ (72)

and substitute it into (71) will be obtain the equation of family of functions that connecting families (70) and (71).

$\tilde{T}=\sqrt{{\left(2*{\tilde{T}}_0*V-h*\sqrt{4-t^2*V^2}+t*V*\left(p-2*\xi \right)\right)}^2+{\left(p*\sqrt{4-t^2*V^2}-h*t*V\right)}^2}/\left(2*V\right)\text{.}$

(73)

The following equation was used to determine the envelope of the family (73):

$\begin{array}{l}{{\tilde{T}}^{\prime }}_t=\left(2*{\tilde{T}}_0*V-h*\sqrt{4-t^2*V^2}+t*V*\left(p-2*\xi \right)\right)\\ *\left(h*t*V+\left(p-2*\xi \right)*\sqrt{4-t^2*V^2}\right)\\ +\left(p*\sqrt{4-t^2*V^2}+h*t*V\right)*\left(h*\sqrt{4-t^2*V^2}-p*t*V\right)=0\text{.}\end{array}$ (74)

To solve this equation is used the numerical Newton’s method.

The method described in the introduction is used to store as a constant the difference between the second derivatives of the envelope and any function of the family (70) at the point of tangency.

The time migrated section may be obtained using Equations (70) and (8) ( Figure 22 ).

We present the Equation (8) again for completeness of presentation:

${\tilde{T}}_0=2*\left(T_0+t*\xi \right)/\sqrt{4-t^2*V^2}\text{.}$

Let’s represent the parameter $T_0$ through the Equation’s (75) other parameters and variables

$T_0={\tilde{T}}_0*\sqrt{4-t^2*V^2}/2-t*\xi .$ (75)

Substituting (75) into (70) we obtain the equation that connects these families:

$\tilde{T}=\sqrt{{\left[\sqrt{4-t^{2*}{V}^2}*\left({\tilde{T}}_0*V-h\right)+t*V*\left(p-2*\xi \right)\right]}^2+{\left(p*\sqrt{4-t^{2*}{V}^2}+h*t*V\right)}^2}/\left(2*V\right)\text{.}$

(76)

To find the envelope of this family, we must solve the following equation:

$\begin{array}{l}{{\tilde{T}}^{\prime }}_t=\left[\sqrt{4-t^2*V^2}*\left({\tilde{T}}_0*V-h\right)+t*V*\left(p-2*\xi \right)\right]\\ *\left[-t*V^2*\left({\tilde{T}}_0*V-h\right)/\sqrt{4-t^2*V^2}*V*\left(p-2*\xi \right)\right]\\ +\left(p*\sqrt{4-t^2*V^2}+h*t*V\right)*\left(-p*t*V^2/\sqrt{4-t^2*V^2}+h*V\right)=0\text{,}\end{array}$ (77)

$t=\sqrt{2*\left(\sqrt{4*{\gamma }^2+{\delta }^2}-\delta \right)/\left(\sqrt{4*{\gamma }^2+{\delta }^2}*V^2\right)}\text{,}$ (78)

where

$\gamma ={\tilde{T}}_0*V*\left(p-2*\xi \right)+2*z*\xi$ and $\delta ={\tilde{T}}_0*V*\left({\tilde{T}}_0*V-2*z\right)+4*\xi *\left(p-\xi \right)$,

and substitute the value $t$ in (76).

The approach described in the introduction is used to save as a constant the difference between the second derivatives of the envelope of the family (74) and any function of the family (68) at the point of their touching.

To obtain a depth migrated section, we can use the equations of families (21) and (70) ( Figure 23 ).

From Equation (21), let’s represent the parameter $T_0$ through the other parameters and variables of this equation:

$T_0=z*\sqrt{4-t^2*V^2}/V-t*\xi .$ (79)

We can obtain the equation of family of functions that connects these families by inserting (79) into (70):

$\tilde{T}=\sqrt{{\left[\left(2*z-h\right)*\sqrt{4-t^2*V^2}+t*V*\left(p-2\xi \right)\right]}^2+{\left(p*\sqrt{4-t^2*V^2}+h*t*V\right)}^2}/\left(2*V\right)\text{.}$

(80)

The envelope of this family can be found from condition ${{\tilde{T}}^{\prime }}_t=0$, that leads to:

$\begin{array}{l}{{\tilde{T}}^{\prime }}_t=\left[\left(2*z-h\right)*\sqrt{4-t^2*V^2}+t*V*\left(p-2\xi \right)\right]\\ *\left[\left(p-2\xi \right)-\left(2*z-h\right)*t*V/\sqrt{4-t^2*V^2}\right]\\ +\left(p*\sqrt{4-t^2*V^2}+2*t*V\right)*\left(h-p*t*V/\sqrt{4-t^2*V^2}\right)=0\end{array}$ (81)

and

$t=2*C/\sqrt{4*A^2+C^2}/V,$ (82)

where $A=\left(2*z-h\right)*\left(p-2*\xi \right)+p*z$; $C=\sqrt{4*A^2+B^2}-B$;

$B={\left(2*z-h\right)}^2+p^2-h^2$.

If $t$ is substituted into (80), we obtain the equation of the envelope of the family (80):

$\widetilde{\stackrel{˜}{T}}=\sqrt{4*\left[{\left(2*z-h\right)}^2+p^2\right]-4*C^2*B/D+16*A^2*C/D}/\left(2*V\right)\text{,}$ (83)

where $D=4*A^2+C^2$.

The VSP data will transform into depth migrated sections by summation along the envelope (83). The approximate approach described in the introduction is used to keep the “touch character” as a constant.

Figure 24 depicts the outcome of testing this procedure.

The geometry of the refraction path is detailed in [20] [48] and others.

The equation of the family of traveltime curves for waves refracted from flat boundary were obtained from these works and has the next form ( Figure 25 ):

$T=\frac{T_0*\sqrt{V_1^2-V_0^2}}{V_1}+\frac{x*\sqrt{4-t^2*V_0^2}}{2*V_1}+\frac{t*x*\sqrt{V_1^2-V_0^2}}{2*V_1}\text{,}$ (84)

where $V_0$ is the velocity of the propagation of reflection waves and $V_1$ is the velocity of the propagation of headwaves.

The refraction arrival point is equal to:

$L=\frac{2*T_0*V_0^2}{\sqrt{4-t^2*V_0^2}*\sqrt{V_1^2-V_0^2}-t*V_0^2}\text{.}$ (85)

To obtain the time section, we express the parameter $T_0$ through the remaining parameters and variables of Equation (7) and substitute it into Equation (84):

$\tilde{T}=\frac{\left({\tilde{T}}_0-t*\xi \right)*\sqrt{V_1^2-V_0^2}}{V_1}+\frac{x*\sqrt{4-t^2*V_0^2}}{2*V_1}+\frac{x*t*\sqrt{V_1^2-V_0^2}}{2*V_1}\text{.}$ (86)

To determine the envelope of Equation (86), condition (5) is used:

$t=\sqrt{\frac{4*{\left(x-2*\xi \right)}^2*\left(V_1^2-V_0^2\right)}{V_0^2*\left[x^2*V_0^2+{\left(x-2*\xi \right)}^2*\left(V_1^2-V_0^2\right)\right]}}$ (87)

and this value is substituted into (86).

To preserve “touch character” as a constant, the approximation approach given in the introduction is utilized.

To obtain the time migrated section, we used Equation (8).

The parameter $T_0$ is expressed through the rest of the parameters and variables of this equation, and the result is plugged into (84):

$\begin{array}{l}\tilde{T}=\frac{\left({\tilde{T}}_0*\sqrt{4-t^2*V_0^2}/2-t*\xi \right)*\sqrt{V_1^2-V_0^2}}{V_1}\\ +\frac{x*\sqrt{4-t^2*V_0^2}}{2*V_1}+\frac{x*t*\sqrt{V_1^2-V_0^2}}{2*V_1}\text{.}\end{array}$ (88)

From condition (5), the value of $t$ is found:

$t=\sqrt{\frac{4*{\left(x-2*\xi \right)}^2*\left(V_1^2-V_0^2\right)}{V_0^2*\left[{\left({\tilde{T}}_0*\sqrt{V_1^2-V_0^2}-x\right)}^2*V_0^2+{\left(x-2*\xi \right)}^2*\left(V_1^2-V_0^2\right)\right]}}\text{.}$ (89)

By substituting the values of $t$ into (88), we find the envelope of the family (88) used to obtain the time migrated section.

To transform refracted waves into depth migrated sections, we may use the equations of families (109) and (21).

Let’s represent the parameter $T_0$ from Equation (21) through the other parameters and variables of this equation:

$T_0=z*\sqrt{4-t^2*V^2}/V-t*\xi .$ (90)

After substituting (90) into (84), we obtain the equation of the family of functions that connect these families.

$\begin{array}{l}\tilde{T}=z*\sqrt{4-t^2*V_0^2}*\sqrt{V_1^2-V_0^2}/\left(V_0*V_1\right)-t*\xi *\sqrt{V_1^2-V_0^2}/V_1\\ +x*\sqrt{4-t^2*V_0^2}/\left(2*V_1\right)+t*\left(x-2*\xi \right)*\sqrt{V_1^2-V_0^2}/\left(2*V_1\right)\text{.}\end{array}$ (91)

From condition (5) let’s find the parameter $t$:

$t=\frac{2*\left(x-2*\xi \right)*\sqrt{V_1{}^2-V_0^2}}{V_0*\sqrt{{\left(2*z*\sqrt{V_1{}^2-V_0^2}+V_0*x\right)}^2+{\left(x-2*\xi \right)}^2*\left(V_1{}^2-V_0^2\right)}}$ (92)

and substitute it into (91).

Obtained in this way the envelope of the family (91) will be used for transforming recorded refracted waves into depth migrated sections.

The approximate method described in the introduction is used to save “touch character” as a constant.

Figure 26 depicts the outcome of testing these procedures.

As an innovative method of transforming seismic data for further analysis, it is proposed to obtain so-called “floating” time sections. These sections differ from CMP time sections primarily in that the velocities used for their realization have a significant effect primarily on the reflection timecurve of the recorded waves, rather than on the reflection amplitude. If the velocity corresponds to the actual velocity, the position of the reflecting horizon on the “floating” time section is correct. The greater the difference between the accepted and true velocities, the greater the difference between the true and accepted recorded wave times. In this case, the dynamics of the reflecting horizons on the “floating” time sections is practically independent of the velocities used. The value of the difference between the recorded reflections from the same horizon on the CMP sections and “floating” time sections can be used to calculate the real propagation velocity of the reflected wave.

The procedure for obtaining a “floating” time section is divided into two stages. First, the CMP seismograms are transformed into “floating” seismograms, in which the family of the reflection traveltime curves is described by a function depending on parameters such as the estimated $\tilde{V}$ and true wave propagation velocities $V$. At the next stage, parameter $V$ will be excluded from the process of determining the envelope used to finally obtain a “floating” time section.

As previously proven, (43) is an equation of the family of traveltimes curves for reflected from the flat reflector waves for the average velocity model recorded on the CMP gathers.

Let’s build a procedure for transformation CMP seismograms into new seismograms, the traveltime curves of reflected waves on which are described by the following equation:

${\tilde{T}}_0=T_0+\left(t+2/V-2/\tilde{V}\right)*\xi +\rho *\left[1+{\left(t+2/V-2/\tilde{V}\right)}^2*V^2\right]/\tilde{V}\text{.}$ (93)

This is the family of functions of two variables $\left(\xi ,\rho \right)$ that depends on three parameters $\left(t,V,\tilde{V}\right)$, where $\tilde{V}$ is the assumed average velocity, $\xi$ is the location of the result seismogram, and $\rho$ is the location of the result trace on this seismogram ( Figure 27 ).

To find an equation for the family of functions that connect families (51) and (93), we represent parameter $T_0$ by other parameters and variables from Equation (93) and substitute it into (51):

$T=\sqrt{{\left\{{\tilde{T}}_0-\left(t+\frac{2}{V}-\frac{2}{\tilde{V}}\right)*\xi -\frac{\rho *\left[1+{\left(t+2/V-\frac{2}{\tilde{V}}\right)}^2*{\tilde{V}}^2\right]}{\tilde{V}}+t*\eta \right\}}^2+\frac{\left(4-t^2*V^2\right)*x^2}{{\tilde{V}}^2}}\text{.}$ (94)

To obtain the envelope of the family (94), two parameters should be excluded using the following conditions:

${{\tilde{T}}^{\prime }}_t=0\text{and}{{\tilde{T}}^{\prime }}_V=0\text{.}$ (95)

$\begin{array}{l}{{\tilde{T}}^{\prime }}_t=\left\{{\tilde{T}}_0-\left(t+2/V-2/\tilde{V}\right)*\xi -\rho *\left[1+{\left(t+2/V-2/\tilde{V}\right)}^2*{\tilde{V}}^2\right]/\tilde{V}+t*\eta \right\}\\ *\left[-\xi -2*\rho *\left(t+2/V-2/\tilde{V}\right)+\eta \right]-t*x^2=0\text{.}\end{array}$ (96)

$\begin{array}{l}{{\tilde{T}}^{\prime }}_V=\left\{{\tilde{T}}_0-\left(t+2/V-2/\tilde{V}\right)*\xi -\rho *\left[1+\left(t+2/V-2/\tilde{V}\right)*{\tilde{V}}^2\right]/\tilde{V}+t*\eta \right\}\\ *\left[\xi +2*\rho *\left(t+2/V-2/\tilde{V}\right)*\tilde{V}\right]*V-2*x^2=0\text{,}\end{array}$ (97)

From Equations (96) and (97), we find values $t$ and $V$:

$t=\left(D-2\right)/V\text{and}V=2*\rho *D^2*\tilde{V}/\left[\left(4*\rho -\xi \right)*D+2*\eta \right]\text{,}$ (98)

where

$D=\left(2+t*V\right)=\frac{-2*B+\sqrt{4*B^2-A*C}}{A},$

where

$A=4*\left({\tilde{T}}_0-\rho /\tilde{V}\right)*\tau *\rho *\tilde{V}+\eta *{\xi }^2+2*\left({\eta }^2-x^2\right)*\left(4*\rho -\xi \right)\text{,}$

$B=\left[{\eta }^2-x^2-\eta *\left(4*\rho -\xi \right)\right]*\tau \text{,}$

$A*{\left(2+t*V\right)}^2+4*B*\left(2+t*V\right)+C=0\text{,}$