CN113568041B - 时移地震三维拖缆采集数据的可重复性分析方法及系统 - Google Patents
时移地震三维拖缆采集数据的可重复性分析方法及系统 Download PDFInfo
- Publication number
- CN113568041B CN113568041B CN202110820603.6A CN202110820603A CN113568041B CN 113568041 B CN113568041 B CN 113568041B CN 202110820603 A CN202110820603 A CN 202110820603A CN 113568041 B CN113568041 B CN 113568041B
- Authority
- CN
- China
- Prior art keywords
- data
- navigation file
- dimensional
- seismic
- base data
- 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.)
- Active
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 34
- 238000012544 monitoring process Methods 0.000 claims abstract description 107
- 210000003746 feather Anatomy 0.000 claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 33
- 238000013508 migration Methods 0.000 claims abstract description 28
- 230000005012 migration Effects 0.000 claims abstract description 28
- 238000012545 processing Methods 0.000 claims abstract description 28
- 238000003384 imaging method Methods 0.000 claims abstract description 16
- 238000004859 neutralization-reionization mass spectrometry Methods 0.000 claims description 17
- 238000004422 calculation algorithm Methods 0.000 claims description 15
- 238000003860 storage Methods 0.000 claims description 10
- 238000004590 computer program Methods 0.000 claims description 9
- 238000006243 chemical reaction Methods 0.000 claims description 3
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 abstract description 6
- 239000003345 natural gas Substances 0.000 abstract description 3
- 239000003209 petroleum derivative Substances 0.000 abstract description 3
- 238000009826 distribution Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 9
- 238000010276 construction Methods 0.000 description 5
- 230000000704 physical effect Effects 0.000 description 5
- 230000008859 change Effects 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 230000008569 process Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 2
- 238000005520 cutting process Methods 0.000 description 2
- 238000009792 diffusion process Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000009469 supplementation Effects 0.000 description 2
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 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/282—Application of seismic models, synthetic seismograms
-
- 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/301—Analysis for determining seismic cross-sections or geostructures
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
技术领域
本发明是关于一种时移地震三维拖缆采集数据的可重复性分析方法及系统，属于石油天然气地震勘探领域。
背景技术
时移地震三维拖缆采集是目前海洋油气勘探最主要的地震数据采集方式，它具有工作效率高、成本低的优点。拖缆采集时，物探船拖拽着数条等间隔排放的长达数千米的接收电缆在海面上匀速航行，电缆排列和拖船之间放置着气枪阵列，通过瞬间释放高压空气产生地震波，地震波向下传播经地层反射后被接收电缆上的水听器接收。受成本考虑约束，到目前为止，三维拖缆采集仍然是海上时移地震勘探的首选采集方式。
时移地震勘探利用油田开发前后两次采集地震数据之间的差异，揭示储层物性变化，预测剩余油的分布。时移地震要求两次地震采集处理保持很好的可重复性，完成匹配处理后，非目的层的两次数据的差异一般接近于零，而差异主要集中在产油层。然而物探船作业时，受风浪、洋流的影响，漂浮在海面以下一定沉放深度的接收电缆常常会随着偏移距的增大逐渐偏移航线，与航线产生一个羽角。有时情况变得更为复杂，匀速航行的拖船不能将水中的采集电缆拉直，电缆在海平面上呈现曲线展布，大大降低时移地震基数据(油田开发前采集的地震数据或该地区前一次采集的地震数据)与监测数据(油田开发生产一段时间后再次采集的地震数据或该地区后一次采集的地震数据)的一致性。
目前，实际生产中对时移地震三维拖缆采集数据的可重复性分析均是基于实际数据的一致性分析，计算标志层位的NRMS值是否小于0.1的要求。这需要在基数据和监测数据都已实际采集后才能完成，它相当于已经完成了油田时移地震资料的处理，处理难度(因为影响因素很多，例如观测系统、震源容量、船速、潮汐)和工作量均很大。
发明内容
针对上述问题，本发明的目的是提供一种处理难度小且工作量小的时移地震三维拖缆采集数据的可重复性分析方法及系统。
为实现上述目的，本发明采取以下技术方案：一种时移地震三维拖缆采集数据的可重复性分析方法，包括：
根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件；
建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据；
对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体；
根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性。
进一步地，所述根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件，包括：
根据勘探地区的实际采集数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件；
将实际采集数据的导航文件作为基数据的导航文件，将生成的导航文件作为监测数据的导航文件，并根据基数据的导航文件和监测数据的导航文件，确定采集时的接收点位置误差。
进一步地，所述根据勘探地区的实际采集数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，包括：
a)读取勘探地区的一实际采集数据的导航文件，解析该导航文件的文件号、炮点号、炮点坐标和对应的接收点坐标，得到解析记录；
b)从解析记录中按照顺序提取出一炮的观测系统记录，以各电缆首末接收点位置为端点，线性插值中间等间距布放接收点的位置坐标，生成与基数据电缆保持相同羽角且呈直线展布的接收电缆；
c)重复所述步骤b)，完成该导航文件全部炮观测系统记录的接收点位置转换；
d)保持原炮点坐标不变，按照标准导航文件格式输出转换后的直电缆采集观测系统；
e)重复所述步骤a)～步骤d)，直至完成勘探地区的全部实际采集数据的导航文件的生成。
进一步地，所述将实际采集数据的导航文件作为基数据的导航文件，将生成的导航文件作为监测数据的导航文件，并根据基数据的导航文件和监测数据的导航文件，确定采集时的接收点位置误差，包括：
将实际采集数据的导航文件作为基数据的导航文件，将生成的炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件作为监测数据的导航文件；
针对每一接收点，分别计算基数据的导航文件和监测数据的导航文件之间的接收点位置误差；
根据每一接收点位置误差，计算全部接收点的平均位置误差；
针对每一炮点，分别计算接收点位置平均误差。
进一步地，所述建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
建立勘探地区的三维地震速度模型；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据。
进一步地，所述针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件，模拟时移地震三维拖缆采集基数据；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据监测数据的导航文件，模拟时移地震三维拖缆采集监测数据。
进一步地，所述根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性，包括：
根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，计算勘探地区储层上方标志层位的NRMS值；
根据计算的NRMS值，分析勘探地区的基数据和监测数据的一致性，若NRMS值小于0.1，则时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据的一致性较好；否则，时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据的一致性较差。
一种时移地震三维拖缆采集数据的可重复性分析系统，包括：
导航文件生成模块，用于根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件；
合成地震数据生成模块，用于建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据；
数据处理模块，用于对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体；
分析模块，用于根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性。
一种处理器，包括计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现上述时移地震三维拖缆采集数据的可重复性分析方法对应的步骤。
一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现上述时移地震三维拖缆采集数据的可重复性分析方法对应的步骤。
本发明由于采取以上技术方案，其具有以下优点：
1、本发明采用射线追踪数值算法定量分析时移地震三维拖缆采集电缆漂移带来的数据的一致性影响，可以为海上时移地震三维拖缆采集前和采集后的分析提供可靠依据。
2、本发明根据实际采集数据的P190导航文件生成监测数据采集观测系统，结合建立的三维地震速度模型，解决了只有一次实际采集数据的时移地震一致性分析难题，可以广泛应用于石油天然气地震勘探领域中。
附图说明
图1是本发明一实施例提供的方法流程图；
图2是本发明一实施例提供的羽角结构示意图；
图3是本发明一实施例提供的实际采集数据电缆羽角实时变化示意图，其中，横坐标为炮号，纵坐标为羽角(度)；
图4是本发明一实施例提供的L油田实际采集时的电缆展布示意图；
图5是本发明一实施例提供的与图4对应炮点的直电缆展布示意图；
图6是本发明一实施例提供的L油田实际采集数据的接收点位置误差分布示意图，其中，横坐标为接收点位置误差，纵坐标为接收点数；
图7是本发明一实施例提供的L油田四个相邻采集工区的接收点平均位置误差示意图；
图8是本发明一实施例提供的L油田单炮接收点平均位置误差平面分布示意图，其中，横坐标为大地横坐标X，纵坐标为大地纵坐标Y；
图9是本发明一实施例提供的弯曲电缆采集的基数据的叠前深度偏移数据体示意图；
图10是本发明一实施例提供的直电缆采集的监测数据的叠前深度偏移数据体示意图；
图11是本发明一实施例提供的弯曲电缆采集基数据与直电缆采集监测数据的一致性对比示意图，其中，横坐标为主测线号，纵坐标为时间(秒)；
图12是本发明一实施例提供的弯曲电缆采集基数据与冗余方法采集监测数据的一致性对比示意图，其中，横坐标为主测线号，纵坐标为时间(秒)。
具体实施方式
下面将参照附图更详细地描述本发明的示例性实施方式。虽然附图中显示了本发明的示例性实施方式，然而应当理解，可以以各种形式实现本发明而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了能够更透彻地理解本发明，并且能够将本发明的范围完整的传达给本领域的技术人员。
应理解的是，文中使用的术语仅出于描述特定示例实施方式的目的，而无意于进行限制。除非上下文另外明确地指出，否则如文中使用的单数形式“一”、“一个”以及“所述”也可以表示包括复数形式。术语“包括”、“包含”、“含有”以及“具有”是包含性的，并且因此指明所陈述的特征、步骤、操作、元件和/或部件的存在，但并不排除存在或者添加一个或多个其它特征、步骤、操作、元件、部件、和/或它们的组合。文中描述的方法步骤、过程、以及操作不解释为必须要求它们以所描述或说明的特定顺序执行，除非明确指出执行顺序。还应当理解，可以使用另外或者替代的步骤。
本发明实施例提供的时移地震三维拖缆采集数据的可重复性分析方法及系统，针对海上时移地震三维拖缆采集的可重复性问题。根据实际采集数据的导航文件构造出保持相同羽角且呈拉直状态的接收电缆采集位置，采用射线追踪数值模拟基数据(真实位置)采集和监测数据(直电缆)采集，通过计算模拟数据标志层位的NRMS值，推断实际采集数据的一致性，为时移地震拖缆采集评价提供可靠依据，解决了只有一次实际采集数据的时移地震一致性分析难题。
实施例1
如图1所示，本实施例提供一种时移地震三维拖缆采集数据的可重复性分析方法，包括以下步骤：
1)根据勘探地区的实际采集数据(即时移地震采集基数据)的P190导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的P190导航文件，其中，P190导航文件是英国勘探委员会UKOOA于20世纪90年代初制定的地震勘探导航格式文件，作为行业标准广泛应用于海上拖缆地震采集定位中，P190导航文件逐炮记录拖缆采集时的炮点和接收点位置信息；羽角(feather angle)是指接收电缆首尾连线与设计测线方位之间的夹角(如图2所示)，海上拖缆采集过程中，受风浪、洋流的影响，电缆羽角是随时变化的，采集电缆也时常呈曲线状态展布。
建立P190导航文件的具体过程为：
1.1)读取勘探地区的一实际采集数据的P190导航文件，解析该P190导航文件的文件号、炮点号、炮点坐标和对应的接收点坐标，得到解析记录。
1.2)从解析记录中按照顺序提取出一炮的观测系统记录，以各电缆首末接收点位置为端点，线性插值中间等间距布放接收点的位置坐标，生成与基数据电缆保持相同羽角且呈直线展布的接收电缆：
海上时移地震三维拖缆采集施工要求选择相同的时间窗口实施，以保证洋流、季风等环境影响因素尽量可以重复，在航线和航向相同的情况下，监测数据采集时电缆羽角与基数据采集时的电缆羽角大致相同，且对电缆的姿态控制精度也越来越高，因此将采集监测数据的观测系统设置为与基数据采集保持相同羽角的直电缆是符合实际情况的。
1.3)重复步骤1.2)，完成该P190导航文件全部炮观测系统记录的接收点位置转换。
1.4)保持原炮点坐标不变，按照标准P190导航文件格式输出转换后的直电缆采集观测系统，其中，海上拖缆采集中对炮点的控制相对容易，炮点位置误差相对于电缆漂移产生的接收点位置误差小很多，可以忽略。
1.5)重复步骤1.1)～步骤1.4)，直至完成勘探地区的全部实际采集数据的P190导航文件的生成。
2)将实际采集数据的P190导航文件作为基数据的P190导航文件，将步骤1)中生成的P190导航文件作为监测数据的P190导航文件，并根据基数据的P190导航文件和监测数据的P190导航文件，确定采集时的接收点位置误差，具体为：
2.1)将实际采集数据的P190导航文件作为基数据的P190导航文件，将步骤1)中生成的炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的P190导航文件作为监测数据的P190导航文件。
2.2)针对每一接收点，分别计算基数据的P190导航文件和监测数据的P190导航文件之间的接收点位置误差，即计算同一炮的基数据的P190导航文件中接收点与监测数据的P190导航文件中对应接收点之间的距离，并绘制接收点位置误差分布直方图：
由于接收点位置误差分布的数量级差别可能很大，直方图的纵坐标可以采用对数坐标。
2.3)根据每一接收点位置误差，计算全部接收点的平均位置误差，它是衡量拖缆采集可重复性的重要指标。
2.4)针对每一炮点，分别计算接收点位置平均误差，并绘制位置误差平面分布分布图：
每一炮激发均由多条电缆上的若干接收器同时接收，因此，首先计算该炮激发时每一接收点的基数据的P190导航文件和监测数据的P190导航文件之间的位置误差，然后计算该炮的全部接收点位置误差的算数平均，得到该炮点的接收点位置平均误差，再重复上述步骤直至计算得到所有炮点的接收点位置平均误差。
上述步骤2.2)至2.4)从不同侧面详细反映基数据采集时的可重复性误差，能够对比监测数据采集后与基数据采集时的一致性改善和提高。
3)建立勘探地区的三维地震速度模型，具体为：
在勘探地区的地震解释成果的层位解释数据的基础上，结合岩性解释数据、测井数据和时深关系，建立勘探地区的深度域三维速度体。
当没有勘探地区的地震解释成果作为参考时，可以采用深度域或时间域的地震成像结果作为勘探地区的构造模型，并采用偏移速度填充为三维地震速度模型。时间域均方根速度转深度域层速度出现异常结果时，应简化勘探地区的构造模型并适度平滑填充速度，得到三维地震速度模型。
4)针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的P190导航文件和监测数据的P190导航文件，模拟时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据，得到合成地震数据，模拟计算中，除P190导航文件指定的接收点位置不同以外，其它模拟参数均保持一致，具体为：
4.1)针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的P190导航文件即步骤1)中实际采集数据的P190导航文件描述的弯曲电缆观测系统，模拟时移地震三维拖缆采集的基数据。
4.2)针对建立的三维地震速度模型，采用三维射线追踪算法，根据监测数据的P190导航文件即步骤1)中生成的P190导航文件描述的直电缆观测系统，模拟时移地震三维拖缆采集的监测数据。
5)采用相同的处理流程和相同的处理参数，对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体，其中，处理流程包括直达波切除、几何扩散补充、道集分选和叠前深度偏移等处理，处理参数包括切除线、能量补充因子、数据域、反假频算子和偏移孔径等参数。
6)根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，计算勘探地区储层上方标志层位的NRMS值，并分析勘探地区的基数据和监测数据的一致性，具体为：
勘探地区储层上方的标志层位不受储层物性变化影响，计算基数据成像剖面和监测数据成像剖面标志层位的NRMS值，评估两次采集地震数据的一致性，NRMS值是监测数据与基数据的差的平均均方根振幅除以两数据的平均均方根振幅和，即：
其中，B为时移地震三维拖缆采集的基数据，M为时移地震三维拖缆采集的监测数据。rms算子定义为：
其中，xi为时窗内的振幅，N为时窗里的样点个数。
NRMS值受相位和振幅差、时移误差和噪音的影响，值越小则说明时移地震三维拖缆采集基数据B和时移地震三维拖缆采集监测数据M的一致性越好。实际生产中，一般以储层上方标志层位的NRMS值小于0.1作为评判标准，NRMS值小于0.1，则说明时移地震三维拖缆采集基数据B和时移地震三维拖缆采集监测数据M的一致性较好，采集位置误差带来的数据一致性影响相对储层物性变化引起的数据差异很小；否则，则说明时移地震三维拖缆采集基数据B和时移地震三维拖缆采集监测数据M的一致性较差，会影响对储层物性变化的正确判断。
下面以中国南海L油田为具体实施例详细说明本发明的时移地震三维拖缆采集数据的可重复性分析方法，其中，该数据于2011年采用双源三缆采集(缆长5100米)：
1)根据L油田的基数据的P190导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的P190导航文件：
P190导航文件逐炮记录了拖缆采集时的炮点和接收点位置信息。受风浪、洋流的影响，海上拖缆采集电缆羽角实时变化，采集电缆常常呈曲线状态展布。如图3所示，展示了L油田某实际炮线的三条采集电缆的羽角随炮点的变化；如图4所示，为某个炮点激发时三条电缆的平面展布，它们与炮线方向形成羽角并发生弯曲；如图5所示，展示了按照步骤1)的内容生成的与图4同一处炮点激发的直电缆展布，它与图4中的三条电缆保持一样的羽角。
2)将实际采集数据的P190导航文件作为基数据的P190导航文件，将步骤1)生成的P190导航文件作为监测数据的P190导航文件，并根据基数据的P190导航文件和监测数据的P190导航文件，确定采集时的接收点位置误差：
如图6所示，展示了L油田实际采集数据的接收点位置误差分布(包含193523次放炮，每次放炮记录1224个接收点位置，超2.3亿个位置误差)，此次实际采集数据的全部接收点平均位置误差为33.64m。如图7所示，为根据与L油田相邻海域的另外四次拖缆采集作业数据统计的全部接收点平均位置误差，其中最小的为21.48m，最大的为26.98m，都比33.64m小，由此可以大致判断它们采集施工时对电缆的姿态控制精度均比L油田2011年采集时控制得更好。
此外，还针对每一炮点，分别计算接收点位置平均误差，绘制位置误差平面分布分布图，如图8所示，为L油田实际采集数据按照炮点位置标注的各炮接收点平均位置误差，突出了接收点位置误差的平面分布和局部特征。上述这些信息有助于了解本次拖缆采集地震数据的可重复性的基本情况。
3)建立勘探地区的三维地震速度模型。
4)针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的P190导航文件和监测数据的P190导航文件，模拟时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据，得到合成地震数据，其中，本实施例采用商业射线追踪模拟软件NOSAR数值模拟地震采集。
5)采用相同的处理流程和相同的处理参数，对合成地震数据进行处理，输出L油田的叠前深度偏移数据体：
步骤4)的合成地震数据不含噪音和多次波，因此只需完成直达波切除和几何扩散补偿后即可结合建立的三维速度模型进行叠前深度偏移处理。如图9和图10所示，为经相同处理流程和处理参数得到的弯曲电缆采集的基数据和直电缆采集的监测数据的叠前深度偏移数据体。仅从inline线和xline线剖面上难以看出基数据和监测数据之间的细微差别，还需要进一步计算二者的NRMS值。
6)根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，计算储层上方标志层位的NRMS值，并分析勘探地区的基数据和监测数据的一致性：
如图11所示，为沿L油田一条xline线截取的叠前深度偏移成像数据体剖面对比，从差剖面上能够看到较强的残留能量，计算该剖面基数据和监测数据的NRMS值为0.27，超过门槛值0.1，这说明此次拖缆实采数据的可重复性较差，如果直接用作时移地震数据分析，很可能会影响对储层物性变化的正确判断，进而导致剩余油预测失败。因此，本发明的方法能够精确确定时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据的一致性，为海上时移地震三维拖缆采集设计和采集后评估提供可靠依据。
如图12所示，为采用一种提高时移地震数据可重复性的采集方法(其原理是利用监测数据的冗余逼近可重复性差的基数据)采集的监测数据，它与原来的基数据的NRMS值降至0.067，低于门槛值0.1，表明因采集施工产生的位置误差带来的数据一致性影响降至可以接受的水平，有效提高了基数据和监测数据的一致性。
实施例2
本实施例提供一种时移地震三维拖缆采集数据的可重复性分析系统，包括：
导航文件生成模块，用于根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件。
合成地震数据生成模块，用于建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据。
数据处理模块，用于对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体。
分析模块，用于根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性。
实施例3
本实施例提供一种与本实施例1所提供的时移地震三维拖缆采集数据的可重复性分析方法对应的处理设备，处理设备可以是用于客户端的处理设备，例如手机、笔记本电脑、平板电脑、台式机电脑等，以执行实施例1的方法。
所述处理设备包括处理器、存储器、通信接口和总线，处理器、存储器和通信接口通过总线连接，以完成相互间的通信。存储器中存储有可在处理器上运行的计算机程序，处理器运行计算机程序时执行本实施例1所提供的时移地震三维拖缆采集数据的可重复性分析方法。
在一些实现中，存储器可以是高速随机存取存储器(RAM：Random AccessMemory)，也可能还包括非不稳定的存储器(non-volatile memory)，例如至少一个磁盘存储器。
在另一些实现中，处理器可以为中央处理器(CPU)、数字信号处理器(DSP)等各种类型通用处理器，在此不做限定。
实施例4
本实施例1的时移地震三维拖缆采集数据的可重复性分析方法可被具体实现为一种计算机程序产品，计算机程序产品可以包括计算机可读存储介质，其上载有用于执行本实施例1所述的时移地震三维拖缆采集数据的可重复性分析方法的计算机可读程序指令。
计算机可读存储介质可以是保持和存储由指令执行设备使用的指令的有形设备。计算机可读存储介质例如可以是但不限于电存储设备、磁存储设备、光存储设备、电磁存储设备、半导体存储设备或者上述的任意组合。
上述各实施例仅用于说明本发明，其中各部件的结构、连接方式和制作工艺等都是可以有所变化的，凡是在本发明技术方案的基础上进行的等同变换和改进，均不应排除在本发明的保护范围之外。
Claims (8)
1.一种时移地震三维拖缆采集数据的可重复性分析方法，其特征在于，包括：
根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件；
建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
建立勘探地区的三维地震速度模型；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件，模拟时移地震三维拖缆采集基数据；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据监测数据的导航文件，模拟时移地震三维拖缆采集监测数据；
对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体；
根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性。
2.如权利要求1所述的一种时移地震三维拖缆采集数据的可重复性分析方法，其特征在于，所述根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件，包括：
根据勘探地区的实际采集数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件；
将实际采集数据的导航文件作为基数据的导航文件，将生成的导航文件作为监测数据的导航文件，并根据基数据的导航文件和监测数据的导航文件，确定采集时的接收点位置误差。
3.如权利要求2所述的一种时移地震三维拖缆采集数据的可重复性分析方法，其特征在于，所述根据勘探地区的实际采集数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，包括：
a)读取勘探地区的一实际采集数据的导航文件，解析该导航文件的文件号、炮点号、炮点坐标和对应的接收点坐标，得到解析记录；
b)从解析记录中按照顺序提取出一炮的观测系统记录，以各电缆首末接收点位置为端点，线性插值中间等间距布放接收点的位置坐标，生成与基数据电缆保持相同羽角且呈直线展布的接收电缆；
c)重复所述步骤b)，完成该导航文件全部炮观测系统记录的接收点位置转换；
d)保持原炮点坐标不变，按照标准导航文件格式输出转换后的直电缆采集观测系统；
e)重复所述步骤a)～步骤d)，直至完成勘探地区的全部实际采集数据的导航文件的生成。
4.如权利要求2所述的一种时移地震三维拖缆采集数据的可重复性分析方法，其特征在于，所述将实际采集数据的导航文件作为基数据的导航文件，将生成的导航文件作为监测数据的导航文件，并根据基数据的导航文件和监测数据的导航文件，确定采集时的接收点位置误差，包括：
将实际采集数据的导航文件作为基数据的导航文件，将生成的炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件作为监测数据的导航文件；
针对每一接收点，分别计算基数据的导航文件和监测数据的导航文件之间的接收点位置误差；
根据每一接收点位置误差，计算全部接收点的平均位置误差；
针对每一炮点，分别计算接收点位置平均误差。
5.如权利要求1所述的一种时移地震三维拖缆采集数据的可重复性分析方法，其特征在于，所述根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性，包括：
根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，计算勘探地区储层上方标志层位的NRMS值；
根据计算的NRMS值，分析勘探地区的基数据和监测数据的一致性，若NRMS值小于0.1，则时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据的一致性较好；否则，时移地震三维拖缆采集基数据和时移地震三维拖缆采集监测数据的一致性较差。
6.一种时移地震三维拖缆采集数据的可重复性分析系统，其特征在于，包括：
导航文件生成模块，用于根据勘探地区的时移地震采集基数据的导航文件，生成炮点位置与基数据相同、接收电缆与基数据电缆保持相同羽角且呈直线状态展布的导航文件，并将其作为时移地震采集监测数据的导航文件；
合成地震数据生成模块，用于建立勘探地区的三维地震速度模型，并根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
建立勘探地区的三维地震速度模型；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件和监测数据的导航文件，模拟基数据和监测数据，得到合成地震数据，包括：
针对建立的三维地震速度模型，采用三维射线追踪算法，根据基数据的导航文件，模拟时移地震三维拖缆采集基数据；
针对建立的三维地震速度模型，采用三维射线追踪算法，根据监测数据的导航文件，模拟时移地震三维拖缆采集监测数据；
数据处理模块，用于对合成地震数据进行处理，得到勘探地区的叠前深度偏移数据体；
分析模块，用于根据勘探地区的叠前深度偏移数据体，对比基数据和监测数据的成像剖面，分析勘探地区的基数据和监测数据的一致性。
7.一种处理器，其特征在于，包括计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现权利要求1-5中任一项所述的时移地震三维拖缆采集数据的可重复性分析方法对应的步骤。
8.一种计算机可读存储介质，其特征在于，所述计算机可读存储介质上存储有计算机程序指令，其中，所述计算机程序指令被处理器执行时用于实现权利要求1-5中任一项所述的时移地震三维拖缆采集数据的可重复性分析方法对应的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110820603.6A CN113568041B (zh) | 2021-07-20 | 2021-07-20 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110820603.6A CN113568041B (zh) | 2021-07-20 | 2021-07-20 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113568041A CN113568041A (zh) | 2021-10-29 |
| CN113568041B true CN113568041B (zh) | 2024-01-30 |
Family
ID=78165809
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110820603.6A Active CN113568041B (zh) | 2021-07-20 | 2021-07-20 | 时移地震三维拖缆采集数据的可重复性分析方法及系统 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113568041B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115510562B (zh) * | 2022-09-30 | 2023-10-24 | 中国海洋大学 | 拖曳式传感器阵列系统的深度预测方法 |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005096018A1 (en) * | 2004-03-17 | 2005-10-13 | Westerngeco Seismic Holdings Ltd. | Marine seismic survey method and system |
| CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
| CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
| CN103645503A (zh) * | 2013-12-17 | 2014-03-19 | 中国海洋石油总公司 | 一种三维时间域照明分析及振幅补偿方法 |
| WO2014207229A2 (en) * | 2013-06-28 | 2014-12-31 | Cgg Services Sa | System and method for estimating repeatability using base data |
| CN104570068A (zh) * | 2013-10-11 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于三维gis的观测系统crp属性分析方法 |
| CN105259570A (zh) * | 2015-11-12 | 2016-01-20 | 中国海洋石油总公司 | 一种海上时移地震接收点位移校正方法 |
| CN105301650A (zh) * | 2015-10-09 | 2016-02-03 | 中国石油天然气集团公司 | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 |
| CN105589102A (zh) * | 2016-02-25 | 2016-05-18 | 中国海洋石油总公司 | 一种海上时移地震记录的潮汐校正方法 |
| CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
| CN106443768A (zh) * | 2016-12-14 | 2017-02-22 | 成都理工大学 | 一种叠前深度域合成地震记录制作方法 |
| CN106908841A (zh) * | 2017-01-17 | 2017-06-30 | 中国海洋石油总公司 | 海上非重复性拖缆时移地震数据高密度采集方法 |
| WO2018140783A1 (en) * | 2017-01-27 | 2018-08-02 | Saudi Arabian Oil Company | Virtual source redatuming using radiation pattern correction |
| CA3056450A1 (en) * | 2017-03-16 | 2018-09-20 | Saudi Arabian Oil Company | Continuous seismic reservoir monitoring using a common focus point method |
| CN109212589A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于gpu/cpu协同并行的地震采集观测系统设计方法 |
| CN110058309A (zh) * | 2018-12-13 | 2019-07-26 | 广州海洋地质调查局 | 一种立体时移地震观测系统及方法 |
| CN112505782A (zh) * | 2020-12-23 | 2021-03-16 | 中国石油大学(北京) | 四维地震中辐射模式校正的干涉基准面重建方法及系统 |
| CN112630825A (zh) * | 2020-12-02 | 2021-04-09 | 中国海洋大学 | 共炮检距域Beam叠前时间偏移成像方法、系统、介质及应用 |
Family Cites Families (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US9116257B2 (en) * | 2012-05-07 | 2015-08-25 | Fugro Norway As | Method and apparatus for sea current aided, enhanced azimuth seismic data acquisition |
| US20140200854A1 (en) * | 2013-01-15 | 2014-07-17 | Westerngeco L.L.C. | Methods and computing systems for survey design and modeling workflow for towed multimeasurement seismic streamer data |
| US10768325B2 (en) * | 2016-10-27 | 2020-09-08 | Exxonmobil Upstream Research Company | Method to estimate 4D seismic acquisition repeatability specifications from high-resolution near-water-bottom seismic images |
-
2021
- 2021-07-20 CN CN202110820603.6A patent/CN113568041B/zh active Active
Patent Citations (20)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2005096018A1 (en) * | 2004-03-17 | 2005-10-13 | Westerngeco Seismic Holdings Ltd. | Marine seismic survey method and system |
| CN1947032A (zh) * | 2004-03-17 | 2007-04-11 | 维斯特恩格科地震控股有限公司 | 海上地震测量方法和系统 |
| CN101013161A (zh) * | 2007-01-15 | 2007-08-08 | 中国石油大港油田勘探开发研究院 | 基于叠前波场模拟的地震勘探层位标定方法 |
| CN101315427A (zh) * | 2007-05-29 | 2008-12-03 | 中国石油天然气集团公司 | 一种复杂地区地震勘探资料处理的方法及系统 |
| WO2014207229A2 (en) * | 2013-06-28 | 2014-12-31 | Cgg Services Sa | System and method for estimating repeatability using base data |
| CN104570068A (zh) * | 2013-10-11 | 2015-04-29 | 中国石油化工股份有限公司 | 一种基于三维gis的观测系统crp属性分析方法 |
| CN103645503A (zh) * | 2013-12-17 | 2014-03-19 | 中国海洋石油总公司 | 一种三维时间域照明分析及振幅补偿方法 |
| CN105301650A (zh) * | 2015-10-09 | 2016-02-03 | 中国石油天然气集团公司 | 一种海上拖缆采集的时移地震数据的质量监测方法和装置 |
| CN105259570A (zh) * | 2015-11-12 | 2016-01-20 | 中国海洋石油总公司 | 一种海上时移地震接收点位移校正方法 |
| CN105589102A (zh) * | 2016-02-25 | 2016-05-18 | 中国海洋石油总公司 | 一种海上时移地震记录的潮汐校正方法 |
| CN106291709A (zh) * | 2016-07-20 | 2017-01-04 | 中国海洋石油总公司 | 一种海上拖缆宽频宽方位地震勘探方法 |
| CN106443768A (zh) * | 2016-12-14 | 2017-02-22 | 成都理工大学 | 一种叠前深度域合成地震记录制作方法 |
| CN106908841A (zh) * | 2017-01-17 | 2017-06-30 | 中国海洋石油总公司 | 海上非重复性拖缆时移地震数据高密度采集方法 |
| WO2018140783A1 (en) * | 2017-01-27 | 2018-08-02 | Saudi Arabian Oil Company | Virtual source redatuming using radiation pattern correction |
| CA3056450A1 (en) * | 2017-03-16 | 2018-09-20 | Saudi Arabian Oil Company | Continuous seismic reservoir monitoring using a common focus point method |
| CN110573911A (zh) * | 2017-03-16 | 2019-12-13 | 沙特阿拉伯石油公司 | 使用共聚焦点法的连续地震储层监测 |
| CN109212589A (zh) * | 2017-07-06 | 2019-01-15 | 中国石油化工股份有限公司 | 一种基于gpu/cpu协同并行的地震采集观测系统设计方法 |
| CN110058309A (zh) * | 2018-12-13 | 2019-07-26 | 广州海洋地质调查局 | 一种立体时移地震观测系统及方法 |
| CN112630825A (zh) * | 2020-12-02 | 2021-04-09 | 中国海洋大学 | 共炮检距域Beam叠前时间偏移成像方法、系统、介质及应用 |
| CN112505782A (zh) * | 2020-12-23 | 2021-03-16 | 中国石油大学(北京) | 四维地震中辐射模式校正的干涉基准面重建方法及系统 |
Non-Patent Citations (6)
| Title |
|---|
| Offshore time-lapse seismic acquisition design and repeatability analysis in L oilfield;Dongchuan X;《SEG 2020 Workshop: Broadband and Wide-azimuth Deepwater Seismic Technology》;第85页Method和Example部分 * |
| Repeatability analysis of land time‐lapse seismic data: CO2CRC Otway pilot project case study;Pevzner, R., Shulakova, V., Kepic, A., & Urosevic, M;《Geophysical prospecting》;第59卷(第1期);66-77 * |
| Time-lapse Seismic Receiver Motion Correction Method and Application;Xue, D., & Zhu, Z;《International Geophysical Conference, Qingdao, China》;729-732 * |
| 合兴场三分量三维地震资料采集观测系统设计;石全科;;《内蒙古石油化工》(第01期);142-144 * |
| 海上拖缆采集地震数据接收点位移校正;薛东川;朱振宇;王小六;何洋洋;;《石油地球物理勘探》;第52卷(第01期);1-4 * |
| 海洋地震垂直缆观测系统的射线照明分析;何勇等;《CT理论研究与应用》;第24卷(第5期);第690-700页正文第1-3节 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113568041A (zh) | 2021-10-29 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| RU2169931C2 (ru) | Способ и устройство для обработки сейсмического сигнала и проведения разведки полезных ископаемых | |
| CN102027390B (zh) | 对地震数据的联合内插及消除虚反射 | |
| CA2964893C (en) | Structure tensor constrained tomographic velocity analysis | |
| CN112883564B (zh) | 一种基于随机森林的水体温度预测方法及预测系统 | |
| WO1997033184A1 (en) | Method for geophysical processing and interpretation using seismic trace difference for analysis and display | |
| US20080137480A1 (en) | Method of Building a Subsurface Velocity Model | |
| US8706462B2 (en) | System and method for providing a physical property model | |
| CN111239827B (zh) | 基于局部相似系数的三维地震数据多次波压制方法 | |
| CN111812730B (zh) | 一种用于滑坡探测的电阻率数据融合三维成像方法及系统 | |
| CN113568041B (zh) | 时移地震三维拖缆采集数据的可重复性分析方法及系统 | |
| EP1292848B1 (en) | Seismic survey system | |
| US20220236435A1 (en) | Low-Frequency Seismic Survey Design | |
| CN113703044A (zh) | 古河道宽度的校正方法、装置、电子设备及存储介质 | |
| CN109387872B (zh) | 表面多次波预测方法 | |
| CN113568042B (zh) | 一种提高时移地震数据可重复性的拖缆采集方法和系统 | |
| CN113777654B (zh) | 一种基于伴随状态法初至波走时层析的海水速度建模方法 | |
| CN115755177A (zh) | 一种海底硫化物矿体结构宽频带地震约束反演成像方法 | |
| Smith et al. | Using topographic slope distributions to infer seafloor patterns | |
| CN113568040B (zh) | 一种时移地震采集数据的可重复性分析方法及系统 | |
| CN111965705A (zh) | 地震单炮记录地质层位的标定方法、装置、设备及介质 | |
| US10067254B2 (en) | Removal of an estimated acquisition effect from a marine survey measurement | |
| EP2784553A2 (en) | Predicting sensitivity to positioning for seismic surveys | |
| CN112305594A (zh) | 非均质储层的油气分布确定方法及系统 | |
| CN111177886B (zh) | 基于物探数据分析的海上测距规划及土厚预测的方法 | |
| CN113484914B (zh) | 海上风浪一致性影响量板制作方法、系统、介质及设备 |
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 | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |