E-mail Alert Rss
 

物探与化探, 2019, 43(2): 257-265 doi: 10.11720/wtyht.2019.1196

地质调查·资源勘查

一种基于随钻VSP地震地质导向的钻井轨迹高效优化方法及其应用

杨敬雅1, 李相文1, 刘永雷1, 徐博1, 高江涛1, 王茂1, 吴江勇2, 张泉1

1. 中国石油集团 东方地球物理勘探有限责任公司,河北 涿州 072750

2. 中国石油集团 塔里木油田分公司勘探开发研究院,新疆 库尔勒 841000

An efficient method for drilling trajectory optimization based on seismological and geological guidance of VSP data during drilling and its application

YANG Jing-Ya1, LI Xiang-Wen1, LIU Yong-Lei1, XU Bo1, GAO Jiang-Tao1, WANG Mao1, WU Jiang-Yong2, ZHANG Quan1

1. Bureau of Geophysical Prospecting ( BGP) Company of CNPC,Zhuozhou 072750,China

2. Exploration and Development Research Institute of Tarim Oilfield Company,Korla 841000,China

责任编辑: 叶佩

收稿日期: 2018-05-22   修回日期: 2019-01-21   网络出版日期: 2019-04-20

基金资助: 国家科技重大专项“大型油气田及煤层气开发,塔里木盆地碳酸盐岩油气田提高采收率关键技术示范工程”.  2016ZX05053-17

Received: 2018-05-22   Revised: 2019-01-21   Online: 2019-04-20

作者简介 About authors

杨敬雅(1986-),女,西南石油大学本科毕业,工程师,主要从事高效解释方法研究和基地管理工作 。

摘要

塔里木盆地发现多个奥陶系碳酸盐岩缝洞型油气田,随着开发的不断推进,中小型储集规模的缝洞体已经成为产能建设的主力目标。勘探开发过程中,由于受到复杂地层速度的影响,缝洞体系偏移成像、归位存在一定误差,常导致钻井的失利。笔者提出一种基于随钻VSP地震地质导向的钻井轨迹高效优化方法,首先应用VSP速度与基于三角网格的新构造模型高效率地校正原始地震偏移速度场,再联合修正后的各向异性参数场对地震数据进行快速偏移成像。其次,对新数据的全层系地层的构造特征再分析,预测优质储层的发育位置,确定最终靶点。最后,根据工程、地质情况,采用最为经济的优化设计钻井轨迹,钻进入靶。应用表明,该技术在大沙漠腹地塔中地区应用取得了较好的效果,地震资料成像品质得到提升。调整轨迹后的储层钻遇率和钻完井直投率分别提升16%和13.2%,有效地促进了油田实施低成本效益开发方案。

关键词： 高效轨迹优化 ; VSP驱动校正速度场 ; 叠前偏移成像 ; 缝洞型碳酸盐岩 ; 塔里木盆地

Abstract

There are several Ordovician carbonate fractured-vuggy type oil and gas fields in the Tarim Basin.With the continuous development and production of the reservoirs,the small and medium sized fractured-vuggy systems have become the main targets of capacity development.In the process of exploration and development,due to the influence of complex formation velocity,there is a certain error in seismic migration data,which often leads to the failure of drilling.In this paper,an efficient drilling track optimization method based on VSP data for seismic data processing is proposed.First,the VSP velocity and the new triangular mesh model are used to effectively correct the original seismic migration velocity field,and then the modified anisotropic parameter field is used for quick migration.Secondly,the structural characteristics of the whole strata are analyzed in the new data,the positions of the reservoirs are predicted,and the targets are determined.Finally,according to the engineering and geological conditions,the most economical optimal design of the drilling trajectory is used to drill into the target.The application shows that the relatively good effect is obtained in Tazhong area of Tarim Oilfield,and the imaging quality of seismic data has been improved obviously.After well track adjustment,the rates of drilled reservoir and the directly put on production are increased by 16% and 13.2% respectively.The low cost and effective development plan of the oilfield is effectively promoted and implemented.

Keywords： efficient drilling track optimization ; VSP driven migration velocity feild ; pre-stack migration imaging ; fractured-vuggy carbonate ; Tarim Basin

PDF (7201KB) 元数据 多维度评价 相关文章 导出 EndNote| Ris| Bibtex  收藏本文

本文引用格式

杨敬雅, 李相文, 刘永雷, 徐博, 高江涛, 王茂, 吴江勇, 张泉. 一种基于随钻VSP地震地质导向的钻井轨迹高效优化方法及其应用. 物探与化探[J], 2019, 43(2): 257-265 doi:10.11720/wtyht.2019.1196

YANG Jing-Ya, LI Xiang-Wen, LIU Yong-Lei, XU Bo, GAO Jiang-Tao, WANG Mao, WU Jiang-Yong, ZHANG Quan. An efficient method for drilling trajectory optimization based on seismological and geological guidance of VSP data during drilling and its application. Geophysical and Geochemical Exploration[J], 2019, 43(2): 257-265 doi:10.11720/wtyht.2019.1196

0 引言

塔里木油田奥陶系碳酸盐岩油气藏属于缝洞型复合油气藏[1],其储集空间主要是以孔、缝、洞为基础且以不同组合形式赋存的非均质储层,缝洞系统结构复杂。近年来,在塔中地区碳酸盐岩油气藏开发中,以“高、大、强”为特点的“串珠”储层已基本被动用,剩下大量中小储集规模缝洞体未被动用,是油田产能建设的主力目标。针对此类较小规模储层的钻探,对储层的空间归位研究尤为重要,是保证钻井成功率的核心资料。

