WO2010118624A1 - 地震波速度场构建中井约束横向可变h-v曲线构建方法 - Google Patents

地震波速度场构建中井约束横向可变h-v曲线构建方法 Download PDF

Info

Publication number
WO2010118624A1
WO2010118624A1 PCT/CN2010/000145 CN2010000145W WO2010118624A1 WO 2010118624 A1 WO2010118624 A1 WO 2010118624A1 CN 2010000145 W CN2010000145 W CN 2010000145W WO 2010118624 A1 WO2010118624 A1 WO 2010118624A1
Authority
WO
WIPO (PCT)
Prior art keywords
velocity
curve
value
depth
well
Prior art date
Application number
PCT/CN2010/000145
Other languages
English (en)
French (fr)
Inventor
李亚林
胡光岷
何光明
巫骏
黄东山
吴秋波
谢万学
Original Assignee
中国石油集团川庆钻探工程有限公司
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by 中国石油集团川庆钻探工程有限公司 filed Critical 中国石油集团川庆钻探工程有限公司
Priority to EP10764037.7A priority Critical patent/EP2420862A4/en
Priority to EA201101493A priority patent/EA020635B1/ru
Publication of WO2010118624A1 publication Critical patent/WO2010118624A1/zh
Priority to US13/273,096 priority patent/US8868348B2/en

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/303—Analysis for determining velocity profiles or travel times
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/66—Subsurface modeling
    • G01V2210/665—Subsurface modeling using geostatistical modeling
    • G01V2210/6652—Kriging

