CN104297790B - 一种起伏海底地震观测***中炮检距分布的定量评估方法 - Google Patents

一种起伏海底地震观测***中炮检距分布的定量评估方法 Download PDF

Info

Publication number
CN104297790B
CN104297790B CN201410482562.4A CN201410482562A CN104297790B CN 104297790 B CN104297790 B CN 104297790B CN 201410482562 A CN201410482562 A CN 201410482562A CN 104297790 B CN104297790 B CN 104297790B
Authority
CN
China
Prior art keywords
geophone
fluctuating
bin
geophone offset
distribution
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
Application number
CN201410482562.4A
Other languages
English (en)
Other versions
CN104297790A (zh
Inventor
翁斌
符力耘
郝振江
糜芳
孙文博
王小六
张云鹏
刘志鹏
杨俊�
刘永江
朱振宇
江南森
薛东川
孙伟家
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
Original Assignee
Institute of Geology and Geophysics of CAS
China National Offshore Oil Corp CNOOC
CNOOC Research Institute Co Ltd
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 Institute of Geology and Geophysics of CAS, China National Offshore Oil Corp CNOOC, CNOOC Research Institute Co Ltd filed Critical Institute of Geology and Geophysics of CAS
Priority to CN201410482562.4A priority Critical patent/CN104297790B/zh
Publication of CN104297790A publication Critical patent/CN104297790A/zh
Application granted granted Critical
Publication of CN104297790B publication Critical patent/CN104297790B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种起伏海底地震观测***中炮检距分布的定量评估方法,其包括以下步骤:给定起伏海底地形下地震观测***中所有炮点和检波点的三维坐标信息,其具体包括炮点的三维空间坐标、检波点的三维空间坐标以及炮点与检波点之间的对应关系;根据各炮检对中点的X坐标和Y坐标,找出属于同一面元的所有炮点和检波点;计算起伏海底地形条件下某个面元所对应的炮检距分布的非均匀性系数;对于某一目的层的所有面元,分别计算各个面元所对应的炮检距分布的非均匀性系数;利用计算得到的炮检距分布的非均匀性系数绘制该目的层炮检距分布的非均匀性系数分布图,其中颜色越深的部分表示该面元处炮检距分布的非均匀性越大。

Description

一种起伏海底地震观测***中炮检距分布的定量评估方法
技术领域
本发明属于石油地震勘探三维观测***优化设计领域,特别是关于一种起伏海底地震观测***中炮检距分布的定量评估方法。
背景技术
海洋是地球上最大的油气资源宝库,随着科学技术的发展和国民生活水平的不断提高,人们对能源的需求量越来越大,同时在陆上油气资源日益减少的形势下,人类对海洋的依赖就更强,所以开发利用海洋油气资源是解决能源接替的根本途径之一。中国海岸线总长度约18000多公里,200海里专属经济区的管辖海域面积约300万平方公里。我国近海蕴藏着丰富的油气资源,据全国二次普查油气资源量测算,仅我国近海域的油气资源中石油和天然气分别达到245.6亿吨和43万亿立方米,具有极大的油气勘探潜力。
对于海底几百米至几千米深处的目标而言,地震勘探方法是目前国内外寻找石油、天然气的主要方法。如今海上三维地震数据采集方法正在飞速发展,一方面,现代地震勘探船所拖的电缆数量在迅速增加;另一方面,由于海上勘探工作条件以及海底的复杂性和特殊性,海洋地震勘探方法还存在许多理论方法和技术上的问题。由于海底地形地震条件和地下地质条件均很复杂,对地震勘探设备和地震勘探技术的要求高,因此许多常规的地震勘探技术无法适用。近20多年来,适应复杂海底地形的地震勘探装备不断研制成功,地球物理技术和计算技术也不断进步,这推动了复杂海底地形下地震勘探技术的迅速发展,为海上油气田的发现发挥了重要的作用。复杂海底三维地震数据采集是地震勘探的基础,野外采集资料的质量直接影响到勘探效果,而采集参数的准确设计是资料采集成功的关键。如今复杂海底下三维地震数据采集设计方法正在飞速发展。一方面,现代地震勘探的震源和接收器数量在迅速增加,但另一方面,由于复杂海底下地震勘探工作条件和海底地形的复杂性和特殊性,现有的基于水平海底地形条件下地震勘探方法还存在许多理论方法和技术上的问题。
在如图1所示的海上起伏海底地形条件下,地震勘探的炮点和检波点具有明显的高程差异,进而造成地震观测***的不规则炮检距分布。这种炮检距的不规则分布并非地震观测***本身的分布特征造成的,而是由起伏海底地形所造成的假象,炮检距的不规则分布直接影响后续地震成像的效果。然而,现有技术中缺乏对起伏海底地形下三维地震观测***炮检距分布进行有效评估的适当方法。
发明内容
针对上述问题,本发明的目的是提供一种起伏海底地震观测***中炮检距分布的定量评估方法,该方法能够克服现有起伏海底地形下地震观测***炮检距分析技术中无法定量评估复杂海底地形影响的缺陷,通过该方法得到炮检距分布的评价依据,从而能够帮助选择更适合于该起伏海底地形的地震观测***,为提高起伏海底地形下地震采集资料的保真度提供保障。
为实现上述目的,本发明采取以下技术方案:一种起伏海底地震观测***中炮检距分布的定量评估方法,其包括以下步骤:1)给定起伏海底地形下地震观测***中所有炮点和检波点的三维坐标信息,其具体包括炮点的三维空间坐标、检波点的三维空间坐标以及炮点与检波点之间的对应关系;2)根据各炮检对中点的X坐标和Y坐标,找出属于同一面元的所有炮点和检波点;3)根据某个面元中的炮点和检波点,计算起伏海底地形条件下该面元所对应的炮检距分布的非均匀性系数;4)对于某一目的层的所有面元,采用与步骤3)相同的方法分别计算各个面元所对应的炮检距分布的非均匀性系数;5)利用计算得到的炮检距分布的非均匀性系数绘制该目的层炮检距分布的非均匀性系数分布图,其中颜色越深的部分表示该面元处炮检距分布的非均匀性越大。
所述步骤2)中,找出属于同一面元的所有炮点和检波点的过程为:如果一个炮检对中点的X坐标和Y坐标位于某个面元内,则该炮检对的三维坐标信息属于该面元。
所述步骤3)中,计算起伏海底地形条件下某个面元所对应的炮检距分布的非均匀性系数,其包括:首先,将炮点在起伏海底上的对应点关于目的层的镜像点与检波点之间的距离作为起伏海底地形下的炮检距,则起伏海底地形下的炮检距为:
R = ( x r - x s ) 2 + ( y r - y s ) 2 + ( z r - z s ′ ) 2 - - - ( 1 )
式中,(xs,ys,zs)为炮点在起伏海底上的对应点的三维坐标,(xs,ys,z's)为炮点在起伏海底上的对应点关于目的层的镜像点,(xr,yr,zr)为检波点的三维坐标;其次,根据式(1)计算得到某个面元所对应的炮检距R1,R2,…Rn;最后,根据某个面元所对应的炮检距R1,R2,…Rn,计算得到该面元所对应的炮检距分布的非均匀性系数为:
C ( R 1 , R 2 , . . . , R n ) = Σ i = 1 n Σ j = 1 , j ≠ i n 1 | R i - R j | - - - ( 2 )
式中,C表示该面元处炮检距分布的非均匀性,C值越小表示炮检距分布越均匀。
本发明由于采取以上技术方案,其具有以下优点:本发明由于通过起伏海底地形下的炮检距公式计算得到某个面元所对应的炮检距,并根据该面元所对应的炮检距计算得到该面元所对应的炮检距分布的非均匀性系数,进而得到某一目的层所有面元对应的炮检距分布的非均匀性系数,绘制分布图,因此本发明能够克服现有起伏海底地形下地震观测***炮检距分析技术中无法定量评估复杂海底地形影响的缺陷,通过该方法得到炮检距分布的评价依据,从而能够帮助选择更适合于该起伏海底地形的地震观测***,为提高起伏海底地形下地震采集资料的保真度提供保障。基于以上优点,本发明可以广泛应用于炮检距分布的评估中。
附图说明
图1是起伏海底地形的示意图;
图2是本发明定量评估方法的流程图;
图3是起伏海底地形条件下炮检距分布的示意图;
图4是采用本发明得到的起伏海底地形下地震观测***中炮检距分布图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图2所示,本发明起伏海底地震观测***中炮检距分布的定量评估方法包括以下步骤:
1)给定起伏海底地形下地震观测***中所有炮点和检波点的三维坐标信息,其具体包括炮点的三维空间坐标、检波点的三维空间坐标以及炮点与检波点之间的对应关系。
2)根据各炮检对中点的X坐标和Y坐标,找出属于同一面元的所有炮点和检波点;其具体过程为:如果一个炮检对中点的X坐标和Y坐标位于某个面元内,则该炮检对的三维坐标信息属于该面元。
3)根据某个面元中的炮点和检波点,计算起伏海底地形条件下该面元所对应的炮检距分布的非均匀性系数;
为了剥离地下介质因素的影响,独立的分析起伏海底地形对面元振幅能量均衡性的影响,对地下介质情况进行简化,假设地下介质为均匀介质,即介质速度为一恒定值,而且仅考虑水平目的层位的情况。
首先,如图3所示,起伏海底地形条件下地震波沿箭头方向传播,将炮点在起伏海底上的对应点关于目的层的镜像点与检波点之间的距离作为起伏海底地形下的炮检距,则起伏海底地形下的炮检距为:
R = ( x r - x s ) 2 + ( y r - y s ) 2 + ( z r - z s ′ ) 2 - - - ( 1 )
式中,(xs,ys,zs)为炮点在起伏海底上的对应点的三维坐标,(xs,ys,z's)为炮点在起伏海底上的对应点关于目的层的镜像点,(xr,yr,zr)为检波点的三维坐标。
其次,根据式(1)计算得到某个面元所对应的炮检距R1,R2,…Rn
最后,根据某个面元所对应的炮检距R1,R2,…Rn,计算得到该面元所对应的炮检距分布的非均匀性系数为:
C ( R 1 , R 2 , . . . , R n ) = Σ i = 1 n Σ j = 1 , j ≠ i n 1 | R i - R j | - - - ( 2 )
式中,C表示该面元处炮检距分布的非均匀性,C值越小表示炮检距分布越均匀。
4)对于某一目的层的所有面元,采用与步骤3)相同的方法分别计算各个面元所对应的炮检距分布的非均匀性系数。
5)如图4所示,利用计算得到的炮检距分布的非均匀性系数绘制该目的层炮检距分布的非均匀性系数分布图。炮检距分布的非均匀性系数分布图中,颜色越深的部分表示该面元处炮检距分布的非均匀性越大。
通过对起伏海底地形下三维地震观测***炮检距分布的定量分析,得到表示炮检距分布的非均匀性的具体数值,并将其作为评价炮检距分布的依据,定量分析起伏海底地形对地震观测***的影响,从而帮助选择更适合于该起伏海底地形的地震观测***,为提高起伏海底地形下地震采集资料的保真度提供保障。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和方法步骤等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (2)