由于塔中I号气田区地处大沙漠区,起伏沙丘对地震波造成强烈的吸收衰减和散射等干扰,致使目的层奥陶系内幕信噪比较低。尽管地震成像技术的发展非常快,地震成像精度也得到了大幅度的提升,但受地震速度影响,地震成像不可能完全准确,仍然存在着一定的误差。尤其是针对目前大批量上钻的中小储集规模、非均质性极强的缝洞储层,这一误差往往会直接导致钻井的失利,需要通过定向侧钻揭开储层或大规模的储层改造,进而大大增加钻井成本。因此,必须在钻井的过程中对地震成像的误差进行修正,从而确保钻井能准确入靶。

目前较成熟的地震导向技术主要是以Schlumberger Seismic VisionTM为代表的随钻地震测量技术[2,3],虽然该技术具有实时测量,不损失钻井时间的优点,但仪器稳定性受井况影响较大,且使用费用高昂,不利于规模化应用。为此,我们应用基于随钻VSP地震数据方法确定靶点位置、优化入靶点轨迹,引导钻井准确入靶,达到地震地质导向的目的。该方法在塔中地区碳酸盐岩油气藏的井位钻探中成功应用与推广,进一步证实该方法的有效性和经济性。

1 随钻VSP地震导向技术

1.1 技术思路

为了获得高品质地震数据,引导钻井入靶,在钻井至目的层中之后实施VSP测井,并利用获取的时深关系数据对原始地震速度场进行修正,进而对地震数据进行重新偏移成像并提高地震资料的分辨率,最终利用新偏移地震数据确定靶点位置、优化入靶点轨迹,实施地震—地质—工程一体化作业,实现经济型钻井入靶点。方法核心就是通过钻井、测井获取准确地质分层实现构造模型的快速更新;通过VSP测量获取更精细的地层速度[4]、Q值[5-7]、甚至是各向异性参数场[8,9,10,11,12],实现各向异性深度偏移处理,获取高品质地震数据;通过精细储层预测确定靶点位置,最终实施轨迹优化调整(详见图1)。

1.2 缝洞储层成像速度敏感性分析

为了更清楚地认识速度精度对缝洞体成像效果的影响,研究中建立两个地质模型,通过正演模拟结果开展分析。设计水平地层条件下,目地层以上覆盖地层速度的总体变化量为2%,如图2所示,速度增加或减小2%即可导致“串珠”状储层的成像效果发生明显变化。其中,上覆地层速度减小2%后深度偏移结果中“串珠”状地震反射两侧出现下拉现象(图2b、图2e),上覆地层速度增加2%后深度偏移结果中“串珠”状地震反射两侧出现上翘现象(图2c、图2f)。同时,目的层及下伏内幕地层同相轴的深度也存在一定误差(图2a、图2b、图2c)。虽然总体速度误差值较小,但对串珠成像及空间归位均有非常大的影响,严重影响靶点的确定和钻井轨迹的优化设计。

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1   随钻VSP地震地质导向技术思路

a—随钻VSP地震地质导向技术总流程;b—TTI各向异性参数场建流程和地质解释;c—轨迹调整流程;流程中的速度为纵波垂直方向传播的速度;δ为变异系数(描述垂向速度的变化程度);ε为纵波各向异性强度(描述水平与垂直速度差异);θ为地层的倾角;Φ为方位角

Fig.1   The research route of VSP seismic geological guidance technology while drilling

a—the general process of VSP seismic geological guidance technology while drilling,including;b—TTI anisotropic parameter field construction process and geological interpretation;c—trajectory adjustment process;the velocity in the process is the velocity of vertical propagation of P-wave;δ:the coefficient of variation (describing the degree of variation of vertical velocity);ε:the intensity of P-wave anisotropy (describing the difference between horizontal and vertical velocity);θ:the dip angle of formation;Φ:the azimuth angle


图2

新窗口打开| 下载原图ZIP| 生成PPT

图2   水平地层时速度误差对“串珠”成像影响的正演分析

a—采用弹性波动方程模拟、逆时偏移处理方法的成果;b—采用与图a相同处理方法所得成果;c—采用与图a相同处理方法所得成果;d—图a对应的各向同性速度场;e—图b对应的各向同性速度场,黑色虚框内速度为图d速度场值的98%;f—图c对应的各向同性速度场,黑色虚框内速度为图d速度场值的102%

Fig.2   Forward analysis of the effect of velocity errors on "beads" imaging in horizontal formation

a—the results of elastic wave equation simulation and reverse time migration processing;b—the results obtained by the same processing method as figure a;c—the results obtained by the same processing method as figure a;d—the results obtained by the same processing method as figure a;e—the isotropic velocity field corresponds to figure b,the velocity in the black virtual frame is 98% of the value of the velocity field in figure d;f—the isotropic velocity field corresponds to figure c,the velocity in the black virtual frame is 102% of the value of the velocity field in figure d


另外,开展针对浅层横向速度变化对串珠偏移影响的正演模拟分析。4个单一参数变化的正演模型中,给定浅层地层速度变化率,即单位纵横向距离时横向速度变化量。如图3所示,图3a横向变化率为2%,图3b横向变化率为9%,图3c横向变化率为17%,图3d横向变化率为24%,其变化率导数为0位置处正下方倾斜地层中发育缝洞体,如图3中A、B、C、D处红方框所示意,其设计深度依次小20 m。图4为图3所示4个模型叠合深部地层缝洞体深度偏移剖面对比图。图4a中A缝洞体成像位置稍有偏离其真实位置,偏离值非常小,B、C和D缝洞体基本没有偏离其真实位置。图4b中A缝洞体成像位置相比图4a偏离值有增大,横向偏离值约约50 m,D缝洞体成像位置出现偏离现象,B、C没有偏离其真实位置。图4c中A缝洞体成像位置相比图4b偏离值又有增大,横向偏离值约100 m,D缝洞体成像位置偏离值的增幅远小于A缝洞体的偏离增幅。图4d中所示A缝洞体横向偏离值达到150 m,纵向偏离约一个相位,约50 m。随着横向速度变化的剧烈程度的增加,缝洞体的地震成像位置偏离其真实位置距离同样在增加,偏离方向为横向速度变小的方向,且缝洞体埋藏深度越大其偏离真实位置的距离值有增大的趋势。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3   针对浅部地层横向速度变化对串珠偏移归位影响分析的正演模型