Definitions

  • the invention relates to a method for constructing a subsurface seismic wave velocity field, and more particularly to a method for constructing a well constrained lateral variable H-V curve for petroleum seismic exploration.
  • the main research includes two parts: First, the calculation accuracy of the superposition speed is improved.
  • the main methods are weighted iterative complex multi-path velocity analysis, weighted iterative tilt velocity analysis, etc.
  • the second is to improve the accuracy of layer velocity calculation, such as the generalized Dix method and Various speed smoothing methods, etc. Due to the many influencing factors of seismic velocity, although the above method improves the accuracy of velocity analysis to a certain extent, it does not guarantee a truly accurate average velocity field.
  • this kind of method is the main research hotspot, and many good research results have been obtained, and it has been widely used.
  • the research in this area relies mainly on the experience of the interpreter to synthesize the seismic velocity and the logging and VSP data.
  • the manual workload is very large, and different interpreters can obtain different seismic velocity fields based on their experience.
  • Accurate subsurface velocity field The velocity of seismic waves propagating in the rock formation is related to the lithology, porosity, burial depth, pore fluid properties, pore fluid pressure and external environmental conditions of the rock formation. The most important of these is the influence of lithology.
  • the strata of different eras differ in lithology and porosity, and therefore there are also differences in speed. Generally, from the new to the old, the layer velocity of the stratum gradually increases.
  • lithology In contemporary strata, especially in marine strata, lithology is interlinked. In general, there are two main factors that cause lateral changes in velocity. First, due to lateral changes in the sedimentary facies, the lithology changes laterally, which causes lateral velocity changes. Theoretically, the layer velocity of a single uniform dense rock is close to a certain value, but the actual rock formation is often not composed of "pure" components, but is composed of different components, such as argillaceous sandstone, gray sandstone, etc., which makes seismic waves The speed also has the nature of a mixed transition and depends on the content of the individual components. In addition, porosity, pore filler type, and pore fluid pressure can also affect speed. Another major factor affecting the lateral change in velocity is depth (pressure). Under normal conditions, the deeper the burial, the greater the compaction, the smaller the porosity and the more turbulent.
  • V is the formation velocity
  • h is the formation depth
  • a and b are the characteristic parameters of the velocity-depth curve (ie, the H-V curve).
  • a layer velocity-depth curve can be obtained according to the logging velocity of multiple wells; a layer velocity-depth curve can be obtained according to the converted velocity layer velocity of multiple wells; we can also obtain a layer velocity according to the seismic velocity.
  • a depth curve Since the seismic velocity and the conversion speed have many influence factors, in these multiple curves, the formation velocity-depth curve based on the logging velocity of multiple wells can best reflect the actual situation in the underground.
  • Another factor that significantly affects the velocity of the seismic layer is the lateral variation of the sedimentary facies.
  • the layer velocity of a single dense rock is close to a certain value, but the actual rock formation is often not composed of "pure" components, but by different components.
  • Mixed such as argillaceous sandstone, gray sandstone, etc., which makes the velocity of the seismic wave also have the nature of mixed transition, and depends on the content of each component.
  • porosity, pore filler type, and pore fluid pressure can also affect speed. If a one-time function is still used to describe the lateral variation of the seismic velocity, we can see that the b-value change is large even in a relatively small work area.
  • the present invention proposes a method for constructing a well-constrained laterally variable HV curve in the construction of a seismic wave velocity field.
  • the application of the present invention to the construction of a seismic wave velocity field ensures that the drilling position and speed are consistent with the drilling conversion speed; In the case of no well, the trend of the velocity is consistent with the trend of the earthquake velocity, thus ensuring a truly accurate average velocity field.
  • the steps for building the curve are as follows:
  • step A For each layer in step A, the coordinates of the well point of each well, with the point (coordinate of the well) as the center, with a certain radius (required by the operator, such as 3000 meters, etc.) a circle that collects all the well points in the circle to form a set. If the number of wells in the set is less than a predetermined value, that is, the minimum number of wells, the radius is expanded until the number of wells in the set is greater than or equal to a predetermined one. a value up; (usually 4)
  • the a value obtained in the C step and the b value obtained in the D step are respectively meshed (the meshing means using a certain meshing method, that is, interpolation), and the whole is obtained.
  • the Kriging interpolation is used to obtain the scale factor of the entire work area
  • the well-constrained laterally variable H-V curve is processed by a speed inversion algorithm, and the processing steps are as follows:
  • the speed inversion algorithm is as follows:
  • the actual reflection time difference between the jth layer and the j-1 layer be ⁇ (calculated from the horizon time obtained by the layer position at that time, if the difference between the theoretical reflection time and the actual reflection time
  • the initial thickness of the jth layer can be the speed of the jth layer of the i-1th CDP.
  • the theoretical calculation and practical application results show that the convergence of the algorithm is very good. The application has not found no convergence.
  • This method is not only suitable for the time-depth conversion of the gentle structure, but also for the time-depth conversion of the high-steep complex structure area.
  • the recursive inversion algorithm can judge and deal with the fault (positive or inverse), which is completely explained by the fault near the fault.
  • the layer model control is not affected by human factors, so the shape of the fault under the disk can be objectively restored.
  • the speed-depth curve adopted by the present invention is a well-constrained laterally variable HV curve constructed, the construction method is novel, and compared with the prior art, the logging acoustic wave curve is used in the construction step C to obtain a value to ensure the speed
  • the depth change trend is in good agreement with the actual situation;
  • the conversion layer velocity of each well in step D ensures that the velocity at the well point is exactly the same as the actual conversion speed, which ensures that the depth obtained by the time-depth conversion is in line with the actual drilling situation;
  • E guarantees that the lateral "value and & value changes are relatively flat; at the same time, because the invention uses the four steps of &, b, c, d to construct the seismic wave velocity field, compared with the prior art, the speed model can be ensured at the drilling position.
  • the speed is consistent with the drilling conversion speed; in the position without the well, the change trend of the speed is consistent with the trend of the seismic speed, thus ensuring a truly accurate average velocity field.
  • the invention can
  • the invention discloses a method for constructing a well-constrained laterally variable HV curve in the construction of a seismic wave velocity field.
  • the velocity-depth ie, the HV curve
  • the steps to construct a constrained laterally variable HV curve are as follows:
  • step A For each layer in step A, the coordinates of the well point of each well, with the point (coordinate of the well) as the center, with a certain radius (required by the operator, such as 3000 meters, etc.) a circle that collects all the well points in the circle to form a set. If the number of wells in the set is less than a predetermined value, that is, the minimum number of wells, the radius is expanded until the number of wells in the set is greater than or equal to a predetermined one. One value up; (usually 4)
  • the Kriging interpolation is used to obtain the scale factor of the entire work area
  • the velocity of seismic waves propagating in the rock formation is related to the lithology, porosity, burial depth, pore fluid properties, pore fluid pressure and external environmental conditions of the rock. The most important of these is the influence of lithology.
  • the strata of different eras differ in lithology and porosity, and therefore there are also differences in speed. Generally, from the new to the old, the layer velocity of the stratum gradually increases.
  • lithology In contemporary strata, especially in marine strata, lithology is interlinked. In general, there are two main factors that cause lateral changes in velocity. First, due to lateral changes in the sedimentary facies, the lithology changes laterally, which causes lateral velocity changes. Theoretically, the layer velocity of a single uniform dense rock is close to a certain value, but the actual rock formation is often not composed of "pure" components, but is composed of different components, such as argillaceous sandstone, gray sandstone, etc., which makes seismic waves The speed also has the nature of a mixed transition and depends on the content of the individual components. In addition, porosity, pore filler type, and pore fluid pressure can also affect speed. Another major factor affecting the lateral variation of velocity is depth (pressure). Under normal conditions, the deeper the burial, the greater the compaction, the smaller the porosity and the higher the velocity.
  • the relationship between velocity and depth can be obtained by statistical methods, that is, the compaction curve.
  • the point data is the horizontal velocity as the horizontal coordinate
  • the depth is the ordinate as the scatter point
  • a curve is used to fit the trend of the velocity with the depth, which constitutes the normal speed-depth curve.
  • Line ie HV curve.
  • the characteristic parameter of a depth curve ie HV curve).
  • a layer velocity-depth curve can be obtained according to the logging velocity of multiple wells; a layer velocity-depth curve can be obtained according to the converted velocity layer velocity of multiple wells; Speed-depth curve. Because of the many factors affecting seismic velocity and conversion speed, we believe that: In these multiple curves, the formation velocity-depth curve based on the logging velocity of multiple wells best reflects the actual situation in the underground.
  • Another factor that significantly affects the velocity of the seismic layer is the lateral variation of the sedimentary facies.
  • the layer velocity of a single dense rock is close to a certain value, but the actual rock formation is often not composed of "pure" components, but by different components.
  • Mixed such as argillaceous sandstone, gray sandstone, etc., which makes the velocity of the seismic wave also have the nature of mixed transition, and depends on the content of each component.
  • porosity, pore filler type, and pore fluid pressure can also affect speed. If a function is still used to describe the lateral variation of the seismic velocity, we can see that even within a relatively small work area, the 6-value change is large.
  • step A For each layer in step A, the coordinates of the well point of each well, with the point (coordinate of the well) as the center, with a certain radius (required by the operator, such as 3000 meters, etc.) a circle that collects all the well points in the circle to form a set. If the number of wells in the set is less than a predetermined value, that is, the minimum number of wells, the radius is expanded until the number of wells in the set is greater than or equal to a predetermined one. a value up; (usually 4)
  • the a value obtained in the C step and the b value obtained in the D step are respectively meshed (the meshing means using a certain meshing method, that is, interpolation), and the whole is obtained.
  • HV curve function iterate the thickness of each layer (where i is the CDP serial number and j is the horizon number), that is, the depth is calculated from the seismic reflection time and the HV curve function, according to the difference between the calculated depth and the ⁇ - ⁇ plate depth. The speed is modified. This process proceeds until the depth error satisfies the accuracy. When the iteration converges, the layer velocity and the bottom boundary depth of each layer at that point are obtained.
  • the Hv curve function is ⁇ .) - / ⁇ ) , assuming that the depth of the ji layer is H, ., ⁇ is known (the depth of the 0th layer is fixed to zero), and the jth layer is calculated.
  • ⁇ £ is a predetermined smaller value, usually ⁇ The 1/2 of the sample rate, at this time is the actual thickness of the jth layer, ⁇ is the actual layer speed of the jth layer.
  • ⁇ 3 ⁇ 4., turn to step (2); if ⁇ - ⁇ >0, it means that the thickness is too large, let h u h Ir ⁇ T ⁇ ' Vij , turn to step (2); to improve the calculation speed of the inversion algorithm, the initial thickness of the jth layer /.
  • the key to constructing the well-constrained laterally variable HV curve in the construction of seismic wave velocity field is to combine the seismic layer velocity and the drilling conversion velocity, that is, to achieve the fusion of different types of velocity.
  • the proportionality coefficient should also be gradually variable.
  • this proportional coefficient is affected by various interference factors, and there may be mutations in individual parts, but the overall trend should be slow.
  • the values caused by various factors such as pores, cracks and fillings can be neglected, and only the lateral changes in velocity caused by the lateral variation of the sedimentary facies are considered.
  • the proportional relationship between the b value of the seismic layer velocity and the b value of the drilling conversion velocity is calculated at the drilling location, ie : scale factor - b value of drilling conversion speed / b value of seismic layer velocity;
  • the Kriging interpolation is used to obtain the scale factor of the entire work area

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Description