1.一种起伏海底地震观测***中炮检距分布的定量评估方法,其包括以下步骤:
1)给定起伏海底地形下地震观测***中所有炮点和检波点的三维坐标信息,其具体包括炮点的三维空间坐标、检波点的三维空间坐标以及炮点与检波点之间的对应关系;
2)根据各炮检对中点的X坐标和Y坐标,找出属于同一面元的所有炮点和检波点;
3)根据某个面元中的炮点和检波点,计算起伏海底地形条件下该面元所对应的炮检距分布的非均匀性系数;
计算起伏海底地形条件下某个面元所对应的炮检距分布的非均匀性系数,其包括:
首先,将炮点在起伏海底上的对应点关于目的层的镜像点与检波点之间的距离作为起伏海底地形下的炮检距,则起伏海底地形下的炮检距为:
R = ( x r - x s ) 2 + ( y r - y s ) 2 + ( z r - z s ′ ) 2 - - - ( 1 )
式中,(xs,ys,zs)为炮点在起伏海底上的对应点的三维坐标,(xs,ys,z's)为炮点在起伏海底上的对应点关于目的层的镜像点,(xr,yr,zr)为检波点的三维坐标;
其次,根据式(1)计算得到某个面元所对应的炮检距R1,R2,…Rn
最后,根据某个面元所对应的炮检距R1,R2,…Rn,计算得到该面元所对应的炮检距分布的非均匀性系数为:
C ( R 1 , R 2 , ... , R n ) = Σ i = 1 n Σ j = 1 , j ≠ i n 1 | R i - R j | - - - ( 2 )
式中,C表示该面元处炮检距分布的非均匀性,C值越小表示炮检距分布越均匀;
4)对于某一目的层的所有面元,采用与步骤3)相同的方法分别计算各个面元所对应的炮检距分布的非均匀性系数;
5)利用计算得到的炮检距分布的非均匀性系数绘制该目的层炮检距分布的非均匀性系数分布图,其中颜色越深的部分表示该面元处炮检距分布的非均匀性越大。
2.如权利要求1所述的一种起伏海底地震观测***中炮检距分布的定量评估方法,其特征在于:所述步骤2)中,找出属于同一面元的所有炮点和检波点的过程为:如果一个炮检对中点的X坐标和Y坐标位于某个面元内,则该炮检对的三维坐标信息属于该面元。
CN201410482562.4A 2014-09-19 2014-09-19 一种起伏海底地震观测***中炮检距分布的定量评估方法 Active CN104297790B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201410482562.4A CN104297790B (zh) 2014-09-19 2014-09-19 一种起伏海底地震观测***中炮检距分布的定量评估方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201410482562.4A CN104297790B (zh) 2014-09-19 2014-09-19 一种起伏海底地震观测***中炮检距分布的定量评估方法

