Hyperbolic Velocity Model

Asymptotically bounded velocity profiles describe the vertical velocity variations in compacted sediments in a more realistic way than unbounded velocity models, and allow presenting the subsurface by a smaller number of thicker layers. The first and the simplest asymptotically bounded model is the Hyperbolic velocity profile proposed by Muscat in 1937, and our paper is an extension of this early study. The Hyperbolic model has an advantage over other bounded models: The velocity increases with depth and approaches the limiting value with a more smooth and gradual rate. We derive the time-depth relationships, forward and backward transforms between the instantaneous velocity profile and the effective models (average, RMS and fourth order average velocities), study the trajectories for pre-critical and post-critical curved rays and derive the equations for traveltime, lateral propagation and arc length. We compare the ray paths obtained with the Hyperbolic model and with the other bounded velocity profiles.


Introduction
The Hyperbolic velocity model was first proposed by Muscat [1] and published in 1937.However, since then the model was not extensively studied and is unjustifiably ignored in the literature.The objective of this research is to extend the original study and to correct the inaccuracies.We show the place of the Hyperbolic model among the other asymptotically bounded models, analyze its basic relationships and attempt to develop a complete theory.
Asymptotically bounded velocity models describe the velocity profile in compacted sediments, where the velocity gradually increases with depth and eventually approaches a limiting value.These models make it possible to describe a vertical velocity profile with a smaller number of intervals as compared to the classical unbounded models, such as linear velocity vs. depth [2,3], unbounded exponent [4], linear slowness [5], "sloth" (linear variation of slowness squared), e.g., [6], parabolic model [7,8], Faust velocity model [9,10] with a reference depth and different root indices.The unbounded models are described by two parameters: the instantaneous velocity at the top interface a V and the vertical velocity gradient a at the same level.The Faust model includes also the root index n, normally .Asymptotically bounded models require an additional parameter: the limiting value of velocity at infinite depth.Two models of this family were studied by Ravve and Koren: the Exponential asymptotically bounded model [11,12] and the Conic model [13].The asymptotically bounded profiles can be used, in particular, as velocity trend functions for the constrained velocity inversion with the best (e.g., least-squares) fit of the input data [14].Examples of asymptotically bounded models are presented below.For each model, we first give the original formulation of the velocity profile as it appears in the original works by the authors, and then we convert it to a canonical form in terms of the "standard" parameters and  The Hyperbolic velocity model by Muscat [1], In our notation, the Hyperbolic profile reads  The Exponential velocity model by Muscat [1], where the parameters are 2 2 2 2 arccosh , .
A detailed review on unbounded and bounded velocity models is given by Kaufman [16].Figure 1 shows graphs of the instantaneous velocity vs. depth for the five asymptotically bounded velocity models mentioned above.
For all models, we assume the same velocity profile parameters: . The vertical gradients of the velocity vs. depth are plotted in Figure 2. It is interesting to note that among the five models presented, the Muscat Hyperbolic model (Equation (1) and grey line on the plot) approaches the limiting value in the slowest and the most gradual manner.The "second slow" is the Conic velocity model (red line), and the "third slow" is the EAB model (blue line).An asymptotically bounded model can be characterized by its gradient-velocity relationship, which is actually the governing differential equation of the velocity model.
This paper is structured as follows.We define the Hyperbolic model using 1) the original Muscat [1] formulation-depth vs. velocity, 2) the physical parameters: maximum gradient R, length scale and vertical shift h, and 3) the "technical" or geophysical parameters: top interface velocity a V , top gradient a and asymptotic velocity .We introduce the dimensionless asymptotic factor Q k V  M that simplifies the transform equations.First we derive the time-depth and the depth-time relationships.
Next we proceed to forward transforms from the instantaneous velocities to the effective models, such as the average, the RMS and the fourth order average velocity.
Then we study the inversion problems, considering the inversion with the instantaneous velocities and gradients, and the inversion with the effective models, i.e., the average or the RMS velocities given vs. time or depth.Next we comment on the two types of curved rays existing in all asymptotically bounded models, depending on the initial take-off angle, and derive the trajectories of the ray paths for the Muscat velocity profile.For both types of the curved rays we derive the lateral propagation, the traveltime and the arc length.

The Hyperbolic Velocity Profile
Muscat [1] defined the Hyperbolic model by where V is the instantaneous velocity, z is depth measured from the top interface, a is the top interface velocity, A is the characteristic distance (scale) that affects the top gradient a , and is the asymptotic velocity.Inverting Equation (1), we obtain The velocity gradient becomes where At , the top gradient is Introduce Equation ( 14) into Equation (13).In our notation, the Hyperbolic profile reads We call values and V  the technical parameters of the profile.At a definite height above the earth surface (above the upper interface), where , a a V k z  , the instantaneous velocity vanishes.According to Equation ( 15), .
Introduce the absolute frame , where the instantaneous velocity vanishes at the origin is the shift between the two frames of reference.In the absolute frame, the velocity profile simplifies to where 2 and .
We call values and the physical parameters of the profile.Note that the linear velocity profile, where the ray trajectories are circular arcs, is a particular case of the Hyperbolic model with and V  , in such way that their product remains a finite value, and parameter becomes the constant velocity gradient of the linear model.The velocity gradient of the Hyperbolic model reads At the absolute origin 0 z   k , the velocity gradient reaches its maximum value max .Comparing Equations (17) and (19), we conclude that Equation ( 20) is the governing differential equation of the Hyperbolic velocity profile.It can be used to plot the gradient-velocity diagram.Such diagrams for several asymptotically bounded velocity models are studied in Appendix A.
Introduce the normalized (dimensionless) velocity v , the normalized gradient and the normalized absolute depth , Note that parameter Q is the reciprocal characteristic length.With these notations, the Hyperbolic velocity profile simplifies to The technical parameters of the velocity profile are related to the physical parameters, 1 The inverse relationship is

