CN105572736B - 一种利用地震资料定量预测坡折带的方法 - Google Patents
一种利用地震资料定量预测坡折带的方法 Download PDFInfo
- Publication number
- CN105572736B CN105572736B CN201610008703.8A CN201610008703A CN105572736B CN 105572736 B CN105572736 B CN 105572736B CN 201610008703 A CN201610008703 A CN 201610008703A CN 105572736 B CN105572736 B CN 105572736B
- Authority
- CN
- China
- Prior art keywords
- msub
- mrow
- gradient
- slope
- break belt
- Prior art date
- Legal status (The legal status 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 status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 230000000694 effects Effects 0.000 abstract description 6
- 238000010276 construction Methods 0.000 description 7
- 230000000875 corresponding effect Effects 0.000 description 6
- 208000035126 Facies Diseases 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000000151 deposition Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 239000004576 sand Substances 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000593 degrading effect Effects 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 230000000994 depressogenic effect Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 239000004744 fabric Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000035772 mutation Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
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)
- Pit Excavations, Shoring, Fill Or Stabilisation Of Slopes (AREA)
- Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)
Abstract
本发明涉及一种利用地震资料定量预测坡折带的方法,A、根据目地质、地震等背景资料,对发育坡折带的层系进行精细三维地震层位解释,得到目的层层位数据;B、求取出目的层的坡度与坡向值;C、根据坡度、坡向平面图上的分布特征初步识别坡折带分布范围;若坡折带边界清晰,则直接求取坡折带几何参数:坡折、坡脚、坡度、高差、水平宽度;若坡折带边界模糊(如断裂影响),则进行D步骤;D、统计并求取已初步识别的坡折带范围内所有坡向值的平均值θ;E、以θ为中心,δ为步长在θ±ε范围内扫描求取最优视坡度φ;F、根据步骤E求取的最优视坡度,识别坡折带边界,并分别求取折带几何参数。本发明实现坡折带坡度平面刻画;提高坡度计算精度。
Description
技术领域
本发明涉及石油勘探技术领域,特别涉及一种利用地震资料定量预测坡折带的方法。
背景技术
坡折带原是地貌学概念,指地形坡度突变的地带,不论在沉积盆地中还是在剥蚀区,都可能发育坡折带。坡折带主要由以下几部分组成:坡折(slope break):指地形由缓变陡的转折点连线,为上凸的脊;陡坡(slope):指其上下的地形坡度都迅速变缓的陡坡地带;坡脚(slope toe):指地形由陡变缓的转折点连线,为下凹的槽。林畅松等通过对半地堑盆地同沉积构造坡折带的研究,认为断陷湖盆中存在的构造坡折带制约着盆地可容空间的变化,控制了层序的发育、沉积体系域及砂体的分布;王英民等通过对准噶尔盆地大型坳陷湖盆坡折带类型的探讨,指出坡折带对层序、沉积和非构造圈闭的发育具有十分重要的控制作用;张善文等提出利用坡折带理论寻找隐蔽油气藏并在济阳坳陷陆相湖盆勘探中取得了显著成果。大量研究表明,坡折带制约着盆地可容纳空间的变化,控制着斜坡区沉积相带的分布和地层层序的发育,对砂体的沉积、横向尖灭和地层岩性圈闭的发育也具有关键控制作用。因此坡折带的识别和几何参数的刻画对斜坡区沉积相划分和地层岩性圈闭发育有利区及圈闭识别都具有重要指导意义。
现有的破折带识别方法概括起来可分为两类,一类是以剖面变化特征来识别,包括地震、地质、沉积相和层序地层等剖面的识别,主要表现为地层厚度的急剧增加以及地层倾角的突然变化;另一类是以坡折带平面变化特征来识别,包括地层厚度、沉积相展布、地层倾角等来识别,其平面变化特征与剖面变化相对应。但是现有方法的不足之处有两点,第一是坡折带的划分多采用人工划分方法,其划分方案影响因素较多,且为对坡折带的定性认识,坡折、坡脚识别不准确;第二是不能刻画坡折带内部及其上下坡度的平面变化。
发明内容
本发明解决上述技术问题所采用的技术方案是:一种利用地震资料定量预测坡折带的方法,包括以下步骤:
A、根据目地质、地震背景资料,对发育坡折带的层系进行精细三维地震层位解释,得到目的层层位数据;
B、根据精细解释的地震层位数据和坡度计算模型求取出目的层的坡度与坡向值;
C、根据坡度、坡向平面图上的分布特征初步识别坡折带分布范围;若坡折带边界清晰,则直接求取坡折带几何参数:坡折、坡脚、坡度、高差、水平宽度;若坡折带边界模糊(如断裂影响),则进行D步骤;
D、统计并求取已初步识别的坡折带范围内所有坡向值的平均值θ;
E、以θ为中心,δ为步长在θ±ε范围内扫描求取最优视坡度φ;
F、根据步骤E求取的最优视坡度,识别坡折带边界,并分别求取折带几何参数。
进一步的,所述步骤E中的最优视坡度φ为已初步识别的坡折带范围内∑ (βt-βr)2最小的视坡度,其中βt为真坡度,βr为视坡度,原理是保持坡折带范围内真坡度减小的最少。
进一步的,所述步骤B中坡度与坡向值的求取公式如下:
其中β为坡度,α为坡向,fx为x方向的偏导数,fy为y方向的偏导数, zi为5×5局部移动窗口内对应的属性值(如高程,地层深度等),w为5×5局部窗口网格间隔,i为5×5局部窗口网格点编号。
本发明的有益效果:(1)实现坡折带从定性识别到定量刻画,对坡折带坡折、坡脚、坡度、水平宽度、高差等几何参数进行定量刻画;(2)实现坡折带坡度平面刻画,对坡折带斜坡内坡度、坡折之上、坡脚之下的坡度进行平面展布刻画;(3)提高坡度计算精度,减小断裂、局部构造变形等非坡折因素的影响,提高坡折带识别和参数刻画精度。
附图说明
图1是坡度、坡向示意图;
图2是5×5局部移动窗口及权值分布图;
图3是实施例中目的层坡折带地震剖面;
图4是实施例中目的层坡度图;
图5是实施例中目的层坡向图;
图6是实施例中最优视坡度图;
图7是实施例中坡折带与断裂系统叠合图;
图8是坡度带参数定量刻画流程图。
具体实施方式
本发明的一种利用地震资料定量预测坡折带的方法,包括以下步骤:A、根据目地质、地震背景资料,对发育坡折带的层系进行精细三维地震层位解释,得到目的层层位数据;B、根据精细解释的地震层位数据和坡度计算模型求取出目的层的坡度与坡向值;C、根据坡度、坡向平面图上的分布特征初步识别坡折带分布范围;若坡折带边界清晰,则直接求取坡折带几何参数:坡折、坡脚、坡度、高差、水平宽度;若坡折带边界模糊(如断裂影响),则进行D步骤;D、统计并求取已初步识别的坡折带范围内所有坡向值的平均值θ;E、以θ为中心,δ为步长在θ±ε范围内扫描求取最优视坡度φ;F、根据步骤E求取的最优视坡度,识别坡折带边界,并分别求取折带几何参数。
优选的实施方案,所述步骤E中的最优视坡度φ为已初步识别的坡折带范围内∑(βt-βr)2最小的视坡度,其中βt为真坡度,βr为视坡度,原理是保持坡折带范围内真坡度减小的最少。
优选的实施方案,所述步骤B中坡度与坡向值的求取公式如下:
其中β为坡度,α为坡向,fx为x方向的偏导数,fy为y方向的偏导数, zi为5×5局部移动窗口内对应的属性值(如高程,地层深度等),w为5×5局部窗口网格间隔,i为5×5局部窗口网格点编号。
上述步骤的流程为图8所示,且原理是利用了坡度计算模型;其中坡度模型中空间曲面的坡度和坡向是互相联系的两个点位函数,坡度反映曲面的倾斜程度;坡向是斜坡面对的方向,一般以北方向为起始方向,并按顺时针方向度量(如图1所示),地表某点的坡度β、坡向α是地形曲面z=f(x,y)分别沿地图坐标x、y方向求取偏导fx、fy的函数,关系式为:
上述公式(1)和(2)中,对偏导数fx、fy的计算是准确求取坡度的关键。 Evan(1980)首次提出局部曲面拟合法,即在规则网格的局部对称窗口(如3×3, 5×5窗口)利用二次曲面函数拟合地表高程,通过函数系数与函数偏导数对应关系求取fx、fy;Zevenbergen和Thornev(1987)等用类似的方法在3×3局部窗口中提出利用不完全四次曲面拟合;Florinsky(2009)在5×5局部窗口提出利用三次曲面拟合,拟合函数:
其中,函数的各项系数与函数偏导数存在以下关系:
在5×5局部窗口(图2)中,设定各点坐标分别为:(-2w,2w,z1),(-w, 2w,z2),(0,2w,z3),(w,2w,z4),(2w,2w,z5),(-2w,w,z6),(-w,w,z7), (0,w,z8),(w,w,z9),(2w,w,z10),(-2w,0,z11),(-w,0,z12),(0,0, z13),(w,0,z14),(2w,0,z15),(-2w,-w,z16),(-w,-w,z17),(0,-w,z18), (w,-w,z19),(2w,-w,z20),(-2w,-2w,z21),(-w,-2w,z22),(0,-2w,z23), (w,-2w,z24)和(2w,-2w,z25)。将各个点坐标带入方程可求解各系数,由于方程个数大于未知数个数,因此采用最小二乘法求解各系数。根据最小二乘法原理,可以得到25个线性方程:
α=Fβ (5)
其中,α为25个已知高程z的列矩阵:
β为包含10个未知数的列向量:
F为已知坐标带入方程后与未知数β、已知数α对应的25×10的常数矩阵:
考虑不同距离的格网点对格网中心点参数提取的不同影响,即权的概念,常采用如下形式的权函数:mi为对应格网点的权值,di为该格网点离中心点的距离,n一般取1或者2。以n=2为例,权矩阵M为25×25 的方阵:
则(5)式可化为:
α=FTMβ (10)
可解得系数β:
β=(FTMF)-1FTMα (11)
以及不含权时,系数β:
β=(FTF)-1FTα (12)
由于计算坡度只需要偏导数因此这里只给出fx、fy表达式,采用反距离平方权:
上述优选的方案是,所述步骤D中的最优视坡度φ为最能反映坡折带边界的视坡度。在实际地质条件中,地层构造往往包含一定的方向性信息,即构造特征在一定范围内是沿着某个方向分布或延伸的,为了突出某一特定方向上的构造信息,在计算地层坡度时引入任意方向的视坡度βr:
βr=βcos(α-θ) (15)
式中,θ为选定的特定角度,是与正北方向顺时针的夹角。定义最能反映地质工作者(或地震解释人员)感兴趣的构造或地质体的视坡度为该构造或体质体的最优视坡度φ,一般情况下,最优视坡度与该构造或地质体的坡向有关,对于单一倾向的构造或地质体,可通过下式确定:
式中,βi、αi分别为该构造或地质体范围内每一点的真坡度和坡向,n为该范围内所有格网点;当Ω取得最小值时,此时的θ为最优视坡度φ。
上述坡度模型的精度分析:由于用实际资料进行精度分析没有明确的可对比对象或真值,因此采用一种与实际地层构造相似的数学曲面进行分析,该数学曲面的坡度真值可通过表达式直接求取,与利用模型求取的坡度对比可进行精度、误差分布等分析。该曲面函数表达式:
将曲面按不同间隔(表1)离散并利用该算法求取坡度,精度指标采用中误差:
式中,βiF为利用(13)、(14)式求取的坡度,βiτ为曲面真坡度,n为样点总数。
表1不同格网间隔中误差
下面结合附图和实施例对发明做进一步详细的说明。
实施例中的研究区是一个形成于晚海西期并长期继承发育的古凸起,整体隆升不均衡;图3该研究区L880测线地震剖面,目的层为北西高、南东低以斜坡形式过渡大型不整合面,其上大部分地层逐层沿该不整合面上超尖灭,其下部分地层直接超覆于石炭系之上并伴随地层削蚀现象;受石炭系古凸起和深部大型逆掩断裂的控制,目的层发育挠曲坡折带。
(1)根据精细的地震解释层位求取目的层坡度和坡向,并根据坡度和坡向平面分布特征初步识别坡折带分布范围;
图4为目的层坡度,图5为目的层坡向;据图4所示,看出坡折带呈“S”型分布贯穿研究区,坡脚线十分清晰,但坡折带中部由于断裂发育,坡折线难以划分;坡折带的坡向分布与坡度分布有一定的相关性,易于区分其他构造形态,同时,坡度和坡向能不同程度的反应断裂的分布。
(2)统计坡折带坡向值平均值θ,θ=89°。
(3)以θ为中心,δ为步长在θ±ε范围内扫描求取最优视坡度,其中δ=3°,ε=15°,得φ=95°。
图6为最优视坡度,是根据(2)、(3)步求取最能反应坡折带坡折和坡脚的视坡度;图7是已识别的坡折带和断裂系统叠加图,对比图3可以看出,最优视坡度减小了坡折带附近断裂的坡度值,而其本身的坡度基本不变,这种坡度的减小与构造的坡向十分相关,与坡折带得坡向差异越大,坡度值减小越多。断裂坡度的减小使得坡折线变得十分清晰,提高了坡折带边界识别精度。同时,坡折带的分布与实际勘探井位十分吻合。
(4)根据以上步骤的成果,求取坡折带几何参数:坡折、坡脚、坡度、高差、水平宽度。根据坡折带形态,将坡折带分为三段,分别为坡折带北段、坡折带中段和坡折带南段,各参数如下表2:
表2坡折带参数表
。
Claims (2)
1.一种利用地震资料定量预测坡折带的方法,其特征在于,包括以下步骤:
A、根据目地质、地震背景资料,对发育坡折带的层系进行精细三维地震层位解释,得到目的层层位数据;
B、根据精细解释的地震层位数据和坡度计算模型求取出目的层的坡度与坡向值;
C、根据坡度、坡向平面图上的分布特征初步识别坡折带分布范围;若坡折带边界清晰,则直接求取坡折带几何参数:坡折、坡脚、坡度、高差、水平宽度;若坡折带边界模糊,则进行D步骤;
D、统计并求取已初步识别的坡折带范围内所有坡向值的平均值θ;
E、以θ为中心,δ为步长在θ±ε范围内扫描求取最优视坡度φ;
F、根据步骤E求取的最优视坡度,识别坡折带边界,并分别求取折带几何参数;
所述步骤E中的最优视坡度φ为已初步识别的坡折带范围内∑(βt-βr)2最小的视坡度,其中βt为真坡度,βr为视坡度,原理是保持坡折带范围内真坡度减小的最少;
其中ε=15°。
2.根据权利要求1所述的一种利用地震资料定量预测坡折带的方法,其特征在于,所述步骤B中坡度与坡向值的求取公式如下:
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>x</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>13968</mn>
<mi>w</mi>
</mrow>
</mfrac>
<mo>{</mo>
<mn>496</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>4</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>24</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>22</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>455</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>21</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>5</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>25</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1960</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>19</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>9</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>7</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>17</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>50</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>15</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>11</mn>
</msub>
<mo>+</mo>
<mn>88</mn>
<mo>(</mo>
<mrow>
<msub>
<mi>z</mi>
<mn>14</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>12</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>152</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>16</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>10</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>20</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mfenced open = "" close = "">
<mtable>
<mtr>
<mtd>
<mrow>
<msub>
<mi>f</mi>
<mi>y</mi>
</msub>
<mo>=</mo>
<mfrac>
<mn>1</mn>
<mrow>
<mn>13968</mn>
<mi>w</mi>
</mrow>
</mfrac>
<mo>{</mo>
<mn>496</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>6</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>10</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>16</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>20</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>455</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>21</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>25</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>1</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>5</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>1960</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>7</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>9</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>17</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>19</mn>
</msub>
<mo>)</mo>
</mrow>
</mrow>
</mtd>
</mtr>
<mtr>
<mtd>
<mrow>
<mo>+</mo>
<mn>50</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>3</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>23</mn>
</msub>
<mo>+</mo>
<mn>88</mn>
<mo>(</mo>
<mrow>
<msub>
<mi>z</mi>
<mn>8</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>18</mn>
</msub>
</mrow>
<mo>)</mo>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mn>152</mn>
<mrow>
<mo>(</mo>
<msub>
<mi>z</mi>
<mn>22</mn>
</msub>
<mo>+</mo>
<msub>
<mi>z</mi>
<mn>24</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>2</mn>
</msub>
<mo>-</mo>
<msub>
<mi>z</mi>
<mn>24</mn>
</msub>
<mo>)</mo>
</mrow>
<mo>}</mo>
</mrow>
</mtd>
</mtr>
</mtable>
</mfenced>
<mrow>
<mi>&beta;</mi>
<mo>=</mo>
<mi>arctan</mi>
<msqrt>
<mrow>
<msubsup>
<mi>f</mi>
<mi>x</mi>
<mn>2</mn>
</msubsup>
<mo>+</mo>
<msubsup>
<mi>f</mi>
<mi>y</mi>
<mn>2</mn>
</msubsup>
</mrow>
</msqrt>
</mrow>
其中β为坡度,α为坡向,fx为x方向的偏导数,fy为y方向的偏导数,zi为5×5局部移动窗口内对应的属性值,w为5×5局部窗口网格间隔,i为5×5局部窗口网格点编号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008703.8A CN105572736B (zh) | 2016-01-07 | 2016-01-07 | 一种利用地震资料定量预测坡折带的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610008703.8A CN105572736B (zh) | 2016-01-07 | 2016-01-07 | 一种利用地震资料定量预测坡折带的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105572736A CN105572736A (zh) | 2016-05-11 |
CN105572736B true CN105572736B (zh) | 2017-12-05 |
Family
ID=55883077
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610008703.8A Expired - Fee Related CN105572736B (zh) | 2016-01-07 | 2016-01-07 | 一种利用地震资料定量预测坡折带的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105572736B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106772586B (zh) * | 2017-02-20 | 2019-04-02 | 长江大学 | 一种基于地震信号奇异性的隐蔽性断裂检测方法 |
CN109979004A (zh) * | 2019-03-20 | 2019-07-05 | 华北理工大学 | 一种量化叠合盆地中帚状构造变形程度的方法和装置 |
CN110703325B (zh) * | 2019-09-11 | 2021-03-05 | 中国海洋石油集团有限公司 | 基于构造差异活动的烃源识别方法 |
CN111505720B (zh) * | 2020-05-21 | 2023-03-28 | 中国海洋石油集团有限公司 | 一种岩性圈闭的刻画方法 |
CN111983678B (zh) * | 2020-07-22 | 2024-02-23 | 中海石油(中国)有限公司深圳分公司 | 一种快速评价深水砂体发育潜力的方法 |
CN113933909B (zh) * | 2021-08-30 | 2024-03-19 | 核工业北京地质研究院 | 一种快速确定裂谷盆地边界的方法 |
CN116047598B (zh) * | 2022-11-07 | 2023-09-19 | 北京派特杰奥科技有限公司 | 一种基于人工智能的石油勘探方法及系统 |
CN116540302B (zh) * | 2023-07-04 | 2023-11-10 | 中海石油(中国)有限公司深圳分公司 | 一种定量评价源上圈闭油气成藏有效性的方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5258960A (en) * | 1992-11-27 | 1993-11-02 | Atlantic Richfield Company | Seismic velocity estimation method |
CN101403798B (zh) * | 2008-11-12 | 2011-01-26 | 重庆大学 | 一种滑坡主滑方向的识别方法 |
CN103364830B (zh) * | 2013-07-24 | 2016-07-06 | 北京师范大学 | 一种基于多因子的地震后崩滑灾害发生位置的预测方法 |
CN105137482B (zh) * | 2015-09-09 | 2017-12-15 | 中国海洋石油总公司 | 一种沉积体古坡度的计算方法 |
-
2016
- 2016-01-07 CN CN201610008703.8A patent/CN105572736B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105572736A (zh) | 2016-05-11 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105572736B (zh) | 一种利用地震资料定量预测坡折带的方法 | |
Sippl et al. | Geometry of the Pamir‐Hindu Kush intermediate‐depth earthquake zone from local seismic data | |
CN102053270B (zh) | 一种基于沉积地层单元的地震相分析方法 | |
Feng et al. | Interpretation of seismic reflection profiling data for the structure of the San Andreas fault zone | |
CN104614763B (zh) | 基于反射率法的多波avo储层弹性参数反演方法及系统 | |
CN101918864B (zh) | 地球物理数据处理系统 | |
CN107966732B (zh) | 基于空间结构导向的地震属性变化率求取方法 | |
CN109143361B (zh) | 一种基于层序地层学的碳酸盐岩地层古地质图的编制方法 | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
Husen et al. | Local and regional minimum 1D models for earthquake location and data quality assessment in complex tectonic regions: application to Switzerland | |
CN106646601A (zh) | 多信息联合约束的浅中深层三维q体建立方法 | |
CN106094030A (zh) | 一种通过地震剖面定量恢复湖盆最大古水深的方法 | |
CN102944896A (zh) | 表层调查数据的模型法静校正方法 | |
McClymont et al. | Characterization of the shallow structures of active fault zones using 3‐D ground‐penetrating radar data | |
CN107748399A (zh) | 利用重力界面反演识别山前带深部构造层方法 | |
CN110231652A (zh) | 一种基于密度的含噪声应用空间聚类的地震相提取方法 | |
Edwards et al. | Corrugated megathrust revealed offshore from Costa Rica | |
Robl et al. | The topography of a continental indenter: The interplay between crustal deformation, erosion, and base level changes in the eastern Southern Alps | |
CN102053263A (zh) | 调查表层结构的方法 | |
Hodge et al. | Evidence from high‐resolution topography for multiple earthquakes on high slip‐to‐length fault scarps: The Bilila‐Mtakataka fault, Malawi | |
Guglielmetti et al. | On the use of gravity data in delineating geologic features of interest for geothermal exploration in the Geneva Basin (Switzerland): prospects and limitations | |
CN107316341A (zh) | 一种多点地质统计学沉积相建模方法 | |
CN107797148B (zh) | 一种基于三维地质建模的航磁异常场分离方法及系统 | |
Saltogianni et al. | The 2014 Cephalonia earthquakes: finite fault modeling, fault segmentation, shear and thrusting at the NW Aegean Arc (Greece) | |
Fisher et al. | Determination of bedrock topography beneath the Greenland ice sheet by three‐dimensional imaging of radar sounding data |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20171205 Termination date: 20210107 |
|
CF01 | Termination of patent right due to non-payment of annual fee |