Publications (2)

Publication Number Publication Date
CN104297790A CN104297790A (zh) 2015-01-21
CN104297790B true CN104297790B (zh) 2017-01-25

Family

ID=52317592

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201410482562.4A Active CN104297790B (zh) 2014-09-19 2014-09-19 一种起伏海底地震观测***中炮检距分布的定量评估方法

Country Status (1)

Country Link
CN (1) CN104297790B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107219550B (zh) * 2017-04-27 2019-03-29 广州威拓电子科技有限公司 一种imf-obs的数据采集提取方法
CN111489410B (zh) * 2019-01-25 2023-11-28 中国石油天然气集团有限公司 观测***炮检点数据绘制方法及装置
CN112444843B (zh) * 2019-08-30 2024-06-18 中国石油化工股份有限公司 一种观测***的调整方法
CN110579799B (zh) * 2019-09-04 2020-10-27 中国科学院地质与地球物理研究所 一种等旅行时间间隔的地震采集观测方法及***
CN113109864B (zh) * 2020-01-09 2023-07-25 中国石油天然气集团有限公司 地震勘探炮点变观方法及装置

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6026059A (en) * 1997-11-14 2000-02-15 Petroleum Geo-Services (Us), Inc. Method of creating common-offset/common-azimuth gathers in 3-D seismic surveys and method of conducting reflection attribute variation analysis
CN101173989A (zh) * 2006-11-03 2008-05-07 中国石油集团东方地球物理勘探有限责任公司 地震采集三维观测***定量分析方法
CN102023310A (zh) * 2009-09-19 2011-04-20 中国石油化工集团公司 三维观测***均匀性分析方法
CN103576213A (zh) * 2012-07-27 2014-02-12 中国石油天然气集团公司 一种炮检距属性均匀性的确定方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6026059A (en) * 1997-11-14 2000-02-15 Petroleum Geo-Services (Us), Inc. Method of creating common-offset/common-azimuth gathers in 3-D seismic surveys and method of conducting reflection attribute variation analysis
CN101173989A (zh) * 2006-11-03 2008-05-07 中国石油集团东方地球物理勘探有限责任公司 地震采集三维观测***定量分析方法
CN102023310A (zh) * 2009-09-19 2011-04-20 中国石油化工集团公司 三维观测***均匀性分析方法
CN103576213A (zh) * 2012-07-27 2014-02-12 中国石油天然气集团公司 一种炮检距属性均匀性的确定方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"三维观测***属性均匀性的定量分析";吴永国 等;《石油地球物理勘探》;20120630;第47卷(第3期);第361-365页 *
"三维观测***采样均匀性分析方法的研究";碗学俭 等;《中国石油勘探》;20091231(第6期);第42-48页 *
"炮检距属性的非均匀性系数分析";赵虎 等;《石油地球物理勘探》;20110228;第46卷(第1期);正文第1-3节,图2 *