Fig.3   A forward model for the effect of lateral velocity variation on migration and relocation of beads in shallow strata


图4

新窗口打开| 下载原图ZIP| 生成PPT

图4   图3所示正演模型叠合深度偏移成果对比

Fig.4   A comparison of the results of forward model stacking depth migration showed in the figure 3


正演分析认为,纵向位置不准的主要原因是纵向速度误差,横向偏移的主要因素是上覆地层速度横向误差。两个正演结果并结合研究区地震资料特点及目的层地质特征,影响缝洞储层成像的速度与VSP测量速度的误差最好控制在2%以内。

1.3 搭建基于三角网格算法的构造模型

由于实施VSP测井是在钻井至目的层中之后,因此基于原地震数据的构造预测值与实钻值误差已经确定,并将该误差快速且准确地修正到构造模型中以便偏移速度场的修正。研究中采用不规则角网(triangulated irregular network,TIN)高效地修正构造模型,以供后期速度场、各向异性参数场、Q值场的调整之需。

不规则角网最常用的生成方法是Delaunay创分方法[13,14,15]。TIN在表示复杂表面上具有许多优越性,被广就应用于数字制用、地用表面的模型化及分析以及LIS中。其中,在所有可能的三角网中,狄洛尼(Delaunay)三角网在地形拟合方面运用的较普遍,因此常被用于TIN的生成[16,17]。该算法对处理平面的几何特征有独特优势,网格顶点光滑分布,

均匀分布,几何相关的各向异性特征保持更好,网格的全局质量更易把控。

应用表明,深度偏移是解决深埋非均质缝洞型储层的空间归位的有效的处理方法之一[18,19,20]。因此,网格的均匀性、地层正确接触关系以及正确处理横向突变关系是叠前深度偏移处理速度建场中关键环节,也是断溶体油藏空间准确成像归位的关键。所以三角网格是一种更适于复杂地质构造建模的方法,建立研究区现今地层格架模型和断层模型,为后续的速度充填提供三维骨架,其质量直接决定着储层模型的质量。图4是基于三角网格构造建模效果与某商业软件常规构造建模效果对比图,图5a的网格均匀性明显优于图5b(蓝色箭头处),图5c对地层与断裂关系的处理效果明显优于图5d(红色虚框处)。另外,基于三角网格的精细构造建模可快速有效地处理断层与地层的交切关系,提高精细构造建模的效率,为后续VSP测量速度约束迭代建场与优化奠定坚实基础。

1.4 VSP驱动迭代速度分析

准确地建立地震数据偏移成像速度场有两个关键环节,即VSP垂向速度约束和地质导向横向速度约束。具体是先要应用VSP速度进行校正原始地震速度场[21,22,23,24,25],其次要根据井点周围各层系地层构造特征约束井周的速度趋势,最后修正各向异性参数场。应用修正后的速度场对地震数据进行快速偏移成像,校正井—震误差。利用基于三角网格精细构造建模与VSP数据驱动优化迭代地震偏移速度场效率明显优于常规手动调整速度效率的特点,使得在有限时间内,通过数次处理解释一体化的优化迭代过程将地震资料品质提升一个台阶,满足油田开发中完井时间段内提出靶点调整与优化轨迹的需求。

如图6a所示,蓝色为VSP测量速度,绿色为校正前地震速度,红色为校正后地震速度,原始地震速度整体变化规律与VSP速度基本一致,但是局部井段的速度细节有较大误差,校正后地震速度与VSP速度基本吻合。对比校正前后地震速度,如图6b和图6c,速度精细度发生明显的变化。同时,如图6d所示,合成记录与新偏移成像的地震数据对比标志层的反射特征一致,全井段波形相似系数为0.88,进一步证实了更新后的地震速度场的全层系地层速度更加准确。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5   基于三角网格算法构造建模效果与Petrel常规构造建模效果对比

a—基于三角网格算法构造建模结果中某一地层平面图格;b—Petrel常规构造建模结果中与图a所示相同位置地层平面图格;c—基于三角网格算法构造建模结果A-A’剖面;d—Petrel常规构造建模结果A-A’剖面

Fig.5   Comparison of modeling effect between triangular mesh algorithm and Petrel conventional construction

a—based triangular mesh algorithm is used to construct a stratigraphic plane grid in the modeling results;b—the conventional structural modeling results using Petrel in the same as the stratigraphic plane grid shown in figure a;c—the constructing A-A’ section based on triangular mesh algorithm;d—conventional structure modeling A-A’ section using Petrel


图6

新窗口打开| 下载原图ZIP| 生成PPT

图6   基于VSP驱动处理前后速度对比图(曲线、剖面)与井震标定剖面

a—VSP速度、校正前地震速度和校正后地震速度叠合对比;b—原始地震速度场连井剖面;c—应用VSP数据修正后的速度场连井剖面;d—基于新偏移成像的地震数据的VSP标定结果

Fig.6   Velocity contrast map (curve,profile) and well-seismic calibration profile based on VSP driving processing

a—the stacking contrast of VSP velocity,pre-correction seismic velocity and post-correction seismic velocity;b—original seismic velocity field cross-well profile;c—the velocity field cross-well profile modified by VSP data;d—the VSP calibration results based on new migration imaging seismic data


1.5 地震—地质—工程一体化钻井轨迹优化调整

应用VSP速度信息更新地震速度场后,对CMP道集进行偏移处理,解释员对新数据进行品质评价,分析全层系地层的构造误差及预测储层,实现对储层、油藏的再认识,最终确定靶点,实现地震—地质—工程一体化钻井轨迹优化调整。根据目前采用的技术流程,到最终提出井轨迹调整方案的时间基本控制在5天以内,与钻井工程安装技术套管同时进行,基本不影响工程进度,具有较好的时效性。