地震波速度场构建中井约束横向可变 H-V曲线构建方法
技术领域
本发明涉及一种地下地震波速度场构建方法， 确切地说涉及一种用于 石油地震勘探中的井约束横向可变 H-V曲线构建方法。
背景技术
复杂地区速度场构建方法研究是石油地球物理勘探中长期存在并在目 前尚未完全解决的前沿性研究课题， 目前研究取得了很多理论和实际应用 成果。 但目前的研究大都集中于以下两个方面：
( 1 ) 研究如何使用地震资料获得更精确的平均速度场
其主要的研究包括两部分内容： 一是提高叠加速度的计算精度， 主要 的方法有加权迭代复数道速度分析、 加权迭代倾斜速度分析等； 二是提高 层速度计算的精度， 如广义 Dix方法和各种速度平滑方法等。 由于地震速 度的影响因素很多， 上述方法虽然在一定程度上提髙了速度分析的精度， 但并不能保证获得真正精确的平均速度场。
( 2 ) 根据地震、 测井和 VSP资料综合研究平均速度场
目前该类方法是主要的研究热点， 已取得了许多较好的研究成果， 获 得了较为广泛的应用。 但这方面的研究主要依靠解释人员的经验进行地震 速度与测井和 VSP资料综合， 人工工作量很大， 而且不同的解释人员根据 他们的经验会得出差异较大的地震速度场， 很难做到准确的地下速度场建 地震波在岩层中传播的速度与岩石的岩性、 孔隙度、 埋藏深度、 孔隙 流体性质、 孔隙流体压力及岩层所处的外部环境条件等因素有关。 其中最 主要的是岩性的影响。不同时代的地层在岩性、孔隙度等方面存在差异， 因 而在速度上也有差异， 一般由新到老， 地层的层速度逐渐增大。
在同时代的地层中， 尤其在海相地层中， 岩性是相通的。 一般说来， 有两种主要因素引起速度的横向变化， 一是由于沉积相带的横向变化， 造 成岩性在横向上发生变化， 这使层速度也发生横向变化。 理论上单一致密 岩石的层速度接近于一个确定值， 但实际岩层往往并非由 " 纯"的成分组 成， 而是由不同成分混合而成， 如泥质砂岩、 灰质砂岩等等， 这使地震波 的速度也具有混合过渡的性质， 并取决于各个组分的含量。 此外孔隙度、 孔隙充填物类型及孔隙流体压力等也会对速度产生影响。 另一种影响速度 横向变化的主要因素是深度 (压力）， 在通常情况下， 埋藏越深， 压实程度 越大， 孔隙度也越小， 速度越髙。
我们可以假定在埋深相同的情况下， 温度、 负载压力等条件是相同的 或近似相同的， 可由统计的办法来求得速度与深度的关系， 即压实曲线， 若将某一岩层各测点的数据以层速度为横坐标， 深度为纵坐标作散点， 再 用一条曲线来拟合速度随深度变化的趋势， 这就构成了正常的速度一深度 曲线 （即 H-V 曲线）。 通常速度随深度的变化可以釆用一次函数来描述， 即认为速度随深度线性增加， 用如下公式来表示：
v=ah+b ( 1 )
其中 V为地层速度， h为地层埋深， a、 b为速度一深度曲线 （即 H-V 曲线） 的特征参数。
由于地下速度的获取方法很多， 因而获得速度一深度曲线的方法也很 多。 如根据多口井的测井层速度可以获得一根层速度一深度曲线； 根据多 口井的换算速度层速度可以获得一根层速度一深度曲线； 根据地震速度我 们也可以获得一根层速度一深度曲线。 由于地震速度和换算速度的影响因 素很多， 因此在这多根曲线中， 根据多口井的测井层速度获得地层速度一 深度曲线最能反映地下实际情况。
上面已经谈到速度随深度的变化可以釆用一次函数来描述， 那么是否 可以用一个速度一深度曲线来描述一个工区内的横向速度变化呢？ 回答显 然是否定的。 这是因为有两种主要因素引起速度的横向变化， 其中速度随 深度 (压力)变化的趋势可以上述讨论的速度一深度曲线中的 a值来描述。 认为在一个相对较小的工区中， 速度一深度曲线中的 a值是相对稳定的， 尤其在海相地层中， a 值的稳定性更强。 因而在一个较小的工区内可以釆 用固定的 a值来描述速度随深度 (压力)变化的趋势， 但在一个较大的工区 内 a值应该是变化的。
另一种显著影响地震层速度的因素是沉积相带的横向变化， 理论上单 一致密岩石的层速度接近于一个确定值， 但实际岩层往往并非由 "纯" 的 成分组成， 而是由不同成分混合而成， 如泥质砂岩、灰质砂岩等等， 这使地 震波的速度也具有混合过渡的性质， 并取决于各个组分的含量。 即使在沉 积相对稳定的海相地层中， 仍然存在成份的差异， 只是相对于陆相地层来 说， 变化很小而己。 此外孔隙度、 孔隙充填物类型及孔隙流体压力等也会 对速度产生影响。 如果仍然采用一次函数来描述地震速度的横向变化， 我 们可以发现即使在一个相对较小的工区内 b值的变化也是很大的。
有上述讨论我们可以得出一个基本的结论： 速度随深度的变化可以釆 用一次函数来描述， 并且 H-V曲线在横向上是变化的， 但曲线中的 a值和 b值变化趋势不同， a值是相对稳定的， 在横向上变化较小； b的影响因素 较多， 在横向上变化也较为剧烈。
应该说上述基本结论是得到普遍认可的， 但问题的关键在于如何将上 述基本结论应用于实际的复杂构造速度构建中， 也就是在实际的速度构建 中如何构建横向可变的 H-V曲线。尤其是复杂构造速度构建结果往往要用 于实际的时深转换， 时深转换的结果要与实际的钻井结果相吻合， 因此必 须采用钻井结果约束横向可变的 H-V曲线的构建。 发明内容
为解决上述技术问题， 本发明提出了一种地震波速度场构建中井约束 横向可变 H-V曲线构建方法，将本发明应用于地震波速度场构建中能确保 在钻井位置， 速度与钻井换算速度一致； 在无井的位置， 速度的变化趋势 与地震速度的变化趋势一致， 从而能保证获得真正精确的平均速度场。
本发明是通过如下技术方案实现的：
一种地震波速度场构建中井约束横向可变 H-V 曲线构建方法， 速度- 深度即 H-V曲线中， 速度随深度的变化采用一次函数来描述， 且速度随深 度线性增加， 用公式 v=ah+b表示， 其中 V为地层速度， h为地层埋深， a， b为速度-深度曲线的特征参数； 其特征在于：所述速度 -深度曲线为井约束 横向可变 H-V曲线， 井约束横向可变 H-V曲线的构建步骤如下：
A、 用声波测井曲线计算出工区内 （工区指地震勘探的工作区域） 各 井各层 （石油钻井的各个地质层位） 的层速度， 用实际的钻井结果 （地层 的实际深度）计算出各井各层（石油钻井的各个地质层位）的换算层速度；
B、对于 A步骤中的每一个层位、每一口井的井点的坐标， 以该点（井 的坐标）为圆心， 以一定的 （需要由操作人员给定， 如 3000米等）半径画 一个圆， 将该圆内的所有钻井点集中起来形成一个集合， 若集合中井的数 量小于预先给定的一个值即最小井数量， 则扩大半径， 直到集合中井的数 量大于或等于预先给定的一个值为止； （通常为 4 )
C、 根据集合中各井的测井层速度， 采用回归分析法拟合速度随深度 变化的一次函数， 将该函数的 a值 （根据公式 v=ah+b确定）， 作为该井点 的 a值；
D、 根据该井点的 a值和由 A步骤中得出该井点的换算层速度， 根据 公式 v=ah+b计算出 b值， 作为该井点的 b值；
E、 用克里金插值算法， 分别对 C步骤中得到的 a值和 D步骤中得到 的 b值进行网格化 （网格化是指采用一定的网格化方法，即插值）， 得到整 个工区内每一个井点的 a值和 b值， 根据公式 v=ah+b形成井约束横向可 变 H-V曲线。
将采用本发明构建的井约束横向可变 H-V 曲线应用于地震波速度场 构建的方法如下：
a、 基于已构建的井约束横向可变 H-V曲线， 基于已建立的井约束横 向可变 H-V曲线， 在钻井位置计算地震层速度的 b值与钻井换算速度的 b 值之间的比例关系， SP : 比例系数-钻井换算速度的 b值 /地震层速度的 b 值；
b、 在横向上 （横向是指平面上延伸的每一个点） 通过克里金插值获 得整个工区的比例系数；
c、用整个工区的地震层速度 b值与所述比例系数相乘，得到整个工区 地震与多井约束的 b值； W d,将 c步骤中得到的 b值应用于上述公式 v=ah+b中，在每一个点上， 根据其 a值、 b值以及深度数据， 可获得每个点的速度， 由此就构建了地 震波速度场。
' 所述的井约束横向可变 H-V曲线是经过速度反演算法处理的，其处理 步骤如下：
(1)应用对比解释的层位模型和 H-V曲线， 按 CDP点用逐层剥离的 方法， 由浅向深逐层递推该点每一层初始深度；
(2) 应用 H-V 曲线函数 (/¾) = /( )迭代每层厚度， 其中 i为 CDP 序号、 j为层位序号， 即由地震反射时间和 H-V曲线函数计算深度， 根据 计算深度与 H-V量板深度的差修改速度， 这一过程进行到深度误差满足精 度为止， 迭代收敛时就取得该点每层的层速度和底界深度；
(3)对每个 CDP点做相同处理， 最终可取得每个 CDP点对应的每层 的层速度和底界深度。
对于第 i个 CDP的第 j层， 其 H-V曲线函数为 ( ) = /(/¾.)， 假设 j-1 层的深度 已知（第 0层的深度固定为零）， 计算第 j层深度和速度的反 演算法如下：
(1) 给定第 j层的初始厚度 -/^.;
(2) 根据第 j层的 H-V曲线函数 / -/^)计算 j层的初始层速度
(3) 计算第 j层与 j-1层的理论反射时间差
Figure imgf000005_0001
(4)设第 j层与 j-1层的实际反射时间差为 ΔΓ (可由层位届时获得的 层位时间计算得到，若理论反射时间与实际反射时间的差 |At一 Δ:Τ|<^(^为 预先给定的一个较小的值， 通常取采样率的 1/2)， 这时的 就是第 j层的 实际厚度， 就是第 j层的实际层速度。 算法结束， 输出反演结果。
(5)若 At- ΔΓ<0， 表明厚度 偏小， 令/ ¾= +|At- ΔΓ|· ，转第（2) 步； 若 Δ - Δ >0， 表明厚度 偏大， 令 |Δ - 转第 （2)步； 为提高反演算法的计算速度，可令第 j层的初始厚度 为第 i-1个 CDP 的第 j层的速度， 理论计算和实际应用的结果表明： 该算法的收敛性很好， 在我们的应用中尚未发现不收敛的情况。 这种方法不仅适合平缓构造的时 深转换， 同时也适合高陡复杂构造区的时深转换， 同时递推反演算法能对 断层 (正或逆)作出判断并处理， 完全由断层附近解释的层位模型控制， 不 受人为因素影响， 因此可客观恢复断层下盘形态。
本发明的优点表现在：
1、 由于本发明采用的速度-深度曲线为构建的井约束横向可变 H-V曲 线， 构建方法新颖， 与现有技术相比， 本构建步骤 C中采用测井声波曲线 得到"值保证了速度随深度变化趋势与实际情况最吻合； 步骤 D中采用各 井的换算层速度保证了在井点处的速度与实际换算速度完全一致， 这就保 证了时深转换得到的深度符合实际钻井情况；步骤 E保证了横向的 "值和 & 值变化比较平缓； 同时， 由于本发明釆用 &、 b、 c、 d四个步骤构建地震波 速度场， 与现有技术相比， 能确保速度模型在钻井位置， 速度与钻井换算 速度一致； 在无井的位置， 速度的变化趋势与地震速度的变化趋势一致， 从而能保证获得真正精确的平均速度场。 采用本发明能.在石油地震勘探过 程中发挥更一步的积极作用。
具体实施方式
实施例 1
本发明公开了一种地震波速度场构建中井约束横向可变 H-V 曲线构 建方法， 速度-深度即 H-V 曲线中， 速度随深度的变化釆用一次函数来描 述， 且速度随深度线性增加， 用公式 v=¾h+b表示， 其中 V为地层速度， h 为地层埋深， a， b为速度-深度曲线的特征参数； 其特征在于： 所述速度- 深度曲线为井约束横向可变 H-V 曲线， 井约束横向可变 H-V 曲线的构建 步骤如下：
A、 用声波测井曲线计算出工区内 （工区指地震勘探的工作区域） 各 井各层 （石油钻井的各个地质层位） 的层速度， 用实际的钻井结果 （地层 的实际深度）计算出各井各层（石油钻井的各个地质层位）的换算层速度；
B、对于 A步骤中的每一个层位、每一口井的井点的坐标， 以该点（井 的坐标）为圆心， 以一定的 （需要由操作人员给定， 如 3000米等）半径画 一个圆， 将该圆内的所有钻井点集中起来形成一个集合， 若集合中井的数 量小于预先给定的一个值即最小井数量， 则扩大半径， 直到集合中井的数 量大于或等于预先给定的一个值为止； （通常为 4)
C、 根据集合中各井的测井层速度， 采用回归分析法拟合速度随深度 变化的一次函数， 将该函数的 a值 （根据公式 v=ah+b确定）， 作为该井点 的 a值；
D、 根据该井点的 a值和由 A步骤中得出该井点的换算层速度， 根据 公式 v=ah+b计算出 b值， 作为该井点的 b值；
E、 用克里金插值算法， 分别对 C歩骤中得到的 a值和 D歩骤中得到 的 b值进行网格化 （网格化是指釆用一定的网格化方法，即插值）， 得到整 个工区内每一个井点的 a值和 b值， 根据公式 v=ah+b形成井约束横向可 变 H-V曲线。
将采用本发明构建的井约束横向可变 H-V 曲线应用于地震波速度场 构建的方法如下：
a、 基于已构建的井约束横向可变 H-V曲线， 基于已建立的井约束横 向可变 H-V曲线， 在钻井位置计算地震层速度的 b值与钻井换算速度的 b 值之间的比例关系， 即： 比例系数=钻井换算速度的 b值 /地震层速度的 b 值；
b、 在横向上 （横向是指平面上延伸的每一个点） 通过克里金插值获 得整个工区的比例系数；
c、用整个工区的地震层速度 b值与所述比例系数相乘， 得到整个工区 地震与多井约束的 b值；
d、将 c步骤中得到的 b值应用于上述公式 v=ah+b中，在每一个点上， 根据其 a值、 b值以及深度数据， 可获得每个点的速度， 由此就构建了地 震波速度场。
实施例 2
作为本发明的另一较佳实施方式如下- 1、 现有技术中， H-V曲线的建立：
地震波在岩层中传播的速度与岩石的岩性、 孔隙度、 埋藏深度、 孔隙 流体性质、 孔隙流体压力及岩层所处的外部环境条件等因素有关。 其中最 主要的是岩性的影响。不同时代的地层在岩性、孔隙度等方面存在差异， 因 而在速度上也有差异， 一般由新到老， 地层的层速度逐渐增大。
在同时代的地层中， 尤其在海相地层中， 岩性是相通的。 一般说來， 有两种主要因素引起速度的横向变化，一是由于沉积相带的横向变化,造成 岩性在横向上发生变化， 这使层速度也发生横向变化。 理论上单一致密岩 石的层速度接近于一个确定值， 但实际岩层往往并非由" 纯"的成分组成, 而是由不同成分混合而成， 如泥质砂岩、 灰质砂岩等等， 这使地震波的速 度也具有混合过渡的性质， 并取决于各个组分的含量。 此外孔隙度、 孔隙 充填物类型及孔隙流体压力等也会对速度产生影响。 另一种影响速度横向 变化的主要因素是深度 (压力）， 在通常情况下， 埋藏越深， 压实程度越大， 孔隙度也越小， 速度越高。
我们可以假定在埋深相同的情况下， 温度、 负载压力等条件是相同的 或近似相同的， 可由统计的办法来求得速度与深度的关系， 即压实曲线,若 将某一岩层各测点的数据以层速度为横坐标， 深度为纵坐标作散点， 再用 一条曲线来拟合速度随深度变化的趋势， 这就构成了正常的速度一深度曲 线 （即 H-V 曲线）。 通常速度随深度的变化可以采用一次函数来描述， 即 认为速度随深度线性增加， 用如下公式来表示： v=ah+b， 其中 V为地层速 度， h为地层埋深， a,b为速度一深度曲线 （即 H-V曲线） 的特征参数。
由于地下速度的获取方法很多， 因而获得速度一深度曲线的方法也很 多。 如根据多口井的测井层速度可以获得一根层速度一深度曲线； 根据多 口井的换算速度层速度可以获得一根层速度一深度曲线； 更具地震速度我 们也可以获得一根层速度一深度曲线。 由于地震速度和换算速度的影响因 素很多， 因此我们认为： 在这多根曲线中， 根据多口井的测井层速度获得 地层速度一深度曲线最能反映地下实际情况。
2、 井约束的横向可变 H-V曲线
上面已经谈到速度随深度的变化可以采用一次函数来描述， 那么是否 可以用一个速度一深度曲线来描述一个工区内的横向速度变化呢？ 回答显 然是否定的。 这是因为有两种主要因素引起速度的横向变化， 其中速度随 深度 (压力)变化的趋势可以用上述讨论的速度一深度曲线中的 a值来描述。 认为在一个相对较小的工区中， 速度一深度曲线中的 a值是相对稳定的， 尤其在海相地层中， a 值的稳定性更强。 因而在一个较小的工区内可以采 用固定的 a值来描述速度随深度 (压力)变化的趋势， 但在一个较大的工区 内 a值应该是变化的。
另一种显著影响地震层速度的因素是沉积相带的横向变化,理论上单 一致密岩石的层速度接近于一个确定值， 但实际岩层往往并非由" 纯"的 成分组成， 而是由不同成分混合而成， 如泥质砂岩、灰质砂岩等等， 这使地 震波的速度也具有混合过渡的性质， 并取决于各个组分的含量。 即使在沉 积相对稳定的海相地层中， 仍然存在成份的差异， 只是相对于陆相地层来 说， 变化很小而已。 此外孔隙度、 孔隙充填物类型及孔隙流体压力等也会 对速度产生影响。 如果仍然釆用一次函数来描述地震速度的横向变化， 我 们可以发现即使在一个相对较小的工区内 6值的变化也是很大的。
有上述讨论我们可以得出一个基本的结论： 速度随深度的变化可以釆 用一次函数来描述， 并且 H-V曲线在横向上是变化的， 但曲线中的 a值和 b值变化趋势不同， a值是相对稳定的， 在横向上变化较小； b的影响因素 较多， 在横向上变化也较为剧烈。
应该说上述基本结论是得到普遍认可的， 问题的关键在于如何将上述 基本结论应用于实际的复杂构造速度构建中， 也就是在实际的速度构建中 如何构建横向可变的 H-V曲线。尤其是复杂构造速度构建结果往往要用于 实际的时深转换， 时深转换的结果要与实际的钻井结果相吻合， 因此必须 采用钻井结果约束横向可变的 H-V曲线的构建。为此我们提出了一种井约 束的横向可变 H-V曲线构建方法， 其基本过程如下-
A、 用声波测井曲线计算出工区内 （工区指地震勘探的工作区域） 各 井各层 （石油钻井的各个地质层位） 的层速度， 用实际的钻井结果 （地层 的实际深度）计算出各井各层（石油钻井的各个地质层位）的换算层速度；
B、对于 A步骤中的每一个层位、每一口井的井点的坐标， 以该点（井 的坐标） 为圆心， 以一定的 （需要由操作人员给定， 如 3000米等） 半径画 一个圆， 将该圆内的所有钻井点集中起来形成一个集合， 若集合中井的数 量小于预先给定的一个值即最小井数量， 则扩大半径， 直到集合中井的数 量大于或等于预先给定的一个值为止； （通常为 4 )
C、 根据集合中各井的测井层速度， 采用回归分析法拟合速度随深度 变化的一次函数， 将该函数的 a值 （根据公式 v=ah+b确定）， 作为该井点 的 a值；
D、 根据该井点的 a值和由 A步骤中得出该井点的换算层速度， 根据 公式 v=ah+b计算出 b值， 作为该井点的 b值；
E、 用克里金插值算法， 分别对 C步骤中得到的 a值和 D步骤中得到 的 b值进行网格化 （网格化是指采用一定的网格化方法,即插值）， 得到整 个工区内每一个井点的 a值和 b值， 这样工区内每个点都有了自己的 H-V 线性关系 v=ah+b， 由于每个点上的 a值和 b值不同， 我们称为横向可变 H-V曲线， 根据公式 v=ah+b形成井约束横向可变 H-V曲线。
3 . 横向可变 H-V曲线约束的速度反演算法
在实际的时深转换中， 我们已知的是地震反射时间， 并不知道反射深 度， 因此上述方法并不能直接引用。 我们必须将地震反射时间转换为反射 深度， 但反射深度正是我们时深转换的目标， 为此我们提出横向可变 H-V 曲线约束的速度反演算法， 方法的思路如下：
( 1 ) 应用对比解释的层位模型和 H-V曲线， 按 CDP点用逐层剥离的 方法， 由浅向深逐层递推该点每一层初始深度；
( 2 ) 应用 H-V曲线函 = 迭代每层厚度 (其中 i为 CDP序 号、 j为层位序号）， 即由地震反射时间和 H-V曲线函数计算深度， 根据计 算深度与 Η-ν量板深度的差修改速度， 这一过程进行到深度误差满足精度 为止， 迭代收敛时就取得该点每层的层速度和底界深度。
( 3 )对每个 CDP点做相同处理， 最终可取得每个 CDP点对应的每层 的层速度和底界深度。
对于第 i个 CDP的第 j层， 其 H-v曲线函数为 ^.) - /^) , 假设 j-i 层的深度 H,.,^已知（第 0层的深度固定为零）， 计算第 j层深度和速度的反 演算法如下- 给定第 j层的初始厚度 根据第 j 层的 H-V 曲线函数 = 计算 j 层的初始层速度 vi/ = /( ); 计算第 j层与 j-l层的理论反射时间差 At = hOy /Vy； 设第 j层与 j-1层的实际反射时间差为 ΔΓ (可由层位届时获得的层位 时间计算得到， 若理论反射时间与实际反射时间的差 |At- ΔΓ|<£ 为预 先给定的一个较小的值， 通常取釆样率的 1/2)， 这时的 就是第 j层的实 际厚度， ^就是第 j层的实际层速度。 算法结束， 输出反演结果。 若 At- ΔΓ<0， 表明厚度 /.偏小， 令 = + |Δ - ΔΓ|·ι¾.， 转第（2)步； 若 Δί-ΔΓ>0, 表明厚度 偏大， 令 hu=hir\ T\'Vij， 转第 （2) 步； 为提高反演算法的计算速度，可令第 j层的初始厚度 /.为第 i-1个 CDP 的第 j层的速度， 理论计算和实际应用的结果表明： 该算法的收敛性很好， 在我们的应用中尚未发现不收敛的情况。 这种方法不仅适合平缓构造的时 深转换， 同时也适合髙陡复杂构造区的时深转换， 同时递推反演算法能对 断层 (正或逆)作出判断并处理， 完全由断层附近解释的层位模型控制， 不 受人为因素影响， 因此可客观恢复断层下盘形态。
4. 地震波速度场构建中井约束横向可变 H-V曲线构建方法
实现地震波速度场构建中井约束横向可变 H-V 曲线构建方法的关键 是要融合地震层速度和钻井换算速度， 即实现不同类型速度的融合。 我们 可以这样来考虑这个问题， 如果消除了地震层速度和钻井换算速度中深度 (压实） 对速度的影响， 引起速度横向变化的主要因素就是沉积相带的横 向变化， 反映在深度一速度关系曲线中就是 &值的变化， 而沉积相带的横 向变化是缓慢的， 因而 &值的变化也应该是缓慢的。如果在地震层速度的 值与钻井换算速度的 δ值之间建立一个比例关系，用一个比例系数来描述， 那么这个比例系数也应该是缓变的。 当然这个比例系数受各种干扰因素的 影响， 在个别部位有可能出现突变， 但总体趋势应该是缓变的。 对于复杂 构造的速度模型构建问题来说，孔隙、裂缝及充填物等多种因素引起的 &值 我们可以忽略不计， 只考虑沉积相带的横向变化引起的速度横向变化。 由 此我们得出进行地震波速度场构建的基本思路如下- 根据多年的工作经验我们认为速度随深度的变化可以采用一次函数来 描述， 即认为速度随深度线性增加的， 可以采用一次函数消除深度（压实） 对速度的影响；
a、 基于已构建的井约束横向可变 H-V曲线， 基于已建立的井约束横 向可变 H-V曲线， 在钻井位置计算地震层速度的 b值与钻井换算速度的 b 值之间的比例关系， 即： 比例系数-钻井换算速度的 b值 /地震层速度的 b 值；
b、 在横向上 （横向是指平面上延伸的每一个点） 通过克里金插值获 得整个工区的比例系数；
c、用整个工区的地震层速度 b值与所述比例系数相乘， 得到整个工区 地震与多井约束的 b值；
d、将 c步骤中得到的 b值应用于上述公式 v=ah+b中，在每一个点上， 根据其 a值、 b值以及深度数据， 可获得每个点的速度， 由此就构建 Γ地 震波速度场。 '

Claims

权 利 要 求 书
1、 一种地震波速度场构建中井约束横向可变 H-V 曲线构建方法， 速 度-深度即 H-V 曲线中， 速度随深度的变化釆用一次函数来描述， 且速度 随深度线性增加，用公式 v=ah+b表示，其中 V为地层速度， h为地层埋深， a, b为速度-深度曲线的特征参数； 其特征在于： 所述速度-深度曲线为井 约束横向可变 H-V曲线， 井约束横向可变 H-V曲线的构建步骤如下-
A、 用声波测井曲线计算出工区内各井各层的层速度， 用实际的钻井 结果计算出各井各层的换算层速度， 所述工区是指地震勘探的工作区域， 所述各井各层是指石油钻井的各个地质层位， 所述实际的钻井结果是指地 层的实际深度；
B、 对于 A步骤中的每一个层位、 每一口井的井点的坐标， 以该坐标 为圆心画一个圆， 将该圆内的所有钻井点集中起来形成一个集合， 若集合 中井的数量小于预先给定的一个值即最小井数量， 则扩大半径， 直到集合 中井的数量大于或等于预先给定的一个值为止， 所述预先给定的一个值为 4；
C、 根据 B步骤中所述集合中各井的测井层速度， 采用回归分析法拟 合速度随深度变化的一次函数， 将该函数的 a值， 作为该井点的 a值；
D、 根据该井点的 a值和由 A步骤中得出该井点的换算层速度， 根据 公式 v=ah+b计算出 b值， 作为该井点的 b值；
E、 用克里金插值算法， 分别对 C步骤中得到的 a值和 D步骤中得到 的 b值进行网格化， 得到整个工区内每一个井点的 a值和 b值， 根据公式 v=ah+b形成井约束横向可变 H-V曲线。
2、 根据权利要求 1所述的地震波速度场构建中井约束横向可变 H-V 曲线构建方法， 其特征在于： 将采用本发明构建的井约束横向可变 H-V曲 线应用于地震波速度场构建的方法如下：
a、 基于已构建的井约束横向可变 H-V曲线， 基于已建立的井约束横 向可变 H-V曲线， 在钻井位置计算地震层速度的 b值与钻井换算速度的 b 值之间的比例关系， 即： 比例系数=钻井换算速度的 b值 /地震层速度的 b 值；
b、 在横向上通过克里金插值获得整个工区的比例系数， 所述横向是 指平面上延伸的每一个点；
c、用整个工区的地震层速度 b值与所述比例系数相乘，得到整个工区 地震与多井约束的 b值；
d、将 c步骤中得到的 b值应用于上述公式 v=ah+b中，在每一个点上， 根据其 a值、 b值以及深度数据， 可获得每个点的速度， 由此就构建了地 震波速度场。 π 3、 根据权利要求 1所述的地震波速度场构建中井约束横向可变 H-V 曲线构建方法， 其特征在于： 所述的井约束横向可变 H-V曲线是经过速度 反演算法处理的， 其处理步骤如下：
(1) 应用对比解释的层位模型和 H-V曲线， 按 CDP点用逐层剥离的 方法， 由浅向深逐层递推该点每一层初始深度；
(2) 应用 H-V曲线函数^(¾ =/( ^^迭代每层厚度， 其中 i为 CDP 序号、 j为层位序号， 即由地震反射时间和 H-V曲线函数计算深度， 根据 计算深度与 H-V量板深度的差修改速度，这一过程进行到深度误差满足精 度为止， 迭代收敛时就取得该点每层的层速度和底界深度；
(3)对每个 CDP点做相同处理， 最终可取得每个 CDP点对应的每层 的层速度和底界深度。 '
PCT/CN2010/000145 2009-04-13 2010-02-02 地震波速度场构建中井约束横向可变h-v曲线构建方法 WO2010118624A1 (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
EP10764037.7A EP2420862A4 (en) 2009-04-13 2010-02-02 Well constrained horizontal variable h-v curve constructing method for seismic wave velocity field construction
EA201101493A EA020635B1 (ru) 2009-04-13 2010-02-02 Способ сейсмической разведки для поиска углеводородов с применением поля средних скоростей сейсмической волны
US13/273,096 US8868348B2 (en) 2009-04-13 2011-10-13 Well constrained horizontal variable H-V curve constructing method for seismic wave velocity field construction

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CN2009100589208A CN101533103B (zh) 2009-04-13 2009-04-13 地震波速度场构建中井约束横向可变h-v曲线构建方法
CN200910058920.8 2009-04-13

Related Child Applications (1)

Application Number Title Priority Date Filing Date
US13/273,096 Continuation US8868348B2 (en) 2009-04-13 2011-10-13 Well constrained horizontal variable H-V curve constructing method for seismic wave velocity field construction

Publications (1)

Publication Number Publication Date
WO2010118624A1 true WO2010118624A1 (zh) 2010-10-21

Family

ID=41103824

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CN2010/000145 WO2010118624A1 (zh) 2009-04-13 2010-02-02 地震波速度场构建中井约束横向可变h-v曲线构建方法

Country Status (5)

Country Link
US (1) US8868348B2 (zh)
EP (1) EP2420862A4 (zh)
CN (1) CN101533103B (zh)
EA (1) EA020635B1 (zh)
WO (1) WO2010118624A1 (zh)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103790569A (zh) * 2012-10-26 2014-05-14 中国石油化工股份有限公司 一种利用vsp时深关系对声波测井资料进行校正的方法
CN110618460A (zh) * 2019-07-22 2019-12-27 中国石油化工股份有限公司 一种结合层位信息的微测井方位加权插值建模方法
CN111856585A (zh) * 2020-07-01 2020-10-30 中国石油天然气集团有限公司 时深转换速度场的构建方法、装置、设备及可读存储介质
CN112255686A (zh) * 2020-10-14 2021-01-22 东北石油大学 一种基于回归算法的断层边部速度建模方法
CN116381782A (zh) * 2023-02-17 2023-07-04 中海石油(中国)有限公司深圳分公司 一种低幅度构造成图方法、设备及存储设备

Families Citing this family (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101533103B (zh) * 2009-04-13 2012-03-21 中国石油天然气集团公司 地震波速度场构建中井约束横向可变h-v曲线构建方法
GB2479172B (en) * 2010-03-31 2016-02-10 Statoil Petroleum As Estimating interval velocities
US10620330B2 (en) 2013-03-19 2020-04-14 Westerngeco L.L.C. Estimating translational data
WO2015137932A1 (en) * 2014-03-11 2015-09-17 Landmark Graphics Corporation Validation of depth-depth curves using time-seismic depth interval velocity
CN106154329B (zh) * 2015-04-16 2018-08-31 中国石油化工股份有限公司 利用vsp速度校正声波均方根速度的方法及其系统
CN105353412B (zh) * 2015-12-14 2017-08-29 中国石油大学(华东) 一种井震联合平均速度场的计算方法及系统
CN108427141B (zh) * 2017-02-15 2019-08-30 中国石油化工股份有限公司 一种在沉积地层中识别与提取周期性波动的方法及系统
CN108562938B (zh) 2018-03-23 2019-09-06 中国石油天然气股份有限公司 一种消除频散效应的方法、装置及系统
CN108828666B (zh) * 2018-08-07 2019-12-27 中国石油大学(华东) 一种正断层附近地层埋藏史的恢复方法
CN111596348B (zh) * 2019-02-21 2023-08-22 中国石油天然气集团有限公司 基于声波时差测井资料的地层层速度校正方法及装置
CN111856567B (zh) * 2019-04-30 2023-09-26 中国石油天然气集团有限公司 基于地震叠加速度的压实系数确定方法及装置
CN112649876A (zh) * 2019-10-12 2021-04-13 中国石油化工股份有限公司 建立地震偏移速度模型的方法和装置
CN112782756A (zh) * 2019-11-08 2021-05-11 中国石油天然气集团有限公司 一种基于自适应构造约束的约束层速度反演方法及系统
CN113589367B (zh) * 2020-04-30 2023-04-28 中国石油化工股份有限公司 基于面积守恒大断裂附近构造趋势矫正方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6502038B1 (en) * 2000-06-30 2002-12-31 Exxonmobil Upstream Research Company Method for non-hyperbolic moveout analysis of seismic data
US6864890B2 (en) * 2002-08-27 2005-03-08 Comoco Phillips Company Method of building and updating an anisotropic velocity model for depth imaging of seismic data
US20070247973A1 (en) * 2006-04-21 2007-10-25 Prism Seismic Inc. Method for converting seismic data from the time domain to the depth domain
WO2009005940A2 (en) * 2007-06-29 2009-01-08 Schlumberger Canada Limited Seismic inversion of data containing surface related multiples
CN101533103A (zh) * 2009-04-13 2009-09-16 中国石油集团川庆钻探工程有限公司 地震波速度场构建中井约束横向可变h-v曲线构建方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5136551A (en) * 1989-03-23 1992-08-04 Armitage Kenneth R L System for evaluation of velocities of acoustical energy of sedimentary rocks
US6185512B1 (en) * 1998-10-13 2001-02-06 Raytheon Company Method and system for enhancing the accuracy of measurements of a physical quantity
FR2856710B1 (fr) * 2003-06-30 2007-11-23 Michel Fonfrede Mur-barrage hydraulique a debit controle
WO2005010797A2 (en) * 2003-07-23 2005-02-03 Lee Wook B Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration
US20060047429A1 (en) * 2004-08-24 2006-03-02 Adams Steven L Method of estimating geological formation depths by converting interpreted seismic horizons from the time domain to the depth domain
CN100487488C (zh) * 2006-05-11 2009-05-13 中国石油集团东方地球物理勘探有限责任公司 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法
CN101329407B (zh) * 2007-06-20 2011-01-12 中国石油天然气集团公司 一种快速p-sv转换波直接模拟以确定地层岩性、岩相变化的方法
CN101354444B (zh) * 2007-07-25 2011-02-09 中国石油天然气集团公司 一种确定地层岩性和孔隙流体的方法
CN101363916B (zh) * 2007-08-10 2012-05-30 中国石油集团东方地球物理勘探有限责任公司 一种高精度折射静校正数据反演方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6502038B1 (en) * 2000-06-30 2002-12-31 Exxonmobil Upstream Research Company Method for non-hyperbolic moveout analysis of seismic data
US6864890B2 (en) * 2002-08-27 2005-03-08 Comoco Phillips Company Method of building and updating an anisotropic velocity model for depth imaging of seismic data
US20070247973A1 (en) * 2006-04-21 2007-10-25 Prism Seismic Inc. Method for converting seismic data from the time domain to the depth domain
WO2009005940A2 (en) * 2007-06-29 2009-01-08 Schlumberger Canada Limited Seismic inversion of data containing surface related multiples
CN101533103A (zh) * 2009-04-13 2009-09-16 中国石油集团川庆钻探工程有限公司 地震波速度场构建中井约束横向可变h-v曲线构建方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
YUAN, JINGJU: "Application of Layer-cake Method to Make Structure Map and Geological Base", GEOPHYSICAL PROSPECTING FOR PETROLEUM, vol. 45, no. 3, May 2006 (2006-05-01), pages 285 - 289, XP008168436 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103790569A (zh) * 2012-10-26 2014-05-14 中国石油化工股份有限公司 一种利用vsp时深关系对声波测井资料进行校正的方法
CN103790569B (zh) * 2012-10-26 2016-12-21 中国石油化工股份有限公司 一种利用vsp时深关系对声波测井资料进行校正的方法
CN110618460A (zh) * 2019-07-22 2019-12-27 中国石油化工股份有限公司 一种结合层位信息的微测井方位加权插值建模方法
CN111856585A (zh) * 2020-07-01 2020-10-30 中国石油天然气集团有限公司 时深转换速度场的构建方法、装置、设备及可读存储介质
CN112255686A (zh) * 2020-10-14 2021-01-22 东北石油大学 一种基于回归算法的断层边部速度建模方法
CN112255686B (zh) * 2020-10-14 2023-08-11 东北石油大学 一种基于回归算法的断层边部速度建模方法
CN116381782A (zh) * 2023-02-17 2023-07-04 中海石油(中国)有限公司深圳分公司 一种低幅度构造成图方法、设备及存储设备

Also Published As

Publication number Publication date
US8868348B2 (en) 2014-10-21
EA201101493A1 (ru) 2012-03-30
EP2420862A1 (en) 2012-02-22
EP2420862A4 (en) 2017-05-03
CN101533103B (zh) 2012-03-21
CN101533103A (zh) 2009-09-16
EA020635B1 (ru) 2014-12-30
US20120053841A1 (en) 2012-03-01

Similar Documents

Publication Publication Date Title
WO2010118624A1 (zh) 地震波速度场构建中井约束横向可变h-v曲线构建方法
EP3028071B1 (en) Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
CN109738945B (zh) 一种利用叠前深度偏移成果直接生成构造图的方法
CN107783187B (zh) 一种将测井速度和地震速度结合建立三维速度场的方法
GB2490278A (en) Multicomponent seismic inversion of VSP data
CN107817526B (zh) 叠前地震道集分段式振幅能量补偿方法及系统
WO2017035104A1 (en) Velocity model seismic static correction
Arnulf et al. Nature of upper crust beneath the Lucky Strike volcano using elastic full waveform inversion of streamer data
CN104516018A (zh) 一种地球物理勘探中岩性约束下的孔隙度反演方法
CN109521474B (zh) 一种三维双控下的叠前地质统计学反演方法
CN111722284B (zh) 一种基于道集数据建立速度深度模型的方法
EP2994781A1 (en) Estimation of q-factor in time domain
CN113031068B (zh) 一种基于反射系数精确式的基追踪叠前地震反演方法
CN108398720A (zh) 一种基于杨氏模量、泊松比的两项式地震叠前反演方法
AU2011318531A1 (en) System and method for characterization with non-unique solutions of anisotropic velocities
CN111077578B (zh) 岩层分布预测方法和装置
CN109143335B (zh) 一种合成地震记录的制作方法、系统、介质及设备
CN109459790B (zh) 针对煤系地层地震速度场建立方法及系统
Petukhin et al. Construction and waveform testing of the crustal and basin structure models for southwest Japan
CN104111474B (zh) 一种基于线性混合的叠前弹性参数反演实现的储层与流体预测方法
CN109839659B (zh) 一种对叠前深度偏移剖面进行迭代优化的方法
Zou et al. Low-frequency impedance modeling based on well logs and seismic attributes
Matheney et al. Seismic attribute inversion for velocity and attenuation structure using data from the GLIMPCE Lake Superior experiment
CN117195511B (zh) 一种初始地壳厚度和伸展系数定量计算方法
CN111596348B (zh) 基于声波时差测井资料的地层层速度校正方法及装置

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 10764037

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: 2010764037

Country of ref document: EP

WWE Wipo information: entry into national phase

Ref document number: 201101493

Country of ref document: EA

WWE Wipo information: entry into national phase

Ref document number: 2400/DELNP/2012

Country of ref document: IN