CN107942388A - 一种山区地表情况下的三角网格逆时偏移方法 - Google Patents
一种山区地表情况下的三角网格逆时偏移方法 Download PDFInfo
- Publication number
- CN107942388A CN107942388A CN201711103028.8A CN201711103028A CN107942388A CN 107942388 A CN107942388 A CN 107942388A CN 201711103028 A CN201711103028 A CN 201711103028A CN 107942388 A CN107942388 A CN 107942388A
- Authority
- CN
- China
- Prior art keywords
- mrow
- msub
- triangle
- computing unit
- polygon
- 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.)
- Pending
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. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
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)
Abstract
本发明涉及一种山区地表情况下的三角网格逆时偏移方法,其特征在于:模型剖分的过程为:将偏移涉及的地质区域利用三角网格进行剖分,三角网格之间无重叠区域和空白区域;剖分完成后,模型变为大量三角网格的集合体,每一个结点可能是多个三角形的共同顶点;具有共同顶点的多个三角网格构成一个计算单元,该计算单元为一个多边形;对计算单元中的每一个三角形取重心,将各三角形的重心连接起来形成一个新的多边形,该新的多边形为计算单元的对偶单元。
Description
技术领域
本发明涉及一种山区地表情况下的三角网格逆时偏移方法,属于油气勘探地震资料偏移成像技术领域。
背景技术
国内西部探区的地震资料品质低,一直是处理的难点。剧烈的地表高程变化和地表以下低降速带的变化使得地震资料明显不同于平原勘探采集的情况。常规的野外高程静校正方法利用小折射或者微测井实验得到近似的近地表速度,从而对地震数据做时差校正,消除高程变化对地震数据旅行时的影响。在高速岩体出露的复杂地表地形条件下,非规则的地表会引起强烈的地震波漫散射,从而使得地震记录信噪比降低,难以识别弱反射。接下来的地震数据处理流程就很受影响。因此,在起伏地表地形的情况下,如何不通过预先静校正就直接实施正演模拟和逆时偏移是很重要的一个研究课题。
逆时偏移(RTM)是基于双程波波动方程对地震波场进行沿时间方向的延拓。因此,RTM最大的优点就是没有对原始的波动方程进行近似,从而保证了波场延拓的精度。在1982年召开的52届SEG年会上,Whtemore等人首次提出了RTM的概念,但受限于当时的计算机计算能力和存储空间,该方法难以应用到实际生产中去。之后,McMechan,Gazdag,Baysal等做了深入研究,在理论方面和实际应用方面都有很多贡献,使其渐渐实用化。随着机群的出现和GPU的应用,机器计算能力大幅度提高,于是RTM自然推广到三维的情况。Hayashi等很多学者在如何减少RTM应用时的计算量和存储空间方面做了很多研究,受到工业界的重视。Valenciano等人在提高RTM波场延拓精度方面,做了很大贡献。在国内,RTM也日渐成为研究的热点,许多国内学者进行了研究,工业界也得到了一定的应用。尧德中、张文生、张关泉、张美根和王妙月等基于不同的算法实现了逆时偏移。但是基于山区地表情况下的RTM尚鲜见论著及应用。
发明内容
针对上述问题,本发明的目的是提供在山区地表地形情况下的三角网格高精度逆时偏移方法,该方法可直接将实际山区起伏地表勘探采集的数据作为输入,从而避开了静校正的难题,保证了成像精度。
为实现上述目的,本发明采用以下技术方案:一种山区地表情况下的三角网格逆时偏移方法,其特征在于:
模型剖分的过程为:将偏移涉及的地质区域利用三角网格进行剖分,三角网格之间无重叠区域和空白区域;剖分完成后,模型变为大量三角网格的集合体,每一个结点可能是多个三角形的共同顶点;具有共同顶点的多个三角网格构成一个计算单元,该计算单元为一个多边形;对计算单元中的每一个三角形取重心,将各三角形的重心连接起来形成一个新的多边形,该新的多边形为计算单元的对偶单元;
在该三角网格化的速度模型中,波场延拓公式表示为:
式中,hi表示计算单元中各三角形的共同顶点与计算单元多边形的第i个顶点的连线长度;si表示对偶单元多边形的第i个边的边长;N表示计算单元中三角网格的数量,Ai是包含公共顶点p0的对偶单元中被三角形网格的边所分割成的第i个小的多边形的面积;
式中,其余符号的含义如下:ui(j)表示第i个网格点处第j个时刻的波场;同理,u0(j-1),u0(j),u0(j+1)分别表示第0个网格点处第j-1、j、j+1时刻的波场;vi为与单元的公共顶点p0相邻的第i个网格点处的介质速度;Δt表示时间间隔。
进一步地,采用了混合边界消除人工边界产生的假反射,使得波场模拟接近实际野外地震勘探的情况。
本发明由于采取以上技术方案,其具有以下优点:1、本发明对速度模型采用三角网格划分,从而可以精确的刻画山区地表地形,以及地下的地层反射界面和奇异地质构造体,这是矩形网格难以做到的,即使对于简单的倾斜界面,矩形网格划分后会变为阶梯状或锯齿状,而网格划分的不精确会导致后续的正演偏移产生假散射等现象,严重影响成像的保真度和精度;此外,三角网格可以随着速度的分别灵活改变尺寸。2、本发明直接应用于山区复杂地表采集的地震数据;通常,对山区采集的地震资料进行处理,最大的瓶颈就是静校正难以做好,从而后续的各种处理都深受影响;而本发明的方法可以避开静校正的难题,直接从起伏地表进行波场延拓与成像,从而保证了精度。相比于有限元算法,本发明更为简洁实用,可操作性强;相比于有限差分算法,它对于起伏地表又有更好的灵活性和适应性;因此,其成像优势十分明显,尤其是在西部山区或黄土塬等复杂地区。
附图说明
图1是本发明对速度模型做三角网格剖分的示意图;
图2是对某油田提供的复杂地表模型做的偏移对比测试;其中,图(a)是山区起伏地表地震模型;图(b)是从地震数据中抽取的某一炮记录;图(c)是利用单程波方程做基于起伏地表延拓的叠前有限差分法偏移的结果;图(d)是本发明利用直接基于山区复杂地表的三角网格逆时偏移结果;
图3是国际通用的山区起伏地表地震模型的三角网格逆时偏移和频空域有限差分法单程波偏移处理结果对比;其中,图(a)是对起伏地表的波场利用直接下延的办法做频率空间域有限差分得到的单程波偏移结果;图(b)是基于山区地表地形的三角网格逆时偏移结果。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
本发明提出了山区地表情况下的三角网格逆时偏移方法,其特征在于:
模型剖分的方案为:将偏移涉及的地质区域利用三角网格进行剖分,三角网格之间无重叠区域和空白区域。剖分完成后,模型变为大量三角网格的集合体,每一个结点可能是多个三角形的共同顶点。具有共同顶点的多个三角网格构成一个计算单元,该计算单元为一个多边形。例如在图1中,具有共同顶点p0的五个三角形构成一个计算单元,即如图所示的五边形。对计算单元中的每一个三角形取重心,将各三角形的重心连接起来形成一个新的多边形,该新的多边形为计算单元的对偶单元。在图1中,某一个计算单元中的五个顶点用Pi表示,其对偶单元也是五边形,顶点用Qi表示。
在该三角网格化的速度模型中,波场延拓公式表示为:
式中,hi表示计算单元中各三角形的共同顶点与计算单元多边形的第i个顶点的连线长度;si表示对偶单元多边形的第i个边的边长;N表示计算单元中三角网格的数量,Ai是包含公共顶点p0的对偶单元中被三角形网格的边所分割成的第i个小的多边形的面积。
式中,其余符号的含义如下:ui(j)表示第i个网格点处第j个时刻的波场;同理,u0(j-1),u0(j),u0(j+1)分别表示第0个网格点处第j-1、j、j+1时刻的波场;vi为与单元的公共顶点p0相邻的第i个网格点处的介质速度;Δt表示时间间隔。
在实际应用中,采用了混合边界,消除了人工边界产生的假反射,从而使得波场模拟接近实际野外地震勘探的情况。
下面通过几个具体实施例来说明本发明的效果。
为了验证本方法的应用效果,本发明分别在某油田提供的复杂地表模型和国际通用的山区起伏地表地震模型上做了测试,取得了满意的效果。同时,也对结果做了分析,验证了本发明的优越性。
图2是某油田提供的复杂地表模型地震模拟数据的处理结果。图2(a)是三角网格RTM用到的地震速度模型,横向展布超过10000米,间隔为50米。纵向深度达到4000米。为了更好地描述薄层,深度间隔取4米。图2(b)是将检波器埋置在复杂山区地表下采集的某一炮的炮记录。总时间长度是2.5秒,采样的时间间隔是0.001秒。图2(d)是本发明得到的三角网格RTM结果。和图2(c)频率空间域的基于复杂地表波场延拓的单程波偏移方法相比,无论是山区地表产生的漫散射,还是中深层的复杂异常构造体,都得到了很好的成像。地表漫散射的影响得到了很好的收敛,提高了偏移的精度,成功避开了静校正的难题。
图3是对国际通用的山区起伏地表地震模型地震数据(Sam Gray and GaryMaclean,1995)的处理结果。该模型描述了加拿大地区特异的地质断层构造特征,成像非常困难。常规的生产处理方法对这样的问题无能为力。本发明直接对未做层析法或射线法静校正的数据做处理,成像效果明显优于其他方法。该加拿大模型横向范围为0-25000米,检波器间隔15米。地质模型深度变化范围为0-10000米,深度采样点为1000个。图3(a)是频率空间域的基于复杂地表波场延拓的单程波偏移方法的结果,图3(b)是本发明做的偏移结果。
可以看到,山地地区由于高程变化及其剧烈,而且近地表横向地质条件非均质性很强,从高速岩体、砂砾石变化到黄土层或沙漠。采集的野外数据中,反射绕射同相轴严重的扭曲,并被来自地表地形产生的漫散射及地下奇异体散射影响,导致地震资料品质差,信噪比低。此时,射线理论不再适用,常规的高程静校正、射线法层析静校正甚至模拟退火法重建近地表模型都难以应用。本发明直接基于复杂地表地形,从而避开了传统处理流程的瓶颈问题—静校正,所以偏移成像结果得到了大大改善,高陡倾角、翻转断层的陡翼都能够清晰地成像。
从图中可以看出,采用本专利所提出的基于山区地表地形的三角网格逆时偏移方法是山区地震勘探资料处理的有效手段,能够解决复杂地表成像问题。这在山区地区勘探、油气田开发中会发挥越来越大的作用。
上述各实施例仅用于说明本发明,其中方法的实施步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。
Claims (2)
1.一种山区地表情况下的三角网格逆时偏移方法,其特征在于:
模型剖分的过程为:将偏移涉及的地质区域利用三角网格进行剖分,三角网格之间无重叠区域和空白区域;剖分完成后,模型变为大量三角网格的集合体,每一个结点可能是多个三角形的共同顶点;具有共同顶点的多个三角网格构成一个计算单元,该计算单元为一个多边形;对计算单元中的每一个三角形取重心,将各三角形的重心连接起来形成一个新的多边形,该新的多边形为计算单元的对偶单元;
在该三角网格化的速度模型中,波场延拓公式表示为:
<mrow>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>-</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>=</mo>
<mn>2</mn>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>+</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
<mo>+</mo>
<mfrac>
<msup>
<mrow>
<mo>(</mo>
<mi>&Delta;</mi>
<mi>t</mi>
<mo>)</mo>
</mrow>
<mn>2</mn>
</msup>
<mrow>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<msub>
<mi>A</mi>
<mi>i</mi>
</msub>
<msubsup>
<mi>v</mi>
<mi>i</mi>
<mn>2</mn>
</msubsup>
</mfrac>
</mrow>
</mfrac>
<munderover>
<mo>&Sigma;</mo>
<mrow>
<mi>i</mi>
<mo>=</mo>
<mn>1</mn>
</mrow>
<mi>N</mi>
</munderover>
<mfrac>
<msub>
<mi>s</mi>
<mi>i</mi>
</msub>
<msub>
<mi>h</mi>
<mi>i</mi>
</msub>
</mfrac>
<mo>&lsqb;</mo>
<msub>
<mi>u</mi>
<mi>i</mi>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>-</mo>
<msub>
<mi>u</mi>
<mn>0</mn>
</msub>
<mrow>
<mo>(</mo>
<mi>j</mi>
<mo>)</mo>
</mrow>
<mo>&rsqb;</mo>
</mrow>
式中,hi表示计算单元中各三角形的共同顶点与计算单元多边形的第i个顶点的连线长度;si表示对偶单元多边形的第i个边的边长;N表示计算单元中三角网格的数量,Ai是包含公共顶点p0的对偶单元中被三角形网格的边所分割成的第i个小的多边形的面积;
式中,其余符号的含义如下:ui(j)表示第i个网格点处第j个时刻的波场;同理,u0(j-1),u0(j),u0(j+1)分别表示第0个网格点处第j-1、j、j+1时刻的波场;vi为与单元的公共顶点p0相邻的第i个网格点处的介质速度;Δt表示时间间隔。
2.如权利要求1所述的一种山区地表情况下的三角网格逆时偏移方法,其特征在于:采用了混合边界消除人工边界产生的假反射,使得波场模拟接近实际野外地震勘探的情况。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711103028.8A CN107942388A (zh) | 2017-11-10 | 2017-11-10 | 一种山区地表情况下的三角网格逆时偏移方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711103028.8A CN107942388A (zh) | 2017-11-10 | 2017-11-10 | 一种山区地表情况下的三角网格逆时偏移方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN107942388A true CN107942388A (zh) | 2018-04-20 |
Family
ID=61933745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711103028.8A Pending CN107942388A (zh) | 2017-11-10 | 2017-11-10 | 一种山区地表情况下的三角网格逆时偏移方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107942388A (zh) |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949421A (zh) * | 2019-03-26 | 2019-06-28 | 中煤航测遥感集团有限公司 | 三角网切割方法及装置 |
CN110376623A (zh) * | 2019-07-08 | 2019-10-25 | 中国空间技术研究院 | 星载gnss-r镜面反射点海洋潮汐修正定位方法和系统 |
CN113367664A (zh) * | 2021-06-26 | 2021-09-10 | 西北工业大学深圳研究院 | 一种基于人体浅表皮肤光声显微成像的声学补偿方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879820A (zh) * | 2012-09-20 | 2013-01-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于三角网格的三维表层模型构建方法 |
CN105386756A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种应用应变量计算脆性地层孔隙度的方法 |
-
2017
- 2017-11-10 CN CN201711103028.8A patent/CN107942388A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102879820A (zh) * | 2012-09-20 | 2013-01-16 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于三角网格的三维表层模型构建方法 |
CN105386756A (zh) * | 2014-09-04 | 2016-03-09 | 中国石油化工股份有限公司 | 一种应用应变量计算脆性地层孔隙度的方法 |
Non-Patent Citations (2)
Title |
---|
孙小东 等: "三角网格有限差分法叠前逆时偏移方法研究", 《地球物理学进展》 * |
孙小东 等: "基于三角网格的有限差分法叠前逆时偏移", 《SPG/SEG深圳2011国际地球物理会议论文集》 * |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109949421A (zh) * | 2019-03-26 | 2019-06-28 | 中煤航测遥感集团有限公司 | 三角网切割方法及装置 |
CN109949421B (zh) * | 2019-03-26 | 2023-06-16 | 中煤航测遥感集团有限公司 | 三角网切割方法及装置 |
CN110376623A (zh) * | 2019-07-08 | 2019-10-25 | 中国空间技术研究院 | 星载gnss-r镜面反射点海洋潮汐修正定位方法和系统 |
CN113367664A (zh) * | 2021-06-26 | 2021-09-10 | 西北工业大学深圳研究院 | 一种基于人体浅表皮肤光声显微成像的声学补偿方法及系统 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104570125B (zh) | 一种利用井数据提高成像速度模型精度的方法 | |
CN107526101B (zh) | 一种获取地震反射波的采集和处理方法 | |
CN103454685B (zh) | 利用测井约束波阻抗反演预测砂体厚度的方法和装置 | |
CN101329407B (zh) | 一种快速p-sv转换波直接模拟以确定地层岩性、岩相变化的方法 | |
Lee et al. | Development of a 3D velocity model of the Canterbury, New Zealand, region for broadband ground‐motion simulation | |
CN105388518A (zh) | 一种质心频率与频谱比联合的井中地震品质因子反演方法 | |
Pamuk et al. | Soil characterization of Bornova Plain (Izmir, Turkey) and its surroundings using a combined survey of MASW and ReMi methods and Nakamura’s (HVSR) technique | |
CN104570101A (zh) | 一种基于粒子群算法的avo三参数反演方法 | |
CN109188506A (zh) | 一种适用于高铁隧底地震ct的纯地表三维观测系统 | |
CN109870719A (zh) | 一种碳酸盐岩致密薄储层分布确定方法、装置及系统 | |
Barakat | Modern geophysical techniques for constructing a 3D geological model on the Nile Delta, Egypt | |
CN107765308A (zh) | 基于褶积思想与精确震源的重构低频数据频域全波形反演方法 | |
CN107942388A (zh) | 一种山区地表情况下的三角网格逆时偏移方法 | |
Chapman et al. | On the geologic structure at the epicenter of the 1886 Charleston, South Carolina, earthquake | |
CN104237946B (zh) | 基于井控的单层反射纵波和反射转换横波的振幅匹配方法 | |
CN104199088B (zh) | 一种提取入射角道集的方法及系统 | |
Gu et al. | Analysis of dynamic sedimentary environments in alluvial fans of some tributaries of the upper Yellow River of China based on ground penetrating radar (GPR) and sediment cores | |
Senna et al. | Modeling of the subsurface structure from the seismic bedrock to the ground surface for a broadband strong motion evaluation in Japan | |
CN102565852A (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
Bala | Modelling of seismic site amplification based on in situ geophysical measurements in Bucharest, Romania | |
CN105259577B (zh) | 一种确定地层界面的角度信息的方法及装置 | |
Parize et al. | A methodological approach to clastic injectites: From field analysis to seismic modeling—Examples of the Vocontian Aptian and Albian injectites (southeast France) | |
Zahoor et al. | Geophysical surveys in the Kashmir valley (J&K Himalayas) part I: estimation of dynamic parameters for soils and investigation of the deep basin structure | |
CN104076395A (zh) | 基于滤波组合的镜面能量提取与成像方法 | |
Dhakal et al. | Validation of the deep velocity structure of the Tokachi basin based on 3-D simulation of long-period ground motions |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180420 |