以塔中地区WELL1井为例,基于VSP数据修正地震速度后的偏移处理数据与老资料对比,目的层段顶面构造总体形态基本一致,局部细节稍有变化。另外,新资料分辩率更高,“串珠”反射形态没有发生较大的变化,只是两套地震资料所揭示“串珠”状储层的空间位置发生了小幅度的偏移,从图7a~图7d可以看到,“串珠”状储层在纵向上偏浅了21 m,在平面上“串珠”中心向南偏移了30 m。通过临井对比分析,认为基于VSP数据修正地震速度后的偏移数据的储层成像位置更准。因此,以此为基础并结合钻井工程狗腿度、方位纠正等实施的可行性分析结论对原始钻井轨迹进行了调整,靶点向南调23 m,深度向上调整24 m。实钻表明,按照调整后轨迹钻进图7a零相位处(海拔 -4 790.72 m)发生放空、漏失,累计放空0.6 m,累计漏失钻井液1 123m3,气测解释TG值达到 99.63%,后5 mm油嘴,试油,油压45.618 MPa获得高产工业油流。WELL1的投产进一步说明,随钻地震地质导向技术是准确可靠的。

2 应用效果

通过文中所提出的方法在塔中地区的应用,提高了储层认识精度,对原有地震数据体进行校正,使“串珠”精确归位,不仅在钻井随钻跟踪中应用效果显著,而且校正过的地震数据体对新井部署(图8)及老井措施建议等应用方面也具有很显著的作用。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7   WELL1井基于VSP驱动处理前后地震偏移剖面、均方根振幅属性对比

a—基于VSP数据驱动地震速度后的偏移处理数据地震剖面;b—老数据地震剖面;c—基于VSP数据修正地震速度后的偏移处理数据的储层预测平面;d—基于老数据的储层预测平面

Fig.7   Seismic migration profile and RMS amplitude attribute contrast map of WELL1 well before and after VSP drive processing

a—seismic profile of migration processing data driven by VSP data;b—old seismic data profile;c—reservoir prediction plane pap of migration processing data modified by VSP data;d—reservoir prediction plane map based on old data


图8

新窗口打开| 下载原图ZIP| 生成PPT

图8   基于多井VSP联合驱动处理前后数据的均方根振幅属性对比

Fig.8   Contrast map of RMS amplitude attributes of data before and after multi-well VSP combined drive processing


从2016年起截至目前,塔中地区完钻井共79口。其中应用VSP数据驱动处理数据的钻井轨迹高效优化方法累计随钻跟踪钻井26口,完钻25口,除W-296井外,均发生不同程度的放空和漏失,储层钻遇率达到96%(表1)。后期开发中直接测试投产井16口,储层改造井6口,完井直投率为72.2%。

表1   塔中地区应用VSP驱动处理数据的随钻跟踪钻井漏失放空及储层改造情况统计

Table 1  Statistical table of lost circulation and emptying and reservoir modification during drilling with VSP-driven data processing in Tazhong area

序号井号年度完钻层位漏失放空酸压/m3
层段/m漏失量/m3井段/m长度/m
1W-511-7X2016O3l5405.94~5415506
2W-2652016O2y20590
3W-43-3X2017O1-2y25452.6~55202994.3
4W-2932017O2y6184.4~62491306
5W-112-2X2017O1-2y16472.6~6545668.86534.95~6535.57
6528.68~6529.35		1.29
6W-102-3X2017O3l6718.93~67516256748.65~6749.9
6750.08~6751		2.17633
7W-7-9X2017O3l5628.31~58131265.8
8W-83-7X2017O3l5786~579173.1
9W-291-5X2017O2y6270.12~6300.121035.56272.28~6273.56
6274.84~6275.38		1.82
10W-1-1X2017O3l5902.75~596011235906.2~5906.80.6
11W-8012017O3l5850.11~5862607.5470
12W-8-2X2017O1-2y16155.63~6195797.09240
13W-2962017O2y740
17W-2952017O2y6078.32~6095.520.56192.47~6192.9
6192.9~6193.01		0.54
18W-102-6X2017O3l6549.98~657079.2740
19W-1132017O1-2y26861.85~68801371
20W-726-1X2018O3l5557.96~55581546.9
21W-83-8X2018O1-2y15622~5635766.9
22W-6-62018O3l5782.3~58077.8
23W-15-7X2018O2y6091.38~6124.361246.2
24W-43-72018O1-2y25948.15~5976456.7
25W-113-22018O1-2y27204.56~7210246.3
26W-45-62018设计O1-2y2—(正钻)—(正钻)

新窗口打开| 下载CSV


常规随钻完井57口,其中11口钻井未发生放空和漏失,储层钻遇率80%。后期开发直接测试投产井34口,储层改造井23口,完井直投率为59%。对比可以看出VSP驱动处理数据可有效地提高储层钻遇率。

3 结论

通过本方法的研究及实际应用效果,有以下几点结论:

1)正演分析认为,储层位置空间偏移归位不准的主要原因是地震偏移速度与真实地层速度存在一定误差,针对本次研究区内影响缝洞储层成像的速度与VSP测量速度的误差应控制在2%以内;

2)基于三角网格的精细构造建模,能正确处理构造、地层横向突变问题,具有精度高、高效率的特点,为迭代更新地震偏移速度场、Q值场奠定基础;

3)VSP驱动处理可有效提升地震资料的成像效果,可实现对储层、油藏特征的再认识,为靶点的确定和井轨迹的优化提供资料依据;

4)地震—地质—工程一体化钻井轨迹调整,基本没影响钻井工程进度,且具有明显的经济优势。

5)该方法同样可以应用于山前复杂构造、超薄储层等的钻探中,在低油价形势下有非常广阔的推广应用前景。

致谢:本文的顺利完成,得益于东方地球物理公司各位领导的指导,公司研究院库尔勒分院的同事帮助,以及塔里木油田分公司的支持,在此一并表示感谢!