Asymptotic Factor
To simplify the equations for velocity transforms, it is suitable to introduce a special parameter M.This parameter can be defined at any point of the profile, and in particular, at the top and the bottom interfaces of an interval, where z  is the interval thickness (the vertical distance between the two interfaces), subscript a is related to the top interface 0 z  , and subscript b is related to the bot- where a and b V are the top and bottom instantaneous velocities, respectively.Next, it follows from Equation (26) that parameter M is the inverse normalized measure of the difference between the velocity at the given depth level and the asymptotic velocity .Equation (26) can be inverted, The velocity gradient is also related to the asymptotic factor, 2 , It follows from Equation (25) that and We use Equations ( 28) and (29) to get the interface gradients, a and , through the increment of the asymptotic factor, , Equations ( 27) and (29) result in the average gradient on the interval, ave , expressed either through the interface asymptotic factors Introduction of Equation (28) into Equation (33) leads to ave .
The average gradient on the interval with the Hyperbolic velocity profile is the geometric average of the top and bottom interface gradients.
Given the velocity and its gradient at one interface, one can calculate these parameters at the other interface.The calculations can be done either in depth or in time.Four problems of this kind are considered in Appendix C.

Depth-Traveltime Relationship
Integrate the slowness to get the vertical traveltime vs. the interval thickness, The traveltime equation can be written in terms of asymptotic factors at the top and bottom interfaces,  25) and (29), we obtain where the top asymptotic factor To get the vertical distance vs. traveltime we should invert Equation (26 Equation (38) should be solved for the unknown bottom asymptotic factor b M , Taking exponent from both sides of Equation (39), we get Equation ( 40) can be solved with the Lambert function, After the bottom asymptotic factor b M or the bottom interface velocity b is found, the interval thickness can be established with Equation (29),

Hyperbolic and Non-Hyperbolic Moveout
In the absence of the intrinsic anellipticity, the hyperbolic parameter W and the non-hyperbolic parameter H on the interval are defined as Introduce the velocity profile from Equation (17).The hyperbolic parameter W becomes The non-hyperbolic parameter H becomes With the use of the top and bottom asymptotic factors, the hyperbolic parameter becomes ln Introducing Equation (37) for the traveltime into Equation (47), we obtain the local RMS velocity U over the interval.By definition, U W  t , so The non-hyperbolic parameter becomes With the use of Equation ( 29), the non-hyperbolic parameter simplifies to When the parameters of the velocity profile are specified, the top asymptotic factor M a is a known value.The bottom asymptotic factor M b can be presented either vs. depth (interval thickness) or vs. traveltime.Thus, the hyperbolic and non-hyperbolic parameters become functions of depth or traveltime, accordingly.The anellipticity induced by the vertically varying velocity is defined as the fractional difference between the fourth-order average velocity and the RMS velocity , Parameter  can be also considered as a function of depth or vertical time.For a particular case of a single infinite layer (half-space) with any vertical velocity profile, The graph for the induced anellipticity  is plotted vs. depth in Figure 3 for three asymptotically bounded velocity models: Exponential, Conic and Hyperbolic.For all the three models, the parameters of the velocity profile are: 3 km/s, a V  1 and .At the surface, the anellipticity is zero as there are yet no accumulated variations of the instantaneous velocity.The induced anellipticity is always positive.It reaches a maximum value a definite depth and then vanishes at the infinity, where the medium velocity is asymptotically constant.

Forward Dix Transform
Consider a package of n layers (vertical intervals), where the nodes (interfaces) are enumerated from zero, and layers are enumerated from 1. Interval n connects nodes (top interface) and n (bottom interface).The nodal average velocity 1,n V , RMS velocity and fourthorder average velocity are .
where n is the one-way interval traveltime, n t  z  is the layer thickness, W n and H n are the interval hyperbolic and non-hyperbolic parameters, respectively.For 1 n  we set 0 in Equation (53).The effective velocities (average, RMS and fourth-order average) can be also defined for any internal point of the interval.0 t 

Inverse Dix Transform
Recall that the Hyperbolic velocity profile on the interval is defined by the three parameters: the top interface instantaneous velocity a , the top interface gradient a , and the asymptotic velocity V  .We consider that the asymptotic velocity is always given a priori.When the two other parameters, a and a k , are also known, then velocity transforms are considered forward.When one or both parameters are unknown (with another data specified instead), we deal with the velocity inversion.There are three groups of inverse transforms studied in Appendices D, E and F.

V V k
Appendix D considers the inversion that does not involve the RMS velocity.These formulations deal with the instantaneous velocity and its gradient only.We solve a problem where the two velocities are given at the interfaces, a and b V , or-alternatively-the two gradients, a and b .Another kind of problem is when the velocity and its gradient are given at the different interfaces of the interval, i.e. the velocity is given at the top interface and the gradient-at the bottom interface, a and b , or vice versa, a k and b V .We solve also a problem where the instantaneous velocity is given at the bottom interface and at the intermediate point of the interval, b and c V .These problems are studied both vs. depth and vs. time.
Appendix E considers the inversion with the RMS ve-locity specified at the interfaces vs. depth or time with a single parameter unknown, either a V or a .We consider also a problem with the traveltime specified vs. the interval thickness, also with a single parameter unknown.Finally, we consider the RMS velocity specified vs. both depth and time, with the two parameters unknown, and .
In Appendix F we study the two-interval inversion.The RMS velocity is given vs. depth or time at the two interfaces and at an internal point of the interval.Alternatively, depth can be specified vs. traveltime at the three points.Both parameters of the velocity profile are unknown.This is a so-called three-point or two-interval inversion.k

Ray Trajectories
In this section we establish the trajectories of non-vertical rays.Due to Snell's law, in 1D medium the horizontal slowness is constant, and the ray angle p  (measured from the vertical axis) becomes Introduce the ray parameter m, where Q and R are the physical parameters of the Hyperbolic velocity profile, is the normalized ray slowness, and m is its inverse value.We call parameter m "eccentricity of the ray trajectory" as it is very similar to the eccentricity of the hyperbolic and elliptic rays of the Conic velocity model [13].With Equation ( 17), the sine of the ray angle becomes where (the conjugate eccentricity squared) may be positive or negative.Introduce the dimensionless coordinates, The tangent of the ray angle becomes Integrating Equation (59), we obtain Copyright © 2013 SciRes.
where c x is the constant of integration.This integral can be reduced to To obtain the integral on the right side of Equation (61), we consider two cases, or two ranges of the eccentricity: m > 1 (pre-critical rays) and m < 1 (post-critical rays), We emphasize that two kinds of rays exist for any monotonously increasing and asymptotically bounded velocity model, and in particular, for the Hyperbolic model.The pre-critical rays that may start on the earth surface, propagate to the infinite depth, and their curvature asymptotically vanishes.The post-critical rays have a limited propagation depth.Their arc-like trajectories have a finite minimum curvature at the turning point, and these rays return to the earth surface.Note that at any point of the trajectory, the ray path curvature  depends on the velocity gradient k only, In particular, the linear velocity model with a constant velocity gradient leads to ray trajectories of constant curvatures, i.e. to the circular arcs.
The critical rays with the unit eccentricity 1 m  are the limit case between the two types of rays.Their takeoff angle (the ray angle at the upper interface) is called the critical angle C  .It follows from Equations ( 54) and (55) that the critical take-off angle is arcsin .
It follows from Equations ( 61) and ( 62) that the trajectories of the pre-critical and post-critical rays are At infinite depth, pre-critical rays become asymptotically straight.Equation (57) leads to However, although the slope of these rays converges to a constant value and their curvature becomes infinitesimal, the pre-critical rays of the Hyperbolic model have no asymptotic straight line, unlike the pre-critical rays of the EAB and the Conic models.The pre-critical, critical and post-critical rays are plotted in Figure 4 for the Hyperbolic, the Conic and the Exponential (EAB) models.Parameters of the velocity profile are the same as above.The three columns of numbers to the right of the plot area are velocities for the three models at the specified depth levels.

Maximum Penetration Depth
Pre-critical rays penetrate to infinite depth.The maximum penetration depth of post-critical rays follows from Equation (59).At the turning point, the ray angle π 2   , and thus, its tangent becomes infinite.This leads to a quadratic equation with a single positive root Recall that is the dimensionless depth measured from the absolute origin (above the upper interface).The maximum penetration depth in units of length, measured from the upper interface reads ẑ where a  and b  are ray angles at the departure and the destination points, respectively.Equations ( 17) and ( 19) make it possible to eliminate depth and to express the gradient through the velocity, Next we apply Snell's law and obtain the vertical gradient vs. the ray angle, The indefinite integral on the right side of this equation essentially depends on the range of the eccentricity m, resulting in For the critical ray, To simplify the notations, we introduce one more function of the ray angles at the endpoints, The normalized lateral propagation becomes The normalized traveltime is The normalized arc length is (84) , and the departure angle is critical, The current depth can be also expressed through the ray angle.It follows from Equation (17) that and therefore In Figure 5 we plot the graphs for the traveltime vs. ray path arc length for the pre-critical, the critical and the post-critical rays of the Hyperbolic velocity profile.
The "trigonometric" solution for the lateral propagation and traveltime of the post-critical rays was obtained (in a different form) by Muscat (1937).However, it was not pointed out in this early study that the solution was related to the post-critical rays only, and that the other, "hyperbolic" solution exists for the pre-critical rays (and a "transient" solution for the critical rays, which are the limit case between the two basic types of rays).
Note that for the vanishing or infinitesimal parameter , the shape of the trajectory, the lateral propagation, the traveltime and the arc length of the Hyperbolic model ray path converge to the corresponding character- istics of the linear velocity profile.In this case the asymptotic velocity V  becomes unbounded, so that the product remains a finite value and converges to a constant gradient of the linear velocity model.The eccentricity m becomes infinitesimal, and Functions I and J from Equations ( 79) and (81) simplify to Equation (66) comes to This is an equation for the circular arc of radius 1 , pR whose center is located at where is the ray slowness.The lateral propagation, Equation (82), becomes

Full Arc of Post-Critical Ray
Consider two points on the earth surface, the transmitter and the receiver, located x  distance apart.The goal is to trace the full arc of the post critical turning ray that connects the two points.Note that due to the symmetry of the arc, the ray angle at the destination point b  is related to the take-off angle a Applying Equations ( 79), ( 81) and (82), we obtain Equation ( 98) should be solved numerically for the unknown eccentricity m.To obtain the initial guess, we assume that the distance x  is small.Then the take-off angle a  approaches π 2 , and according to Equation (97), the eccentricity exceeds the sine of the critical angle only slightly.We assume sin , where m  is a small positive value.Next we expand Equation (98) into the Taylor series and neglect the high order terms, The cubic Equation (100) has a single positive root.For example, for the velocity profile 3 km/s a V  , .The take-off angle becomes a .The arcs are plotted in Figure 6 for the three asymptotically bounded velocity models.In the shallow region, the Hyperbolic model has a smaller vertical gradient (and thus, a smaller curvature) than the Conic and the EAB models, and thus, the Hyperbolic model yields a smaller take-off angle.The ray path arc of the Hyperbolic model passes above the Conic and the EAB arcs.0.67638 7.67

Boundary Value Ray Tracing
Given data are the departure point  , p are positive.Assume also b a , which is also not a limitation (one can reverse the endpoints otherwise).Since the ray tracing equations depend on the type of ray, we need to determine, whether the ray is pre-critical or post-critical.For this, one can plot a critical path that starts at the departure point at the critical take-off angle If the destination point lays to the left from the critical trajectory, then the ray path is pre-critical.The ray path is post-critical if the destination point lays to the right.The critical lateral propagation x  is delivered by Equation ( 63), which can be rearranged as Given the vertical coordinates of the source and the receiver, a and b , we calculate the critical lateral propagation, and then apply the criterion z z pre-critical ray, critical ray, post-critical ray.
The velocities at the end points of the trajectory and are known values.It follows from Snell's law that the ray angles at the end points of the trajectory are the functions of the eccentricity alone, Note that for the pre-critical rays and for the post critical rays before the turning point, the ray angle is acute, while for the turning rays after the turning point the ray angle is obtuse, arcsin , before turning point, π arcsin , after turning point.
Equation ( 82) relates the lateral propagation x  to the ray angles at the endpoints, which, in turn, depend on the eccentricity according to Equations ( 103) and (104), where

Conclusion
The Hyperbolic asymptotically bounded exponential velocity model has been studied and compared to other asymptotically bounded models, in particular, the Exponential and the Conic.The forward and the inverse velocity transforms are derived.The Hyperbolic model allows a better representation of the vertical velocity variations in compacted sediments, especially in the case of thick layers.An advantage of the Hyperbolic model is that the instantaneous velocity reaches the asymptotic value in a more slow and gradual fashion, as compared to other asymptotically bounded models.Ray tracing equations have been derived.The ray trajectories, traveltimes and arc lengths have been studied analytically, and the boundary value ray tracing problem have been solved.We have tried to present a complete theory for both vertical and non-vertical rays propagating through the Hyperbolic model.Application of the Hyperbolic velocity distribution enables us to present realistic geological models using fewer parameters, as compared to the classical linear velocity function.We showed that the linear velocity function is a limiting particular case of the Hyperbolic model.

Appendix A. Gradient-Velocity Diagrams
In this appendix we derive the gradient-velocity diagrams for the five asymptotically bounded velocity models.In all cases we pass to a shifted frame , in which the governing equations are essentially simplified.The value of the vertical shift h is different for all models.For all models the "absolute" origin corresponds to a point of maximum gradient.For all models except the Exponential slowness, this is also the point of a vanishing instantaneous velocity.As we show below, for the Exponential slowness model, the absolute origin corresponds to the half-limiting velocity

A.1. The Hyperbolic Muscat Model
The velocity profile is given by Establish the depth level where the velocity vanishes.This point is located above the earth surface, , Introduce the shifted frame, .The velocity profile becomes The vertical gradient is At the absolute origin , the gradient is maximal, Introduction of Equation (A-5) into (A-3) and (A-4) leads to Finally, elimination of depth from the two equations of Equation Set (A-6) results in

A.2. The Exponential Muscat Model
The velocity model is described by The vertical shift is , and in the shifted frame, The velocity gradient is Eliminate the absolute depth from Equations (A-9) and (A-11) and obtain the gradient-velocity relationship,

A.3. The Exponential Slowness Model
The profile equation reads Unlike the other asymptotically bounded models mentioned in the introduction, the velocity in the Exponential slowness model does not vanish at a finite negative depth.The velocity vanishes at and approaches to the asymptotic value V z    at .The gradient of the velocity is The gradient   k z vanishes at both remote ends, , and z     k z has a single critical point: the maximum of the gradient occurs at , We emphasize that the logarithm in Equation (A-15) may prove to be both positive and negative.At the point z z   , the maximum gradient and the velocity are max , 4 2 Note that in case when the depth of the maximum gradient is positive, , this point really exists underground, and the gradient first increases, then accepts the maximum value Below this point, the gradient begins to decay, and eventually vanishes at the infinite depth.In case when the depth of the maximum gradient is negative, this point is above the earth surface, and throughout the whole depth range 0 the velocity gradient Next we assume the shift , and pass to the shifted frame, .The gradient accepts now a maximum value at the origin.Rearrange Equation (A-13), Note that the velocity profile in Equation (A-19) can be set in an alternative way, Finally, we eliminate depth from Equation (A-22) and obtain max 1 4

A.4. The Exponential Asymptotically Bounded Velocity Model (EAB)
In this case the velocity profile is The velocity vanishes above the earth surface at , z  ln .
In the shifted frame , the velocity profile simplifies to Finally, we eliminate depth from Equation (A-27) and obtain the governing equation, max 1.
Note that only for the EAB velocity model the diagram Equation (A-28) is linear.

A.5. Conic Velocity Model
For the Conic profile, the velocity and its gradient in the absolute frame are given by This equation can be rearranged as Next we eliminate the absolute depth from Equation (A-29) and get the governing differential equation of the Conic velocity model, The Conic velocity profile, Equation (A-30), can be also set in an equivalent form, through a hyperbolic and an inverse hyperbolic function, The gradient-velocity diagrams for the five asymptotically bounded models are plotted in Figure 7.Note that in the original frame of reference the asymptotically bounded models are described by the three parameters.As we mentioned, in the shifted frame where the velocity vanishes at the origin (or the vertical gradient accepts a maximum value at the origin), only two parameters are needed.These two parameters may be the maximum gradient max and the asymptotic velocity V k  as in Equation Set (A-33).The constant value that appears upon the integration of each equation is not a new parameter as it should be adjusted to match the maximum values and V .max  We emphasize that only for the EAB model the gradient-velocity relationship is linear: Derivative of an exponent is proportional to the same exponent.For all models except the Exponential slowness [8], the gradient decreases with the increase of velocity (and depth).In case of the Exponential slowness, the vertical gradient increases along with the velocity, until the velocity reaches one half of the asymptotic value 2 V  t this point the gradient reaches its maximum value max k then begins to decrease with depth.The point of maximum gradient in the Exponential slowness model may really exist in the subsurface, or it may be an imaginary point located above the earth surface (or above the upper interface of a layer).This point is real in case when the initial velocity a V does not exceed the halflimiting value, .A and . Furthermore, we comment on the special central symmetry between the two Muscat (1937) models: Hyperbolic (HM) and Exponential (EM), see Equation (A-33) and the two corresponding plots in Figure 7, where HM F and EM F are the corresponding gradientvelocity functions in Equation (A-33) for these two models.Although these two models are described by essentially different vertical velocity profiles, there is an apparent similarity in the gradient-velocity diagrams.

Appendix B. Lambert Function
The Lambert function [17]    for branch minus one.The value range is     for branch zero and 1 y     for branch minus one.In particular, this means that for a positive argument x only branch zero exists, while for a negative argument both branches do exist.Therefore, in the latter case, the branch index should be specified to avoid ambiguity.The derivative of the Lambert function is and for the infinitesimal argument Comment.A general comment is related to Appendices C to F. The transform equations are formulated in the dimensionless form, with the unknown top and bottom asymptotic factors, a M and b M .After the transform equation or equation set is resolved, we apply Equation ( 27) to find the top and bottom instantaneous velocities, a and b .If the transform is formulated in depth (i.e., the interval thickness is specified), we apply Equation (31) to find the top and bottom gradients of velocity, a and b k .If the transform is formulated in time (i.e., the interval traveltime is specified), then we first apply Equation (37) to establish the interval thickness , and then Equation (31) to find the top and bottom gradients.

Appendix C. Swapping Interfaces
In this appendix, we find the instantaneous velocity and its gradient at the bottom interface, given these parameters at the top interface, and vice versa, and consider these problems both vs. depth and vs. time.
Problem C1.Given the velocity a and its gradient a at the top interface and the layer thickness Problem C2.The instantaneous velocity and its gradient are specified at the bottom interface and the interval thickness is given.Velocity and gradient should be found at the top interface.Thus, parameters are given, and parameters are to be found.In this case calculate the bottom asymptotic factor V k M with the second equation of Equation Set (26) and apply the second equation of Equation Set (32) to get the top get a M .Problem C3.The velocity and its gradient at the top interface, a and a are given, the interval traveltime is known, and the bottom interface parameters b V and b should be found.Combining the first equation of Equation Set (32) and Equation (36), we eliminate the interval thickness and obtain where the top asymptotic factor a M is known from Equation (26).Equation (C-1) should be solved for the unknown bottom asymptotic factor b M .We use a cubic approximation for Equation (C-1) to get the initial guess, assuming the increment of the asymptotic factor Problem C4.The velocity and its gradient are specified at the bottom interface, along with the interval traveltime, and the profile parameters should be found at the top interface.Thus, parameters and , b b V k t  are given, while parameters a V and a are to be established.For this, we combine the second equation of Equation Set (32) and Equation (36), The bottom asymptotic factor b M is known from Equation (26), and Equation (C-3) should be solved for the unknown top asymptotic factor a M .The initial guess for

Appendix D. Inversion with Instantaneous Velocity
In this appendix we consider inversion problems that do not involve the effective models (average and RMS velocity).In this inversion group, one or both parameters at the top interface, a V and a are unknown, with the other data given instead.The group includes four problems vs. depth (interval thickness) and four similar problems vs. interval traveltime.
Problem D4.Gradients vs. time.Given data are the top and bottom interface vertical gradients, a and b , the asymptotic velocity  and the interval traveltime Find the top interface instantaneous velocity a .Solution.Introduce solution (D-1) into the traveltime Equa- tion (36).This leads to a nonlinear equation vs. the unknown interval thickness where .
Equation (D-2) should be solved numerically.It is suitable to normalize the gradients and the interval thick- The normalized equation becomes To get an initial guess, we expand Equation (D-4) into a power series and obtain a cubic approximation, After the interval thickness is found, apply solution (D-1).
Problem D5.Velocity and gradient vs. depth at different interfaces.Given data are the top interface velocity a , the bottom interface vertical gradient b , the asymptotic velocity V  and the interval thickness  , one can establish the single physical root of the quadratic equation, to get a cubic approximation, Problem D8.Given data are the top interface gradient a , the bottom interface instantaneous velocity b , the asymptotic velocity k V V  and the interval traveltime t  .Find the top interface velocity a V .Solution.Apply Equation (26) to calculate the bottom asymptotic factor b M .To calculate the unknown top asymptotic factor, we solve the nonlinear Equation (C-1).To get the initial guess, assume Introduce parameter , Equation (D-15) becomes

Appendix E. Two-Point Effective Model Inversion
In this appendix we consider the two-point (single-interval) inversion where the RMS velocity is specified at the interfaces of an interval vs. depth or traveltime, or alternatively, depth is specified vs. traveltime instead of the RMS velocity.One of the two parameters at the top interface-either the top velocity a , or the top gradient a -is a known value, while the other one is unknown and should be established.In all cases, the asymptotic velocity V k V  and the top interface absolute traveltime a are assumed known values.We consider also a special case when the RMS velocity is specified vs. both depth and time, and both parameters, and , are unknown.
Recall that a and b are one-way absolute top and bottom interface traveltimes (measured from the earth surface), t t t  is the one-way interval traveltime, and W is the hyperbolic parameter through the interval.Equation (E-1) can be arranged as .
Introduction of Equation (37) for the traveltime t  and Equation (47) for the hyperbolic parameter W into Equation (E-2) results in where A and B are known dimensionless parameters, The top asymptotic factor a M is delivered by Equation (26), and the nonlinear Equation (E-4) should be solved for the unknown bottom asymptotic factor b M .To get the initial guess, we assume a small increment of the asymptotic factor on the interval, , and expand Equation (E-3) into a power series.The cubic approximation reads where Equation (E-3) should be solved for We solve Equation Set (E-7) for the unknown asymptotic factors a M and b M .To obtain the initial guess for a M , assume that the increment of the asymptotic factor M  is small, and linearize Equation Set (E-7).This leads to Problem E3.RMS vs. time with unknown gradient.Given data are the RMS velocities at the top and bottom interfaces, 2,a and 2,b V , the interval traveltime t  and the top interface velocity a .Find the top gradient a .Solution.First we apply Equation (E-1) and calculate the hyperbolic parameter W. At this time, W is a known value.Next we apply Equation (48), where C is a known dimensionless parameter, the normalized local RMS velocity, and U is the non-normalized local RMS velocity on the interval.The top asymptotic factor a M is delivered by Equation (26), and the bottom asymptotic factor b M is established from Equation (E-9).To solve this nonlinear equation, an initial guess is needed.Assume that the increment of the asymptotic factor M  on the interval is small, and expand Equation (E-9) into a power series.The cubic approximation reads   .
Problem E4.RMS vs. time with unknown velocity.Given data are the RMS velocities at the top and bottom interfaces, 2,a and 2,b V , the interval traveltime V t  and the top interface gradient a .Find the top interface velocity a .Solution.We solve a set consisting of two equations: the first is Equation (E-9) combined with Equation (C-1), and the second is To get the initial guess, we linearize Equation Set (E-12) for a small increment of the asymptotic factor M  and obtain where D is a known dimensionless parameter-the normalized interval velocity, To get the initial guess for Equation (E-14), we expand it into a power series for a small increment of the asymptotic factor M  ,     A A B and c are known values; they follow from Equation (E-4), where the asymptotic factors , a M M and c M are to be found, , , Divide the second equation of Equation Set (F-4) over its first equation, Thus, we solve Equation Set (F-1) together with Equation (F-5), to obtain the three unknown asymptotic factors.To obtain the initial guess, we assume that the layer is thin, i.e. that the differences in the asymptotic factors Problem F2.RMS vs. time with the unknown top interface velocity a and top gradient a .Given data are the RMS velocities at the top and bottom interfaces, 2,a and 2,b , and at the inner point of the interval, 2,c , the full interval traveltime b (between the top and bottom interfaces) and the partial traveltime c (between the top interface and the intermediate level).Find the top velocity a and the gradient a .Solution.The resolving equation set follows from Equation (E-8), Note that for the small increments of the asymptotic factor Equation Set (F-16) should be solved numerically along with Equation (F-5).To obtain the initial guess, we linearize Equation Set (F-16) for the small increments of the asymptotic factors,

Appendix G. Numerical Examples of BVRT
In this appendix, we present three examples of the boundary value ray tracing, Figure 9.The parameters of the velocity profile are the same as above.For each case, the source point A is located at the origin of the frame.
The destination point is different for each case: B 1 (x = 2 km, z = 3 km), B 2 (x = 4 km, z = 2 km), and B 3 (x = 8 km, z = 2 km).As we will show, the first case corresponds to the pre-critical ray path, the second corresponds to the post-critical path that does not include the turning point, and the last case leads to the post-critical path with the turning point.

Figure 1 .Figure 2 .
Figure 1.Asymptotically bounded velocity models: Muscat hyperbolic model, Muscat exponential model, the Exponential slowness model, the Exponential asymptotically bounded model and the Conic model: Velocity vs depth.For all models, the profile parameters are: the top interface velocity V a = 3 km/s, the top gradient k a = 1 s −1 and the asymptotic velocity V ∞ = 6 km/s.
the use of Equations ( equation of Equation Set (32).The interval velocity (local average velocity) through the layer between the interfaces becomes

Figure 4 .
Figure 4. Pre-critical (red lines), critical (green lines) and post-critical (blue lines) ray trajectories for Hyperbolic (lines 1, 4, 7), Conic (lines 2, 5, 8) and EAB (lines 3, 6, 9) velocity models.Pre-critical rays propagate to an infinite depth and become asymptotically straight.Critical ray propagate to an infinite depth, and the ray angle approaches π/2 at large depth, but this ray has no asymptote.Postcritical rays pass the turning point and return to the earth surface, their trajectories are symmetric arcs.The Conic rays pass above the Hyperbolic model rays because their curvature is larger.For the same reason, the EAB rays pass above the Conic rays.For all models, parameters of the medium are the same as in Figure 1.The take-off angle of the pre-critical rays is 3π/24, that of the critical rays π/6, and that of the post-critical rays 5π/24.
Let a and b  be the ray angles at the start point and the destination point of the ray path, respectively.Equation Set (76) can be re-arranged as follows  For the pre-critical rays, the ray angles at the endpoints of the path.The following identities were used,

Figure 5 .
Figure 5. Traveltime vs. arc length of ray path for the three kinds of rays of the Hyperbolic velocity model: the precritical ray α a = 22.5˚, the critical ray α a = α C = 30˚ and the post-critical ray α a = 37.5˚.

where 2 1
m   m is the conjugate eccentricity of the post-critical ray path.Recall that sin sin .

7 ,
and Equation (99) yields the initial guess 3 m  .Solving Equation (98) with the Newton method, we obtain the eccentricity m

Figure 6 .
Figure 6.Full arcs of post-critical rays for the EAB, the Conic and the Hyperbolic velocity profiles.

Figure 7 .
Figure 7. Diagrams "Gradient-Velocity" for asymptotically bounded velocity models.Only for the EAB model the diagram is linear.For all models, except the Exponential slowness model, the gradient decreases with the increase of velocity.For the Exponential slowness model the gradient reaches maximum when the velocity becomes one half of the asymptotic value.Note the central symmetry between the two Muscat models, Hyperbolic and Exponential.

Figure 8 .
Figure 8. Lambert function y = L(x) is the solution of the transcendent equation y•exp(y) = x for a given value x and an unknown value y.The function has two branches.Branch zero is plotted in red, and branch minus one is plotted in blue.
can establish the corresponding parameters at the bottom interface.For this, we calculate the top asymptotic factor a M with the first equation of Equation Set (26).The bottom as-ymptotic factor b M can be obtained with the first equation of Equation Set (32).

k
Problem D1.Instantaneous velocities vs. depth.Given data are the top and bottom interface instantaneous velocities, a and b V , the asymptotic velocity V V  and the interval thickness z  .Find the top interface gradient a .Solution.Apply Equation (26) to calculate the top and bottom asymptotic factors, Instantaneous velocities vs. time.Given data are the top and bottom interface instantaneous velocities, a and b V , the asymptotic velocity V V  and the interval traveltime t  .Find the top interface gradient a .Solution.Apply Equation (26) to calculate the top and bottom asymptotic factors, k a M and b M .Problem D3.Gradients vs. depth.Given data are the top and bottom interface vertical gradients, and , .IJG the asymptotic velocity V  and the interval thickness .Find the top interface instantaneous velocity a V .Solution.Solve Equation Set (32) for the unknown asymptotic factors, the top interface gradient a .Solution.Apply Equation (26) to calculate the top asymptotic factor

6 )
With this notation, the second equation of Equation Set (32) be-comes 0.
Velocity and gradient vs. depth at different interfaces.Given data are the top interface gradient a , the bottom interface instantaneous velocity b , the asymptotic velocity V k V and the interval thickness z  .Find the top interface velocity a .Solution.Calculate the bottom asymptotic factor Velocity and gradient vs. time at different interfaces.Given data are the top interface velocity a , the bottom interface vertical gradient b , the asymptotic velocity V V k and the interval traveltime t  .Find the top interface gradient a .Solution.Apply Equation (26) to calculate the top asymptotic factor a k M .Then solve the nonlinear Equation (C-3) for the unknown bottom asymptotic factor b M .To get the initial guess, we assume 3) for a small increment of the asymptotic factor M 

