Abstract:
A well constrained horizontal variable height-velocity curve constructing method for seismic wave velocity field construction involves the steps of: a) calculating horizon velocity of each horizon by a sonic logging curve, and calculating the conversion horizon velocity of each horizon; b) drawing a circle to collect well points; c) calculating the characteristic parameter values of the height-velocity curve by the horizon velocity and the conversion horizon velocity of each well; d) calculating the characteristic parameter values by Kriging interpolation.

Description:
CROSS-REFERENCE TO RELATED APPLICATIONS 
       [0001]    This application is a continuation of International Application No. PCT/CN2010/000145, filed on Feb. 2, 2010, which claims the benefit of Chinese Patent Application No. 200910058920.8, filed Apr. 13, 2009. The disclosures of the above applications are incorporated herein by reference. 
     
    
     TECHNICAL FIELD 
       [0002]    The invention relates to a method for constructing seismic wave velocity field; more particularly, the invention relates to a method for constructing well-constrained transverse variable H-V curve in petroleum geophysical exploration. 
       BACKGROUND OF THE INVENTION 
       [0003]    Research for the method of constructing velocity field in a complex area is a frontier research subject that long exists in petroleum geophysical exploration, and has not been completely solved. Currently, many theoretical and practical application achievements have been obtained, but the research mostly focuses on the following two aspects: 
         [0004]    (1) Research on How to Obtain More Accurate Average Velocity Field with Seismic Data 
         [0005]    The research mainly includes two aspects: firstly, improving calculation accuracy of iterative velocity, with methods of iterative weighted complex trace velocity analysis, iterative weighted inclination velocity analysis, etc.; secondly, improving calculation accuracy of interval velocity, with methods such as generalized Dix method and various velocity smoothing methods, etc. Since there are many factors that affect seismic velocity, the above methods though improve accuracy of velocity analysis to some extent, but cannot guarantee accurate average velocity field. 
         [0006]    (2) Comprehensive Study on Average Velocity Field in Accordance with Seismic, Logging and VSP Data 
         [0007]    At present, said methods are main hotspots of research, and have achieved great research results and relatively wide applications. However, the research of such aspect mainly depends on an interpreter&#39;s experience, summarizing seismic velocity, logging and VSP data, which requires very heavy human workload. In addition, different interpreters would obtain seismic velocity field differing greatly from their own experience, and it is difficult to accurately constructing a subterranean velocity field. 
         [0008]    The propagation velocity of a seismic wave in strata is related to factors of lithology of rock, porosity, burial depth, pore-fluid property, pore fluid pressure and external environment of the strata, and the like, wherein the lithology is the most important factor. Strata of different eras present difference in terms of lithology and porosity, and thus in terms of velocity. Generally, from new eras to old ones, the interval velocity of a stratum gradually increases. 
         [0009]    In strata of the same era, particularly in marine bed, the lithology is common. Generally, transverse variation of velocity is caused by two factors. One is transverse variation of sedimentary facies, resulting in the lithology changes transversely, so that the interval velocity also changes transversely. Theoretically, the interval velocity of a single compact rock is close to a determined value, but actual strata are not composed of “single” component, but formed by mixing different components, such as argillaceous sandstone, calcareous sandstone, and the like, which makes the velocity of a seismic wave have mixed transitional nature, depending on contents of various components. Besides, porosity, style of pore fillings and pore fluid pressure and the like also affect velocity. The other factor affecting the transverse variation of velocity is depth (pressure). Generally, the deeper burial, the more compact, the smaller porosity, and the faster velocity. 
         [0010]    It is assumed that in the case the burial depth is the same, and the temperature, load pressure and the like conditions are the same or similar, relationship between velocity and depth, i.e. a compaction curve, can be calculated by a statistical method. When taking interval velocity of data of the respective measured points in a certain terrane as the X-axis, a depth thereof as the Y-axis, a normal velocity-depth curve (i.e. H-V curve) is obtained by fitting tendency that the velocity varies with the depth of the scattered points. Generally, the variation of velocity with depth can be described by a function of first degree, that is, the velocity is deemed to increase linearly with depth, represented by the following formula: 
         [0000]        V=ah+b    (1)
 
         [0011]    wherein v is interval velocity, h is strata burial depth, and a, b are characteristic parameters of the velocity-depth curve (i.e. H-V curve). 
         [0012]    Since there are lots of methods for acquiring subterranean velocity, there are lots of methods for acquiring a velocity-depth curve. For example, an interval velocity-depth curve can be obtained according to logging velocities of multiply wells; a velocity-depth curve can be obtained according to conversion interval velocities of multiple wells; or an interval velocity-depth curve can be obtained according to seismic velocity. Since there are many factors that affect the seismic velocity and the conversion velocity, among the plurality of curves, the interval velocity-depth curve obtained according to the logging interval velocities of multiple wells can best represent real subterranean conditions. 
         [0013]    As mentioned above, the variation of velocity with depth can be described with a function of first degree, then whether the one velocity-depth curve can describe transverse velocity variation in a survey or not? The answer is no. This is because there are two main factors affecting transverse variation of velocity. The tendency that velocity changes with depth (pressure) can be described by value a of the above-mentioned velocity-depth curve. In a relatively small survey, the value a in the velocity-depth curve is relatively constant, particularly in the marine bed, the value a is almost constant. Thus in a relatively small survey, a fixed value a can be used to describe the tendency that velocity changes with depth (pressure). However, in a relatively large survey, the value a should be variable. 
         [0014]    The other factor prominently affecting seismic interval velocity is transverse variation of sedimentary facies. Theoretically, the interval velocity of a single compact rock is close to a determined value, actually, strata are not composed of “single” component, but formed by mixing different components, such as argillaceous sandstone, calcareous sandstone, and the like, which makes the velocity of a seismic wave have mixed transitional nature, depending on contents of various components. Even if in a marine bed with relatively stable sedimentary facies, its ingredients are still different, but the velocity varies less relative to land strata. Besides, porosity, style of pore fillings and pore fluid pressure and the like also affect velocity. If the same function of first order is used to describe transverse variation of seismic velocity, the result is, even in a relatively small survey, the variation of value b is very large. 
         [0015]    A basic conclusion is obtained from the above discussion, that is, variation of velocity with depth can be described with a function of first order, and the H-V curve varies transversely, but the variation tendency of value a and that of value b in the curve are different, wherein the value a is relatively constant, varies little transversely; whereas the value b is affected by many factors, varies violently transversely. 
         [0016]    The above basic conclusion is commonly accepted, but the key point is how to apply the basic conclusion in practical complex velocity construction, i.e. how to construct a transverse variable H-V curve in the practical velocity construction. Especially, result of the complex velocity construction is generally to be used in practical time-depth conversion, the result of which should correspond to practical drilling result; therefore, construction of a transverse variable H-V curve has to be constrained by the drilling result. 
       SUMMARY OF THE INVENTION 
       [0017]    In order to solve the above technical problems, the invention provides a method for constructing a well-constrained transverse variable H-V curve in the construction of seismic wave velocity field. Applying the invention to the construction of the seismic wave velocity field can guarantee the consistence between the velocity and a drilling conversion velocity at the location of a drilling well; and at a location without any well, the velocity has a variation tendency that is consistent with that of a seismic velocity, thus guaranteeing that an accurate average velocity field is obtained. 
         [0018]    The invention is implemented by the following solutions. 
         [0019]    A method for constructing a well-constrained transverse variable H-V curve in the construction of seismic wave velocity field, wherein in the velocity-depth curve, i.e., the H-V curve, variation of velocity with depth is described with a function of first degree, and the velocity linearly increases with depth, represented by a formula v=ah+b, wherein v is interval velocity, h is strata burial depth, and a, b are characteristic parameters of the velocity-depth curve; characterized in that, the velocity-depth curve is the well-constrained transverse variable H-V curve which can be constructed by the following steps: 
         [0020]    A. calculating interval velocity of each strata (i.e. various geologic strata at a well location) in each well in a survey (i.e. the work area of seismic exploration) with an acoustic logging curve, and calculating a conversion interval velocity of each strata (i.e. various geologic strata at the well location) in each well with actual drilling result i.e. actual depth of the strata); 
         [0021]    B. with the coordinates of the well location in each strata in step A, drawing a circle with the point (coordinates of the well) as the center, with&#39;a certain radius (set by an operator, for example 3000 meters, etc.), and forming a set by all the well locations in the circle; if the number of wells in the set is smaller than a preset value, i.e. the minimal number of wells, enlarging the radius until the number of wells in the set is larger than or equals to the preset value, (generally the preset value is 4); 
         [0022]    C. fitting a function of first degree regarding the variation of velocity with depth by regression analysis, according to the logging interval velocity of each well in the set, taking the value a (determined from the formula v=ah+b) of the function as the value a of the well location; 
         [0023]    D. calculating value b from the formula v=ah+b according to the value a of the well location and the conversion interval velocity derived from step A, taking it as the value b of the well location; 
         [0024]    E. gridding (i.e. using certain gridding method, that is interpolation) the value a obtained from step C and value b from step D, respectively, according to Kriging interpolation algorithm, so as to obtain value a and value b of each well location in the whole survey, and forming the well-constrained transverse variable H-V curve in accordance with the formula v=ah+b. 
         [0025]    The method for applying the well-constrained transverse variable H-V curve constructed in the invention to the construction of seismic wave velocity field comprises the following steps: 
         [0026]    a. calculating proportional relationship between the value b of seismic interval velocity and the value b of drilling conversion velocity at the well location, i.e. proportional coefficient=value b of drilling conversion velocity/value b of seismic interval velocity, based on the constructed well-constrained transverse variable H-V curve and the established well-constrained transverse variable H-V curve,; 
         [0027]    b. obtaining a proportional coefficient of the whole survey by Kriging interpolation in transverse direction (which means each point extends in a plane); 
         [0028]    c. multiplying the value b of seismic interval velocity of the whole work area by said proportional coefficient, to obtain seismic and multiple-well-constrained value b of the whole work area; 
         [0029]    d. applying the value b obtained in step c in the above formula v=ah+b, velocity of each point can be obtained according to its value a, value b and depth data at each point, thus constructing the seismic wave velocity field. 
         [0030]    Said well-constrained transverse variable H-V curve is processed with velocity inversion algorithm, comprising the following steps: 
         [0031]    (1) using the comparatively interpreted horizon model and H-V curve, recursively deducing the initial depth of each CDP, by stripping layer by layer from shallow to deep in accordance with CDP point; 
         [0032]    (2) iterating thickness of each strata with H-V curve function v(h ij )=f(h ij ) (wherein i is serial number of CDP, j is serial number of layer), i.e. calculating the depth from the seismic reflection time and the H-V curve function, modifying velocity according to difference between the calculated depth and H-V template depth; which process continues until the depth error satisfies its accuracy, the interval velocity and bottom depth of each layer at the point is obtained when the iteration converges; 
         [0033]    (3) processing each CDP point in the same way, finally acquiring the interval velocity and bottom depth of each layer corresponding to each CDP. 
         [0034]    For layer j of CDP i, the H-V curve function of which is v(h ij )=f(h ij ). Provided that depth H i,j-1  of layer j−1 is known (the depth of layer 0 is fixed to zero), the inversion algorithm for calculating depth and velocity of layer j is as follows. 
         [0035]    (1) Giving the initial thickness h ij =ho ij  of layer j; 
         [0036]    (2) calculating the initial interval velocity v ij =f(ho ij ) of layer j from the function v(h ij )=f(h ij ) of H-V curve of layer j; 
         [0037]    (3) calculating theoretical reflection time difference Δt=ho ij /v ij  between layer j and layer j−1; 
         [0038]    (4) assuming actual reflection time difference between layer j and layer j-1 is ΔT (which can be calculated from travel time in that strata); if the difference between the theoretical reflection time and actual reflection time |Δt−ΔT|&lt;ε (ε is a small preset value, commonly taking ½ sampling rate), h ij  is the actual thickness of layer j, v ij  is the actual interval velocity of layer j; when the process ends, inversion result is output; 
         [0039]    (5) if Δt−ΔT&lt;0, indicating thickness h ij  is relatively small, making h ij =h ij +|Δt−ΔT|.v ij , and turning to step (2); otherwise, if Δt−ΔT&gt;0, indicating thickness h ij  is relatively large, making h ij =h ij −|Δt−ΔT|.v ij , and turning to step (2). 
         [0040]    In order to improve calculation efficiency of the inversion algorithm, the initial thickness h ij  of layer j can be the velocity of layer j of CDP i−1. Result of theoretical calculation and practical application shows that the algorithm has good convergence, and un-convergent condition is not met in our application. The method is not only suitable for time-depth conversion in gentle construction area, but also suitable for time-depth conversion in a highly steep and complex construction area. In addition, recursive inversion algorithm can determine and deal with fault (normal or reverse), totally controlled by the interpreted horizon model around the fault rather than affected by human factors, thus it can objectively recover shape of a footwall. 
         [0041]    Advantages of the invention lie in that, 
         [0042]    1. The velocity-depth curve constructed in the invention is a well-constrained transverse variable H-V curve, and the construction method is novel. Compared with the prior art, step C of the construction guarantees the variation tendency of velocity with depth to best correspond to actual conditions by using the value a obtained from a logging acoustic wave curve; use of conversion interval velocities of various wells in step D guarantees the velocity of a well location is totally identical to the actual conversion velocity, which thus guarantees the depth obtained from time-depth conversion conforms with actual drilling status; step E guarantees variation of value a and value b in the transverse direction is relatively gentle; in addition, the invention constructs seismic wave velocity field with steps a, b, c and d, compared with prior art, the invention guarantees consistence between the velocity and a drilling conversion velocity at the location of a drilling well in the velocity model; and at a location without a well, the velocity has a variation tendency that is consistent with that of a seismic velocity, thus guaranteeing that an accurate average velocity field is obtained. The invention further plays an active part in petroleum geophysical exploration. 
     
    
     DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 
     Embodiment 1 
       [0043]    The invention discloses a method for constructing well-constrained transverse variable H-V curve in the construction of seismic wave velocity field, wherein in the velocity-depth curve, i.e., the H-V curve, variation of velocity with depth is described with a function of first degree, and the velocity linearly increases with depth, represented by a formula v=ah+b, wherein v is interval velocity, h is strata burial depth, and a, b are characteristic parameters of the velocity-depth curve; characterized in that, the velocity-depth curve is the well-constrained transverse variable H-V curve which can be constructed by the following steps: 
         [0044]    A. calculating interval velocity of each strata (i.e. various geologic strata at a well location) in each well in a survey (i.e. the work area of seismic exploration) with an acoustic logging curve, and calculating a conversion interval velocity of each strata (i.e. various geologic strata at the well location) in each well with actual drilling result (i.e. actual depth of the strata); 
         [0045]    B. with coordinates of the well location in each strata in step A, drawing a circle with the point (coordinates of the well) as the center, with a certain radius (set by an operator, for example 3000 meters, etc.), and forming a set by all the well locations in the circle; if the number of wells in the set is smaller than a preset value, i.e. the minimal number of wells, enlarging the radius until the number of wells in the set is larger than or equals to the preset value, (generally the preset value is 4); 
         [0046]    C. fitting a function of first degree regarding the variation of velocity with depth by regression analysis, according to the logging interval velocity of each well in the set, taking the value a (determined from the formula v=ah+b) of the function as the value a of the well location; 
         [0047]    D. calculating value b from the formula v=ah+b according to the value a of the well location and the conversion interval velocity derived from step A, taking it as the value b of the well location; 
         [0048]    E. gridding (i.e. using certain gridding method, that is interpolation) the value a obtained from step C and value b from step D, respectively, according to Kriging interpolation algorithm, so as to obtain value a and value b of each well location in the whole survey, and forming the well-constrained transverse variable H-V curve in accordance with the formula v=ah+b. 
         [0049]    The method for applying the well-constrained transverse variable H-V curve constructed in the invention to the construction of seismic wave velocity field comprises the following steps: 
         [0050]    a. calculating proportional relationship between the value b of seismic interval velocity and the value b of drilling conversion velocity at the well location, i.e. proportional coefficient=value b of drilling conversion velocity/value b of seismic interval velocity, based on the constructed well-constrained transverse variable H-V curve and the established well-constrained transverse variable H-V curve; 
         [0051]    b. obtaining proportional coefficient of the whole survey a by Kriging interpolation in transverse direction (which means each point extends in a plane); 
         [0052]    c. multiplying the value b of seismic interval velocity of the whole work area by said proportional coefficient, to obtain seismic and multiple-well-constrained value b of the whole work area; 
         [0053]    d. applying the value b obtained in step c in the above formula v=ah+b, velocity of each point can be obtained according to its value a, value b and depth data at each point, thus constructing the seismic wave velocity field. 
       Embodiment 2 
       [0054]    Another preferred embodiment of the invention is described hereafter. 
         [0055]    1. In the prior art, establishment of an H-V curve is as follows. 
         [0056]    Propagation velocity of a seismic wave in strata is related to factors of lithology of rock, porosity, burial depth, pore-fluid property, pore fluid pressure and external environment of the strata, and the like, wherein the lithology is the most important factor. Strata of different eras present difference in terms of lithology and porosity, and thus in terms of velocity. Generally from young era to old ones, the interval velocity of a stratum gradually increases. 
         [0057]    In strata of the same era, particularly in marine bed, the lithology is common. Generally, transverse variation of velocity is caused by two factors. One is transverse variation of sedimentary facies, resulting in the lithology changes transversely, so that the interval velocity also changes transversely. Theoretically, the interval velocity of a single compact rock is close to a determined value, but actually strata are not composed of “single” component, but formed by mixing different components, for example argillaceous sandstone, calcareous sandstone, and the like, which makes the velocity of seismic wave have mixed transitional nature, depending on contents of various components. Besides, porosity, style of pore fillings and pore fluid pressure and the like also affect velocity. The other factor affecting the transverse variation of velocity is depth (pressure). Generally, the deeper burial, the more compact, the smaller porosity, and the faster velocity. 
         [0058]    It is assumed that in the case the burial depth is the same, and the temperature, load pressure and the like conditions are the same or similar, relationship between velocity and depth, i.e. compaction curve, can be calculated by a statistical method. When taking interval velocity of data of the respective measured points in a certain terrane as the X-axis, a depth thereof as the Y-axis, a normal velocity-depth curve (i.e. H-V curve) is obtained by fitting tendency that the velocity varies with the depth of the scattered points. Generally the variation of velocity with depth can be described by a function of first degree, that is, the velocity is deemed to increase linearly with depth, represented by the following formula: v=ah+b, wherein v is interval velocity, h is strata burial depth, and a, b are characteristic parameters of the velocity-depth curve (i.e. H-V curve). 
         [0059]    Since there are lots of methods for acquiring subterranean velocity, there are also lots of methods for acquire a velocity-depth curve. For example, an interval velocity-depth curve can be obtained according to logging interval velocities of multiply wells; an interval velocity-depth curve can be obtained according to conversion interval velocities of multiple wells; or an interval velocity-depth curve can be obtained according to seismic velocity. Since there are many factors that affect seismic velocity and conversion velocity, among the plurality of curves, the interval velocity-depth curve obtained according to the logging velocities of multiple wells can best represent real subterranean conditions. 
         [0060]    2. Well-Constrained Transverse Variable H-V Curve 
         [0061]    As mentioned above, the variation of velocity with depth can be described with a function of first degree, then whether the one velocity-depth curve can describe transverse velocity variation in a work area or not? The answer is no. This is because there are two main factors affecting transverse variation of velocity. The tendency that velocity changes with depth (pressure) can be described by value a of the above-mentioned velocity-depth curve. In a relatively small work area, the value a in the velocity-depth curve is relatively constant, particularly in the marine bed, the value a is almost constant. Thus in a relatively small work area, a fixed value a can be used to describe the tendency that velocity changes with depth (pressure). However, in a relatively large work area, the value a should be variable. 
         [0062]    The other factor prominently affecting seismic interval velocity is transverse variation of sedimentary facies. Theoretically, the interval velocity of a single compact rock is close to a determined value, but actual strata are not composed of “pure” components, but formed by mixing different components, for example argillaceous sandstone, calcareous sandstone, and the like, which makes the velocity of seismic wave have mixed transitional nature, depending on contents of various components. Even if in a marine bed with relatively stable sedimentary facies, its ingredients are still different, but the velocity varies less relative to land strata. Besides, porosity, style of pore fillings and pore fluid pressure and the like also affect velocity. If the same function of first order is used to describe transverse variation of seismic velocity, the result is even in a relatively small work area, the variation of value b is very large. 
         [0063]    A basic conclusion is obtained from the above discussion, that is, variation of velocity with depth can be described with a function of first order, and the H-V curve varies transversely, but the variation tendency of value a and that of value b in the curve are different, wherein the value a is relatively constant, varies little transversely; whereas the value b is affected by many factors, varies violently transversely. 
         [0064]    The above basic conclusion is commonly accepted, but the key point is how to apply the basic conclusion in practical complex velocity construction, i.e. how to construct a transverse variable H-V curve in the practical velocity construction. Especially, result of the complex velocity construction is generally to be used in practical time-depth conversion, the result of which should correspond to practical drilling result; therefore, construction of a transverse variable H-V curve has to be constrained by the drilling result. Thus we propose a method for constructing well-constrained transverse variable H-V curve, comprising the following steps: 
         [0065]    A. calculating interval velocity of each strata (i.e. various geologic strata at a well location) in each well in a survey (i.e. the work area of seismic exploration) with an acoustic logging curve, and calculating a conversion interval velocity of each strata (i.e. various geologic strata at the well location) in each well with actual drilling result (i.e. actual depth of the strata); 
         [0066]    B. with the coordinates of the well location in each strata in step A, drawing a circle with the point (coordinates of the well) as the center, with a certain radius (set by an operator, for example 3000 meters, etc.), and forming a set by all the well locations in the circle; if the number of wells in the set is smaller than a preset value, i.e. the minimal number of wells, enlarging the radius until the number of wells in the set is larger than or equals to the preset value, (generally the preset value is 4) 
         [0067]    C. fitting a function of first degree regarding the variation of velocity with depth by regression analysis, according to the logging interval velocity of each well in the set, taking the value a (determined from the formula v=ah+b) of the function as the value a of the well location; 
         [0068]    D. calculating value b from the formula v=ah+b according to the value a of the well location and the conversion interval velocity derived from step A, taking it as the value b of the well location; 
         [0069]    E. gridding (i.e. using certain gridding method, that is interpolation) the value a obtained from step C and value b from step D, respectively, according to Kriging interpolation algorithm, so as to obtain value a and value b of each well location in the whole work area. Thus H-V linear relation, i.e. v=ah+b, is obtained for each point in the survey. Because the value a and the value b at each location are different, the H-V curve is named traverse variable H-V curve, and the well-constrained transverse variable H-V curve is formed in accordance with the formula v=ah+b . 
         [0070]    3. Velocity Inversion Algorithm Constrained by a Transverse Variable H-V Curve 
         [0071]    In actual time-depth conversion, seismic refection time is known, but reflection depth is unknown, thus the above method cannot be directly cited. The seismic reflection time has to be transformed into reflection depth which is exactly the target of time-depth conversion. Therefore the invention provides a velocity inversion algorithm constrained by a transverse variable H-V curve, comprising the following steps: 
         [0072]    (1) using comparatively interpreted horizon model and H-V curve, recursively deducing the initial depth of each CDP, by stripping layer by layer from shallow to deep in accordance with CDP point; 
         [0073]    (2) iterating thickness of each strata with H-V curve function v(h ij )=f(h ij ) (wherein i is serial number of CDP, j is serial number of layer), i.e. calculating the depth from the seismic reflection time and the H-V curve function, modifying velocity according to difference between the calculated depth and H-V template depth; which process continues until the depth error satisfies its accuracy, the interval velocity and bottom depth of each layer at the point is obtained when the iteration converges; 
         [0074]    (3) processing each CDP point in the same way, finally acquiring the interval velocity and bottom depth of each layer corresponding to each CDP point. 
         [0075]    For layer j of CDP i, the H-V curve function of which is v(h ij )=f(h ij ). Provided that depth H ij−1  of layer j−1 is known (the depth of layer 0 is fixed to zero), the inversion algorithm for calculating depth and velocity of layer j is as follows. 
         [0076]    Giving the initial thickness h ij =ho ij  of layer j; 
         [0077]    calculating the initial interval velocity v ij =f(ho ij ) of layer j from the function v(h ij )=f(h ij ) of H-V curve of layer j; 
         [0078]    calculating theoretical reflection time difference Δt=ho ij /v ij  between layer j and layer j=1; 
         [0079]    assuming actual reflection time difference between layer j and layer j−1 is ΔT (which can be calculated from travel time in that strata); if the difference between the theoretical reflection time and actual reflection time |Δt−ΔT|&lt;ε (ε is a small preset value, commonly taking ½ sampling rate), h ij  is the actual thickness of layer j, v ij  is the actual interval velocity of layer j; when the process ends, inversion result is output; 
         [0080]    if Δt−ΔT&lt;0, indicating thickness h ij  is relatively small, making h ij =h ij +|Δt−ΔT|.v ij , and turning to step (2); otherwise, if Δt−ΔT&gt;0, indicating thickness h ij  is relatively large, making h ij =h ij −|Δt−ΔT|.v ij , and turning to step (2). 
         [0081]    In order to improve calculation efficiency of the inversion algorithm, the initial thickness h ij  of layer j can be the velocity of layer j of CDP i−1. Result of theoretical calculation and practical application shows that the algorithm has good convergence, and un-convergent condition is not met in our application. The method is not only suitable for time-depth conversion in a gentle construction area, but also suitable for time-depth conversion in a highly steep and complex construction area. In addition, recursive inversion algorithm can determine and deal with fault (normal or reverse), totally controlled by the interpreted horizon model around the fault rather than affected by human factors, thus it can objectively recover shape of a footwall. 
         [0082]    4. A method for Constructing Well-Constrained Transverse Variable H-V Curve in the Construction of Seismic Wave Velocity Field 
         [0083]    The key to realize the method for constructing well-constrained transverse variable H-V curve in the construction of seismic wave velocity field is to integrate seismic interval velocity and drilling conversion velocity, i.e. to integrate various velocities. Considering the problem like this, if the effect of depth (compactness) on velocity in seismic interval velocity and drilling conversion velocity is eliminated, main factor causing velocity variation in transverse direction is the transverse variation of sedimentary facies, which is variation of value b in the depth-velocity curve. The transverse variation of sedimentary facies is slow, thus variation of value b is also slow. If a proportional relationship between value b of seismic interval velocity and value b of drilling conversion velocity is described by a proportional coefficient, the proportional coefficient should vary slowly. Of course, the proportional coefficient is affected by various interference factors, and there may be mutation at certain position, but the total tendency is slowly variable. Regarding complex velocity model construction, the value b caused by pores, cracks, fillings and the like can be ignored, and only transverse variation of velocity caused by transverse variation of sedimentary facies is considered. It is derived from this the basic process for constructing seismic wave velocity field is as follows. 
         [0084]    According to years of experience, it is considered the variation of velocity with depth can be described by a function of first order, i.e. velocity increases linearly with depth, and the function of first order can eliminate the effect of depth (compactness) on velocity: 
         [0085]    a. based on the constructed well-constrained transverse variable H-V curve and the established well-constrained transverse variable H-V curve, calculating proportional relationship between the value b of seismic interval velocity and the value b of drilling conversion velocity at the well location, i.e. proportional coefficient=value b of drilling conversion velocity/value b of seismic interval velocity; 
         [0086]    b. obtaining a proportional coefficient of the whole survey through Kriging interpolation in transverse direction (which means each point extends in a plane); 
         [0087]    c. multiplying the value b of seismic interval velocity of the whole work area by said proportionality coefficient, to obtain seismic and multiple-well-constrained value b of the whole work area 
         [0088]    d. applying the value b obtained in step c in the above formula v=ah+b, to obtain velocity of each point according to its value a and value b and depth data at each point, thus constructing the seismic wave velocity field.