The authors have declared that no competing interests exist.
作者已声明无竞争性利益关系。

参考文献

杜金虎, 周新源, 李启明 , 等.

塔里木盆地碳酸盐岩大油气区特征与主控因素

[J]. 石油勘探与开发, 2011,38(6):652-661.

DOI:10.1016/S1876-3804(12)60002-0      URL     Magsci     [本文引用: 1]

A comprehensive analysis of hydrocarbon accumulation conditions and typical reservoirs reveals that a large carbonate petroleum province is formed around the Manjiaer depression in the Tarim Basin and it is consisted of Ordovician non-structural petroleum reservoirs. Multi-type litho-stratigraphic reservoirs are stacked vertically and connected laterally along the paleo-uplift slopes, with a large area. Reservoirs are intensively heterogeneous, extensively developed and distributed in layers along the top of different-scale sequence boundaries. Hydrocarbon phase and physical properties are diverse and the distribution pattern is “lower gas, upper oil”. The oil/gas/water output alters and declines quickly, and 20%1730% high-efficiency wells provide more than 70% of the production. The stable paleo-uplift, widely developed heterogeneous reservoir, multi-stages of hydrocarbon accumulation are the basic forming and distribution conditions of the large carbonate petroleum province. The paleo-uplift near the hydrocarbon kitchen is developed stably and controls the reservoir distribution, migration and accumulation. The special characteristics of carbonate reservoirs result in the wide and discontinuous distribution of irregular fracture-cavity hydrocarbon pools. The enrichment of hydrocarbon is controlled by reservoirs, which are the main controlling factor of the formation and distribution of the large petroleum province. Multi-stages of accumulation and adjustment of hydrocarbon are universal and result in the large area and diversity of the oil and gas.

Du J H, Zhou X Y, Li Q M , et al.

Characteristics and controlling factors of the large carbonate petroleum province in the Tarim Basin,NW China

[J]. Petroleum Geophysical Exploration, 2011,38(6):652-661.

Magsci     [本文引用: 1]

王绍先 .

随钻地震技术的研究及应用

[D]. 北京:中国科学院地质与地球物理研究所, 2008.

[本文引用: 1]

Wang S X .

Research and application of seismic technology while drilling

[D].Beijing:Institute of Geology and Geophysics, Chinese Academy of Sciences, 2008.

[本文引用: 1]

Celniker G W . Performing geoscience interpretation with simulated data[P].US, US6256603B1, 2001.

[本文引用: 1]

王成礼, 叶古清 .

利用非零井源距VSP资料计算层速度

[J]. 石油地球物理勘探, 1988,23(5):616-619.

URL     Magsci     [本文引用: 1]

问题的提出随着垂直地震剖面法的广泛应用,用它计算层速度的方法也日益增多。由于对射线路径的假设条件不同,因而导出了各种计算层速度的方法。其中以直达波直射线传播路径为假设条件的计算层速度方法用得较多。

Wang C L, Ye G Q .

Calculating formation velocity using VSP data of non-zero well spacing

[J]. Petroleum Geophysical Exploration, 1988,23(5):616-619.

Magsci     [本文引用: 1]

朱光明, 高静怀 .

用VSP资料估算Q值

[C]//中国地球物理学会学术年会, 1991.

[本文引用: 1]

Zhu G M, Gao J H .

Estimating Q value with VSP data

[C]//Annual Meeting of the Chinese Geophysical Society, 1991.

[本文引用: 1]

张大伟, 孙赞东, 王学军 , 等.

利用零井源距VSP资料进行品质因子反演

[J]. 石油地球物理勘探, 2011,46(s1):47-52.

URL    

品质因子Q值是描述地震波衰减的一个重要物理参数,它描述了介质非完全弹性特征,是地层内部本质特征,所以研究和估算Q值在地震勘探中的作用越来越重要。谱比法和质心频移法是现在常用的两种方法。本文详细介绍了两种方法的原理并用无噪声理论模型验证了两种方法的可行性。反演结果显示在无噪声的条件下两种方法都能很好地估算介质的品质因子。通过对含有不同强度噪声的理论数据、全波场理论数据和实际零井源距VSP数据的处理,对两种方法进行了比较。结果表明,质心频移法反演的结果比谱比法更精确、更可信。

Zhang D W, Sun Z D, Wang X J , et al.

Quality factor inversion using zero well spacing VSP data

[J]. Petroleum Geophysical Exploration, 2011,46(s1):47-52.

王正和, 崔永福, 向东 .

井控处理中的真振幅恢复与Q补偿方法及应用

[J]. 物探与化探, 2008,32(4):434-437.

URL     Magsci     [本文引用: 1]

<p>在地震资料处理中,可以使用井控地震处理技术,即通过井资料求取球面扩散补偿因子和<em>Q补偿因子。在获得的Q值附近进行扫描,将不同Q补偿的偏移剖面与井资料进行匹配分析,最终获得精确的Q值。实际资料真振幅恢复和Q补偿结果表明,该方法有效地实现了地震资料真振幅恢复和Q</em>补偿,并取得了满意的结果。</p>

Wang Z H, Cui Y F, Xiang D .

The method and application of true amplitute recovery and Q-compensation in well-control processing

[J]. Geophysical and Geochemical Exploration, 2008,32(4):434-437.

Magsci     [本文引用: 1]

Bakulin A, Slater C, Bunain H, 等.

在有横向非均匀性情况下由多方位变井源距VSP估算方位各向异性和裂缝参数

[C]//SEG年会论文概要, 2000.

[本文引用: 1]

Bakulin A, Slater C, Bunain H , et al.

In the case of transverse heterogeneity,estimating azimuthal anisotropy and fracture parameters by multi-azimuthal variable well spacing VSP

[C]//SEG Annual Conference Paper Summary, 2000.

[本文引用: 1]

侯爱源, 耿伟峰, 张文波 .