(D- 12 )
Problem D9.Given data are the instantaneous velocity at the bottom interface, and at an intermediate b V Copyright © 2013 SciRes.IJG level inside the interval, c , and the asymptotic velocity .Two vertical distances are specified: the full interval thickness b (the distance between the top and bottom interfaces) and the partial interval thickness distance between the top interface and the intermediate level).Find the top interface velocity and gradient, a V and a .Solution.Use Equation (26) to calculate the bottom asymptotic factor b k M and the intermediate asymptotic factor c M .Apply Equation (29) for the full interval and for the partial interval, .Given data are the instantaneous velocity at the bottom interface, and at an intermediate level c inside the interval, and the asymptotic velocity .Two vertical traveltimes are specified: the full interval traveltime b (the traveltime between the top and bottom interfaces) and the partial traveltime between the top interface and the intermediate level).Find the top interface velocity and gradient, a and a .Solution.Use Equation (26) to calculate the bottom asymptotic factor b V k M and the intermediate asymptotic factor c M .Apply Equation (28) for the full interval and for the partial interval, parameter becomes, Equation (D-16) can be solved with the Lambert function, branch zero, RMS vs. depth with unknown gradient.Given data are the RMS velocities at the top and bottom interfaces, 2,a and 2,b , the interval thickness V V z  and the top interface velocity a .Find the top gradient a .Solution.It follows from the definition of the hyperbolic parameter, Equation (44), E2. RMS vs. depth with unknown velocity.Given data are the RMS velocities at the top and bottom interfaces, 2,a and 2,b , the interval thickness V V z  and the top interface gradient .Find the top interface a k velocity a .Solution.Use Equation (E-3) and the first equation of Equation Set (32),