Also Published As

Publication number Publication date
CN104297790A (zh) 2015-01-21

Similar Documents

Publication Publication Date Title
Giuffrida et al. Fracture simulation parameters of fractured reservoirs: Analogy with outcropping carbonates of the Inner Apulian Platform, southern Italy
CN104297790B (zh) 一种起伏海底地震观测***中炮检距分布的定量评估方法
CN107526101B (zh) 一种获取地震反射波的采集和处理方法
CN103529475B (zh) 一种识别和解释碳酸盐岩古岩溶储层三维结构的方法
CN102012524B (zh) 一种海上三维地震观测***羽状漂移定量评估方法
Anderson et al. An integrated approach to marine electromagnetic surveying using a towed streamer and source
CN102012521B (zh) 一种地震储层预测中叠前裂缝的检测方法
CN103645503B (zh) 一种三维时间域照明分析及振幅补偿方法
CN112883564B (zh) 一种基于随机森林的水体温度预测方法及预测***
CN104316965B (zh) 一种裂缝方位和强度的预测方法及***
CN105116453A (zh) 一种冻土带天然气水合物的瞬变电磁勘探方法及装置
Ali et al. Geophysical investigation of Al Jaww Plain, eastern Abu Dhabi: Implications for structure and evolution of the frontal fold belts of the Oman Mountains
Kortekaas et al. Improved definition of faults in the Groningen field using seismic attributes
Park et al. CSEM sensitivity study for Sleipner CO2-injection monitoring
Alcocer et al. Reducing uncertainty by integrating 3D CSEM in the Mexican deep-water exploration workflow
Peirce et al. Constraints on crustal structure of adjacent OCCs and segment boundaries at 13° N on the Mid-Atlantic Ridge
CN103293553A (zh) 一种复杂海底上下缆地震采集数据边界元延拓校正方法
CN103454681A (zh) 评价三维地震观测***成像效果的方法和设备
CN104570095A (zh) 一种基于Radon变换消除斜缆虚反射的方法
LONG et al. Seismic geothermal resource exploration based on CPU/GPU collaborative parallel prestack time migration
CN104977615A (zh) 一种基于模型统计拾取的深水obc资料多次波压制方法
CN104391332A (zh) 浅海双频可控源电磁勘探方法
CN103492909A (zh) 动态地震勘探方法
CN111158050B (zh) 数据采集***、方法及隧道地震波超前预报方法
Sacchi et al. The use and beauty of ultra-high-resolution seismic reflection imaging in Late Quaternary marine volcaniclastic settings, Napoli Bay, Italy

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CP01 Change in the name or title of a patent holder
CP01 Change in the name or title of a patent holder

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Co-patentee after: Institute of Geology and Geophysics, Chinese Academy of Sciences

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC Research Institute

Patentee before: China National Offshore Oil Corporation

Co-patentee before: Institute of Geology and Geophysics, Chinese Academy of Sciences

TR01 Transfer of patent right
TR01 Transfer of patent right

Effective date of registration: 20191210

Address after: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee after: CNOOC research institute limited liability company

Patentee after: China Offshore Oil Group Co., Ltd.

Address before: 100010 Beijing, Chaoyangmen, North Street, No. 25, No.

Co-patentee before: CNOOC research institute limited liability company

Patentee before: China Offshore Oil Group Co., Ltd.

Co-patentee before: Institute of Geology and Geophysics, Chinese Academy of Sciences