Walkaway VSP井旁各向异性速度模型反演方法

[J]. 石油地球物理勘探, 2016,51(2):301-305.

DOI:10.13810/j.cnki.issn.1000-7210.2016.02.013      URL     Magsci     [本文引用: 1]

针对一个井源距VSP初至反演的井旁速度模型不能用于同一口井多个不同井源距VSP数据的偏移成像这一问题,本文介绍了一种使用Walkaway VSP的初至旅行时反演井旁各向异性模型的层速度和各向异性参数的方法。该方法根据声波测井速度的变化趋势对井旁一维速度模型进行分层,以零井源距VSP初至时间反演的速度作为初始速度模型,用一种接近沉积岩特性的多层各向异性介质模型描述和定义井旁速度模型的各向异性参数,利用不同井源距VSP的多条VSP初至时距曲线反演井旁水平层状各向异性速度模型的层速度和各向异性参数,得到Walkaway VSP偏移成像的井旁各向异性速度模型。把该模型用于Walkaway VSP实际数据的偏移成像处理,取得了较好的效果。

Hou A Y, Geng W F, Zhang W B .

Inversion method of anisotropic velocity model near wellbore by Walkaway VSP

[J]. Petroleum Geophysical Exploration, 2016,51(2):301-305.

Magsci     [本文引用: 1]

黄翼坚 .

多井源距VSP速度分析及逆时偏移

[D]. 西安:长安大学, 2010.

[本文引用: 1]

Huang Y J .

Velocity analysis and reverse-time migration of multi-well source spacing VSP

[D]. Xi’an:Chang’an University, 2010.

[本文引用: 1]

罗小明, 凌云, 牛滨华 , 等.

井地联合提取VTI介质各向异性参数

[J]. 地学前缘, 2005,12(4):576-580.

DOI:10.3321/j.issn:1005-2321.2005.04.028      URL     [本文引用: 1]

在新疆塔里木盆地某地区的宽方位三维地震数据中,VTI介质火成岩的各向异性和垂向速度变化引起非双曲线时差.通过对实际资料的系统分析,我们不能利用DMO反演法,移动震源VSP反演法和多偏移距VSP初至旅行时反演法提取Thomsen参数用于各向异性地震资料的处理.因此本文提出一个新的方法:运用地震资料和零偏移距VSP井求取火成岩的各向异性参数δ;在比较精确确定火成岩地层深度和厚度的情况下,利用偏移距VSP井,求取了各向异性参数ε,从而获得火成岩的瞬时各向异性参数η;充分考虑VTI介质垂向速度变化对时差曲线的影响,获得了各向异性参数ηeff.最后,利用所提取的各向异性参数,进行VTI介质各向异性速度分析和成像,成像质量得到明显改善.

Luo X M, Ling Y, Niu B H , et al.

Joint extraction of anisotropic parameters of VTI media

[J]. Geoscience Frontier, 2005,12(4):576-580.

[本文引用: 1]

Liu Y M, Liang X H, Yin X Y , et al.

Estimation of Thomsen’s anisotropic parameters from walkaway VSP and applications

[J]. Applied Geophysics, 2014,11(1):23-30.

DOI:10.1007/s11770-013-0410-z      URL     [本文引用: 1]

Estimation of Thomsen’s anisotropic parameters is very important for accurate time-to-depth conversion and depth migration data processing. Compared with other methods, it is much easier and more reliable to estimate anisotropic parameters that are required for surface seismic depth imaging from vertical seismic profile (VSP) data, because the first arrivals of VSP data can be picked with much higher accuracy. In this study, we developed a method for estimating Thomsen’s P-wave anisotropic parameters in VTI media using the first arrivals from walkaway VSP data. Model first-arrival travel times are calculated on the basis of the near-offset normal moveout correction velocity in VTI media and ray tracing using Thomsen’s P-wave velocity approximation. Then, the anisotropic parameters and are determined by minimizing the difference between the calculated and observed travel times for the near and far offsets. Numerical forward modeling, using the proposed method indicates that errors between the estimated and measured anisotropic parameters are small. Using field data from an eight-azimuth walkaway VSP in Tarim Basin, we estimated the parameters and and built an anisotropic depth-velocity model for prestack depth migration processing of surface 3D seismic data. The results show improvement in imaging the carbonate reservoirs and minimizing the depth errors of the geological targets.

谭仁春, 杜清运, 杨品福 , 等.

地形建模中不规则三角网构建的优化算法研究

[J]. 武汉大学学报信息科学版, 2006,31(5):436-439.

DOI:10.3321/j.issn:1671-8860.2006.05.015      URL     [本文引用: 1]

对利用数字线划图(DLG)重建地形三维模型进行了研究,提出了基于三角形拓扑关系的自动联结三角网算法,并构造了复杂地形的三维数据模型。基于此模型,设计实现了构网优化算法,实现了复杂地形的自动三维重建。试验证明,本文提出的构网优化算法是稳健、高效的。

Tan R C, Du Q Y, Yang P F , et al.

Study on the optimization algorithm of irregular triangular network construction in terrain modeling

[J]. Journal of Wuhan University Information Science Edition, 2006,31(5):436-439.

[本文引用: 1]

花卫华, 邓伟萍, 刘修国 , 等.

一种改进的不规则三角网格曲面切割算法

[J]. 地球科学, 2006,31(5):619-623.

DOI:10.3321/j.issn:1000-2383.2006.05.008      URL     [本文引用: 1]

三角网切割是实现三维地质建模和模型分析的关键算法，它的效率直接决定了建模算法的效率．通过建立三角网的方向包围盒（oriented bounding box，OBB）树实现曲面间的碰撞检测，然后对发生相交的三角形对统一计算交点，通过对顶点坐标归一化完成曲面切割后的快速重构，并针对不同的曲面类型采用不同的切割分类方法进行2侧切割结果的划分．阐述了算法的实现过程并展示了切割后的图形效果．