.
assumption allows linearization of Equation Set (F-1), which, in turn, results in We solve Equation (F-7), choose a proper root, and this yields the initial guess for a M .Then we use Equation (F-6) to obtain the initial guess for b M and c M .Next we solve Equation Set (F-1), (F-5) numerically.
c and b are the normalized local RMS velocities, for the partial interval and the full interval, respectively.These are the known dimensionless values,

14 )
13) is similar to (F-5); however, (F-5) is exact and valid for any layer thickness, while (F-13) is an approximation for a thin layer only.Introduction of (F-13) into (F-11) results in a quadratic equation,The equation has two positive roots, but only one of them exceeds 1.The solution is .we solve Equation Set (F-8), (F-10) numerically and deliver the three asymptotic factors.Problem F3.Depth vs. time with the unknown top interface velocity a and top gradient a .Given data are the traveltime and layer thickness for the full and the partial intervals, respectively: Find the top velocity a and the gradient a .Solution.The resolving equation set follows from Equation (E-14), b and b are the normalized interval velocities for the full and the partial intervals,

Figure 9 .
Figure 9. Boundary value ray tracing with the Hyperbolic velocity model: AB 1 -pre-critical ray, AB 2 -post-critical ray without the turning point, AB 3 -post-critical turning ray, green line-critical ray, grey line-auxiliary post-critical ray with the turning point at the given depth.Introduction of Equation Set (F-18) into Equation (F-5) results in    1 . 1 c a a c b a a b D M M z D M M z        (F-19)Equation (F-19) yields the initial guess for the top asymptotic factor a M ,

3 AB
First we calculate the critical lateral propagation with Equation (101): Δx C = 1.587 km for the destination point depth z b = 2 km, and Δx C = 2.646 k km.Next we compare the horizontal offset m for z = 3 , while rays AB 2 and AB 3 are post-critical.We solve Equation (include the turning point or not.If the postcritical ray includes a turning point, then the ray angle at the destination is obtuse, and this should be taken into account in Equation (104).For this we establish an auxiliary post-critical ray, whose maximum penetration to the destination depth of points B 1Compare Equation (G-1) with Equation (22).It is interesting to note that the eccentricity of a post-critical ray with the turning point at the given depth is equal to the normalized velocity at this de e obtain the eccen-  and the top interface velocity a .Find the top gradient a .Solution.Use Equation (26) to find the top asymptotic factor