Hua W H, Deng W P, Liu X G , et al.

An improved algorithm for irregular triangular mesh surface cutting

[J]. Geoscience, 2006,31(5):619-623.

[本文引用: 1]

朱庆, 陈楚江 .

不规则三角网的快速建立及其动态更新

[J]. 武汉测绘科技大学学报, 1998,23(3):204-207.

DOI:10.2190/K5DG-LENJ-AMDR-591P      URL     [本文引用: 1]

针对线路勘测设计的特点，介绍了一种适合于线路ＣＡＤ应用的快速建立和在线动态更新不规则三角网（ＴＩＮ）的方法。该方法根据带状数据分布快速生成“类Ｄｅｌａｕｎａｙ三角形”，并实现了ＴＩＮ的局部更新，从而为线路ＣＡＤ提供了基本的数学基础。最后用实例验证了该算法的正确性。

Zhu Q, Chen C J .

Rapid establishment and dynamic updating of irregular triangular networks

[J]. Journal of Wuhan University of Surveying and Mapping Science and Technology, 1998,23(3):204-207.

[本文引用: 1]

吴宇晓, 张登荣 .

生成三角网的快速合成算法

[J]. 浙江大学学报理学版, 2004,31(3):343-348.

DOI:10.3321/j.issn:1008-9497.2004.03.023      URL     [本文引用: 1]

合成算法结合了传统的递归分割法和逐点插入法的优点,兼顾空间和时间性能.然而,该算法不可避免地继承了两种传统算法的不足,在执行效率上受到限制.为了解决执行效率问题,提出了快速合成算法,对合成算法进行了改进和优化.该算法基于面积坐标的点定位算法和简化的高效空外接圆判断算法,从而大大提高算法的整体执行效率;同时充分考虑平面点集的任意性,适用于对任意平面点集构建Delaunay三角网.

Wu Y X, Zhang D R .

Fast synjournal algorithm for generating triangular networks

[J]. Journal of Zhejiang University, Science Edition, 2004,31(3):343-348.

[本文引用: 1]

郭良辉, 孟小红, 郭志宏 , 等.

地球物理不规则分布数据的空间网格化法

[J]. 物探与化探, 2005,29(5):438-442.

DOI:10.3969/j.issn.1000-8918.2005.05.016      URL     Magsci     [本文引用: 1]

<p>不规则分布数据的网格化处理是地球物理数据处理和解释的基础问题, 是保证许多地球物理数据处理方法得以成功实施的前提。文中介绍了线性插值法、多元二次函数法、普通克里格法、反插值法等4种空间网格化法的原理, 并利用这些网格化方法进行了理论模型和实际航磁数据的试验分析。</p>

Guo L H, Meng X H, Guo Z H , et al.

Gridding methods of geophysical irregular data in space domain

[J]. Geophysical and Geochemical Exploration, 2005,29(5):438-442.

Magsci     [本文引用: 1]

彭更新, 但光箭, 郑多明 , 等.

塔里木盆地哈拉哈塘地区三维叠前深度偏移与储层定量雕刻

[J]. 中国石油勘探, 2011,16(5):49-56.

DOI:10.3969/j.issn.1672-7703.2011.05-06.008      URL     [本文引用: 1]

塔里木盆地哈拉哈塘奥陶系碳酸盐岩缝洞型储层十分发育,裂缝溶洞 的展布主要受断裂、古水系控制,油层厚度小,水体活跃.缝洞储层具有大小规模不同、空间分布形态不规则、充填性质各异的缝洞组合特征.如何准确刻画碳酸盐 岩缝洞储层空间位置、缝洞储层容积、油气水分布特征,已成为碳酸盐岩油藏勘探开发井位部署的主要难点之一.在高品质的高密度全方位采集的地震资料的保障 下,提出了叠前深度偏移、正反演结合定量雕刻、叠前地震各向异性分析和叠前弹性参数反演预测方法.实际地震资料的应用表明,叠前深度偏移可以准确刻画缝洞 体空间位置,叠前各向异性分析能有效预测裂缝发育强度及方位,结合井资料叠前弹性反演能有效预测缝洞储层空间形态和储层的含油气性.多种手段的应用,有效 地解决了碳酸盐岩缝洞油藏研究的难题,为哈拉哈塘碳酸盐岩油藏的勘探开发井位部署提供了有力支持.

Peng G X, Dan G J, Zheng D M , et al.

3D Prestack depth migration and quantitative characterization of carbonate reservoirs in Halahatang Area,Tarim Basin

[J]. China Petroleum Exploration, 2011,16(5):49-56.

[本文引用: 1]

黄元溢, 罗仁泽, 王进海 , 等.

几种叠前深度偏移技术效果的对比

[J]. 物探与化探, 2011,35(6):798-803.

URL     Magsci     [本文引用: 1]

<p>在回顾叠前深度偏移发展历程的基础上,概述了当前几类叠前深度偏移方法:克希霍夫积分法、波动方程法、逆时偏移。在分析叠前深度偏移的波场成像方法、适用性、近地表条件敏感性及存在问题的基础上,对比偏移剖面,比较不同偏移方法的成像效果,指出叠前深度偏移的发展方向,为以后的叠前深度偏移的研究提供参考建议。</p>

Huang Y Y, Luo R Z, Wang J H , et al.

A comparison of the application effects of several prestack depth migration techniques

[J]. Geophysical and Geochemical Exploration, 2011,35(6):798-803.

Magsci     [本文引用: 1]

常建华, 张玉娥, 林卫 , 等.

三维叠前深度偏移处理技术应用

[J]. 物探与化探, 1999,23(4):271-276.

DOI:10.3969/j.issn.1000-8918.1999.04.005      URL     [本文引用: 1]

近年来，三维叠前深度偏移地震资料处理技术得到了一定的发展，但是由于三维叠前深度偏移处理计算量大，要求有精确的深度－层速度模型，一直未同工业化生产阶段。本文是在具备大型计算机的基础上，着重对如何求取精确的深度－层速度模型进行了探讨，提出了三维层速度相干反演和三维深度－层速度模型的优化处理方法。

Chang J H, Zhang Y E, Lin W , et al.

Application of 3-D prestack depth migration processing technology

[J]. Geophysical and Geochemical Exploration, 1999,23(4):271-276.

[本文引用: 1]

胡英, 姚逢昌 .

用于叠前深度偏移速度建模的一种新方法

[J]. 石油勘探与开发, 2000,27(2):62-64.

DOI:10.3321/j.issn:1000-0747.2000.02.017      URL     [本文引用: 1]

叠前深度偏移是解决复杂地质构造成像问题的重要技术。高质量的叠前深度偏移成像依赖于高精度的深度 速度模型。常规叠前深度偏移速度建模方法不仅要求提供高质量的叠前地震数据 ,而且迭代建模过程较长 ,从而限制了这项技术的推广应用。针对我国东部老油区钻井较多的特点 ,提出测井约束建模法 ,此方法利用已知声波测井、VSP等资料中的低频速度信息 ,帮助处理人员建立初始的宏观层速度模型 ;再根据叠前深度偏移之后 ,深度剖面上的主要层位应与测井中相应层位的深度值一致的约束条件来修正层速度模型。通过理论模型试算和实际应用证明 ,该方法实用有效 ,能够在缩短叠前深度偏移周期、节省机时费用的基础上提高叠前深度偏移处理

Hu Y, Yao F C .

A new method for prestack depth migration velocity modeling

[J]. Petroleum Exploration and Development, 2000,27(2):62-64.

[本文引用: 1]

陈生昌 .

波动方程叠前深度偏移成像与速度建模方法研究

[D]. 上海:同济大学, 2002.

[本文引用: 1]

Chen S C .

Wave equation prestack depth migration imaging and velocity modeling method

[D]. Shanghai:Tongji University, 2002.

[本文引用: 1]

马彦彦, 李国发, 张星宇 , 等.

叠前深度偏移速度建模方法分析

[J]. 石油地球物理勘探, 2014,49(4):687-693.

URL     Magsci     [本文引用: 1]

基于层位的层析反演和基于网格的层析反演是实际地震资料叠前深度偏移的两种主要速度建模方法。根据歧口凹陷典型地质模型及其正演模拟的地震波场，对两种速度建模方法进行了实验分析。实验结果表明：基于层位的层析反演方法具有较高的稳定性，但只能得到速度模型的低频分量;基于网格的层析反演方法可以获得速度场的高频分量，但受初始模型的影响较大，不容易收敛到实际速度模型。为此联合应用两种方法，即先利用基于层位的层析反演方法获得速度场的低频分量，再利用基于网格的层析反演方法获得速度场的高频分量，这样能够提高速度模型精度，改善叠前深度偏移的成像质量。

Ma Y Y, Li G F, Zhang X Y , et al.

Prestack depth migration velocity modeling method analysis

[J]. Petroleum Geophysical Exploration, 2014,49(4):687-693.

Magsci     [本文引用: 1]

高书琴, 冯许魁, 郭向宇 , 等.

波动方程叠前深度偏移处理技术开发与应用

[J]. 天然气工业, 2007,27(s1):205-207.

URL     [本文引用: 1]

针对地质结构复杂、岩性变化大、勘探目的层深、速度横向变化大地区的地震数据准确偏移成像,已成为目前石油工业界所面临的一大难题;传统的Kirchhoff积分法偏移固有其很强的适应性和高效的运算能力,但其射线追踪计算旅行时和振幅保持等方面的局限性越来越难以满足这些构造偏移成像的要求;为满足复杂油气田勘探开发的需要,对具有更高成像精度的波动方程叠前偏移成像技术的开发及应用必将成为地震偏移成像技术发展的必然趋势。

Gao S Q, Feng X K, Guo X Y , et al.

Development and application of wave equation prestack depth migration processing technology

[J]. Natural Gas Industry, 2007,27(s1):205-207.

[本文引用: 1]

谢会文, 许永忠, 郑多明 , 等.

火成岩速度反演及三维速度场建立的关键技术及效果分析

[J]. 物探与化探, 2013,37(6):1071-1079.

DOI:10.11720/j.issn.1000-8918.2013.6.22      URL     Magsci     [本文引用: 1]

<p>以地震反演为主，结合测井、岩石物理、地震相和正演模拟等技术综合研究火成岩的速度场特征。岩石物理实验是研究火成岩物性参数有效的工具，实验数据的分析能帮助预测火成岩岩性和速度信息，为反演等工作利用地震资料描述火成岩提供依据。原始测井曲线的环境校正和标准化处理是高品质合成记录和反演效果的前提条件。稀疏脉冲反演、基于模型的反演和神经网络反演这三种方法结果的分辨率较原始地震数据有了很大提高，从中可识别出火成岩的空间展布，但基于模型的反演分辨率相对较高。伽马属性反演的分辨率在识别火成岩方面比上述三种方法要高得多，但是单一的属性又很难满足火成岩特征描述的要求，因此多参数融合是一种行之有效的方法。结合地震相的应用，最终得到符合地质特征的火成岩三维速度场，并通过误差分析和正演模拟验证了速度场精度，为变速度成图和低幅构造勘探提供了科学的依据。</p>

Xie H W, Xu Y Z, Zheng D M , et al.

An analysis of the key techniques and effects of the igneous rock velocity inversion and the establishment of 3D velocity field

[J]. Geophysical and Geochemical Exploration, 2013,37(6):1071-1079.

Magsci     [本文引用: 1]

/

〈 〉

京ICP备05055290号-3
版权所有 © 2021《物探与化探》编辑部
通讯地址：北京市学院路29号航遥中心 邮编：100083
电话：010-62060192；62060193 E-mail：whtbjb@sina.com