CN102236104B - 三维密度变速地震反射空校方法 - Google Patents
三维密度变速地震反射空校方法 Download PDFInfo
- Publication number
- CN102236104B CN102236104B CN201010169276.4A CN201010169276A CN102236104B CN 102236104 B CN102236104 B CN 102236104B CN 201010169276 A CN201010169276 A CN 201010169276A CN 102236104 B CN102236104 B CN 102236104B
- Authority
- CN
- China
- Prior art keywords
- velocity
- density
- seismic
- reflection
- dimensional
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
重力-地震勘探的数据处理中三维密度变速地震反射空校方法,利用地震速度数据体计算各反射层的厚度和底界海拔埋深,利用重力三维剩余密度体数据计算厚度和底界深度,对公共区域的重力三维剩余密度数据进行插值,各地震测点下的各反射层进行密度校正,根据新的层速度计算反射层厚度,按坐标得到每一地震层校正前后的反射层速度、平均速度和埋深,获得到经过三维密度变速地震反射空校的结果。本发明改进了地震速度的效果,提高了反射层速度的精度,解决了复杂区实际地质问题。
Description
技术领域
本发明属于地球物理勘探方法,涉及重力-地震勘探的数据处理技术,具体是一种三维密度变速地震反射空校方法。
背景技术
地震勘探是利用地下地层层速度差异引起相应的反射波特征。现行地震勘探中反射层速度观测和研究是基于地震勘探方法建立的,即通过地震反射能量速度扫描等技术实现,因而在地震资料品质差的地区,无法进行有效的速度扫描,观测困难,甚至造成速度失真。在一些山区或构造复杂地区,地震勘探的反射层的速度受种种因素的影响,精度降低,导致地质模型的畸变,其结果与实际的差异较大,无法满足石油勘探对复杂地区的需要。
重力勘探是油气勘探中的另一种方法,它是针对地下地层密度差异引起对应的重力异常进行测量和分析,用于研究地质构造。
目前采用的三维密度分布是重力勘探获得的三维反演结果,比以往的二维反演技术更进一步。但在以往的地震速度勘探中很少利用重力勘探获得的三维反演结果,也没有利用三维密度对地震的反射层速度进行校正,而提高重力资料利用价值和提高地震勘探成功率的技术。
发明内容
本发明的目的在于提供一种重力与地震联合勘探,提高反射层速度的精度,解决复杂区实际地质问题的三维密度变速地震反射空校方法。
本发明采用以下技术步骤实现:
1)野外采集的重力数据和地震数据,分别处理得到三维剩余密度数据体和地震速度数据体;
步骤1)所述的三维剩余密度数据体为坐标X、坐标Y、坐标Z、剩余密度Δσ,坐标Z为地下半空间某一剖分长方体单元的中心海拔,剩余密度Δσ为该剖分单元的相对正负值,单位为g/cm3。
步骤1)所述的地震速度数据体为坐标X、坐标Y、反射层名LName、反射层的反射时间T,单位为毫秒、反射层的层速度V、平均速度Va、反射层相对于统一基准面Height的底界埋深D。
2)利用地震速度数据体数据计算各反射层的厚度和底界海拔埋深;
步骤2)所述的地震速度数据体计算是确定数据体的坐标点数,每个坐标点的坐标及该位置下面的地震反射层数目,循环地震速度数据体。
步骤2)按照下式计算各层的厚度Thicki:
Thicki=Ti*Vi/2000;
按照下式计算各层的厚度相加得到各层的深度Di:
Di=(Thick1+Thick2+…+Thicki);
按照下式计算各层的海拔埋深Di:
Di=Height-Di;
式中:地震各层的反射时间为Ti,层速度为Vi,统一基准面为Height。
3)利用重力三维剩余密度体数据计算厚度和底界深度;
步骤3)所述的计算是首先确定数据体的记录数和求取数据体在空间三个方向上的循环数,再循环读取重力三维剩余密度体数据。
步骤3)所述的厚度为剖分单元的厚度GThick,底界深度为该剖分单的底界海拔GD,按照下式计算:
GThick1=Z1-Z2,GD1=Z1-(Z1-Z2)/2;GThicki=2*(GDi-1-Zi),GDi=Zi-GThicki/2
式中Z1,Z2,……,Zi,分别为三维剩余密度数据体中地下第一,第二,……,第i层剖分长方体单元的中心海拔。
4)对公共区域的重力三维剩余密度数据进行插值;
步骤4)所述的插值是依据地震速度数据体的坐标位置,对重力三维剩余密度数据体进行插值,计算出地震数据体所在位置的三维剩余密度。
步骤4)所述的插值为一元三次自然样条插值:先沿x方向以纵向坐标Y为变量对三维剩余密度数据Δσ进行插值得到Δσx,再沿y方向以横坐标X为变量对Δσx进行插值,得到地震数据体所在位置的三维剩余密度Δσx y。
5)对公共区域内的各地震测点下的各反射层进行密度校正;
步骤5)密度校正是首先由下式计算出层密度σ0:σ0=0.31*Vi 1/4;式中:Vi层速度;再用层密度σ0加重力三维剩余演密度Δσ(Δσ为地震反射层中的各重力剖分单元的重力三维剩余演密度Δσxy按厚度加权平均所得)得到新的层密度σ1;将σ1再代入上式得到校正后的新层速度NewVi(即:NewVi=(σ1/0.31)4)。
6)根据新的层速度按照下式计算反射层厚度NewThicki:
NewThicki=Ti*NewVi/2000;
按照下式计算反射层界面埋深NewDi:
NewDi=(NewThick1+NewThick2+…+NewThicki);
按照下式计算反射层的新的平均速度NewVa:
NewVa=NewDi*2/Ti;
7)按坐标得到每一地震层校正前后的反射层速度、平均速度和埋深,获得到经过三维密度变速地震反射空校的结果。
步骤7)所述的按坐标得到每一地震层校正前后的反射层速度、平均速度和埋深地震速度数据体为坐标X,坐标Y,反射层名LName,反射层的反射时间T、校正前反射层的层速度V、平均速度Va、底界埋深D,校正后反射层的层速度NewV、平均速度NewVa、底界埋深NewD。
本发明具有改进地震速度的效果,应用于重力三维剩余密度数据体对地震的反射层速度进行校正,提高了反射层速度的精度,解决了复杂区实际地质问题。
附图说明
图1本发明使用的原始构造区三维剩余密度分布图(单位为g/cm3);
图2本发明处理后的构造区变速等深对比图,a.地震速度变速等深图,b.三维密度变速等深图;
图3图2a与b的叠合对比图;
图4本发明处理后的构造区T8反射层平均速度对比图,a.地震T8反射层平均速度,b.三维密度变速后T8反射层平均速度,c.a与b的叠合对比图。
具体实施方式
本发明使用的原始三维剩余密度体,包括坐标X,坐标Y,坐标Z,剩余密度Δσ(如图1)和地震速度数据体,包括坐标X,坐标Y,反射层名LName、反射层时间T、反射层速度V、平均速度Va、反射层界面埋深H等数据。
1)使用地震速度数据体:即坐标X,坐标Y,反射层名LName、及各反射层时间T、反射层速度V、平均速度Va、反射层界面埋深H等数据,并计算各反射层的厚度D。
2)使用重力三维剩余密度体数据:坐标X,坐标Y,坐标Z,三维剩余密度Δσ,并计算重力三维剩余密度数据的底界深度GD、厚度GThick。
3)依据地震速度数据体的坐标位置,对该坐标点用一元三次自然样条法作剩余密度插值,计算出地震数据体所在位置的三维剩余密度Δσxy。
4)分别计算第一层到第六层各层的新层速度。先由地震的层速度(Vi,i=1、2、3、4、5、6)按σ0=0.31*Vi0.25公式计算出各层的层密度σ0;对地震层中各重力剖分单元的三维剩余密度Δσxy按厚度加权平均得到Δσ;新的层密度为σ1=σ0+Δσ;计算校正后的新层速度,即NewVi=(σ1/0.31)4,i分别等于1、2、3、4、5、6。
5)根据校正后的新层速度计算出反射层厚度Thicki=Ti*NewVi/2000,(i分别等于1、2、3、4、5、6),第六反射层新的底界埋深NewD6=(NewThick1+NewThick2+NewThick3+NewThick4+NewThick5+NewThick6),第六反射层的新的平均速度Va6=NewD6*2/T6。
6)得到各个坐标点的计算结果:经过校正后的第六反射层的底界埋深NewD6、平均速度Va6,见图2、图3、图4。
Claims (8)
1.一种三维密度变速地震反射空校方法,其特征在于采用以下技术步骤实现:
1)野外采集的重力数据和地震数据,分别处理得到重力三维剩余密度数据体和地震速度数据体;
2)利用地震速度数据体计算各反射层的厚度和海拔基准的底界深度;
3)利用重力三维剩余密度数据体计算厚度和底界深度;
所述的计算是首先确定数据体的记录数和求取数据体在空间三个方向上的循环数,再循环读取重力三维剩余密度数据体;
4)对公共区域的重力三维剩余密度数据体进行插值;
所述的插值是依据地震速度数据体的坐标位置,对重力三维剩余密度数据体进行插值,计算出地震速度数据体所在位置的三维剩余密度;
5)对公共区域内的各地震测点下的各反射层进行密度校正;
6)根据新的层速度按照下式计算反射层厚度NewThicki:
NewThicki=Ti*NewVi/2000;
按照下式计算反射层界面新的底界深度NewDi:
NewDi=(NewThick1+NewThick2+…+NewThicki);
按照下式计算反射层的新的平均速度NewVa:
NewVa=NewDi*2/Ti;
式中:Ti为地震各层的反射时间,NewVi为校正后新的层速度;
7)按坐标得到每一地震层校正前后的反射层速度、平均速度和底界深度,获得经过三维密度变速地震反射空校的结果。
2.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤1)所述的重力三维剩余密度数据体为坐标X、坐标Y、坐标Z、剩余密度Δσ,坐标Z为地下半空间某一剖分单元的中心海拔,剩余密度Δσ为该剖分单元的相对正负值,单位为g/cm3。
3.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤1)所述的地震速度数据体为坐标X、坐标Y、反射层名LName、反射层的反射时间T、反射层的层速度V、平均速度Va、反射层相对于统一基准面Height的底界深度D,T的单位为毫秒。
4.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤2)按照下式计算各层的厚度Thicki:
Thicki=Ti*Vi/2000;
按照下式计算各层的厚度相加得到各层的底界深度Di:
Di=(Thick1+Thick2+…+Thicki);
按照下式计算各层的海拔埋深Depi:
Depi=Height-Di;
式中:地震各层的反射时间为Ti,层速度为Vi,统一基准面为Height。
5.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤3)所述的厚度为剖分单元的厚度GThick,底界深度为剖分单元的底界海拔GD,按照下式计算:
GThick1=Z1-Z2,GD1=Z1-(Z1-Z2)/2;GThicki=2*(GDi-1-Zi),GDi=Zi-GThicki/2
式中Z1,Z2,……,Zi,分别为三维剩余密度数据体中地下第一,第二,……,第i层剖分单元的中心海拔。
6.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤4)所述的插值为一元三次自然样条插值:先沿x方向以纵向坐标Y为变量对三维剩余密度Δσ进行插值得到Δσx,再沿y方向以横坐标X为变量对Δσx进行插值,得到地震速度数据体所在位置的三维剩余密度Δσxy。
7.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤5)密度校正是首先由下式计算出层密度σ0:σ0=0.31*Vi 1/4;式中:Vi为层速度;再用层密度σ0加剩余密度Δσ,得到新的层密度σ1;将σ1再代入σ0=0.31*Vi 1/4,得到校正后的新层速度NewVi;新层速度NewVi=(σ1/0.31)4。
8.权利要求1的三维密度变速地震反射空校方法,其特征在于步骤7)所述的按坐标得到每一地震层校正前后的反射层速度、平均速度和底界深度,地震速度数据体格式为坐标X,坐标Y,反射层名LName,反射层的反射时间T、校正前反射层的层速度V、平均速度Va、底界埋深D,校正后反射层的层速度NewV、平均速度NewVa、底界埋深NewD。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010169276.4A CN102236104B (zh) | 2010-05-06 | 2010-05-06 | 三维密度变速地震反射空校方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201010169276.4A CN102236104B (zh) | 2010-05-06 | 2010-05-06 | 三维密度变速地震反射空校方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102236104A CN102236104A (zh) | 2011-11-09 |
CN102236104B true CN102236104B (zh) | 2014-12-24 |
Family
ID=44886937
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201010169276.4A Active CN102236104B (zh) | 2010-05-06 | 2010-05-06 | 三维密度变速地震反射空校方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102236104B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104678440A (zh) * | 2015-02-15 | 2015-06-03 | 山东科技大学 | 井约束二维地震变速度场非线性误差校正方法 |
CN106569260B (zh) * | 2015-10-10 | 2019-01-08 | 中国石油化工股份有限公司 | 针对高陡构造地区二维地震变速空校成图方法 |
CN107765341A (zh) * | 2017-10-20 | 2018-03-06 | 中国地质大学(北京) | 一种确定地层剩余密度的方法 |
CN112485841A (zh) * | 2019-09-11 | 2021-03-12 | 中国石油天然气集团有限公司 | 深层地层岩性识别方法及装置 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
CN1869734A (zh) * | 2005-05-25 | 2006-11-29 | 中国石油集团东方地球物理勘探有限责任公司 | 重力勘探数据处理变密度地形校正方法 |
CN101382598A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种真三维地震数据线性噪音的压制方法 |
CN101393270A (zh) * | 2007-09-21 | 2009-03-25 | 中国石油集团东方地球物理勘探有限责任公司 | 一种连续介质地区的表层高精度静校正方法 |
CN101625421A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 三维重磁复式采集方法 |
CN101625419A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 利用井口时间提高静校正量精度的方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6868037B2 (en) * | 2002-08-20 | 2005-03-15 | Saudi Arabian Oil Company | Use of drill bit energy for tomographic modeling of near surface layers |
-
2010
- 2010-05-06 CN CN201010169276.4A patent/CN102236104B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
CN1869734A (zh) * | 2005-05-25 | 2006-11-29 | 中国石油集团东方地球物理勘探有限责任公司 | 重力勘探数据处理变密度地形校正方法 |
CN101382598A (zh) * | 2007-09-03 | 2009-03-11 | 中国石油天然气集团公司 | 一种真三维地震数据线性噪音的压制方法 |
CN101393270A (zh) * | 2007-09-21 | 2009-03-25 | 中国石油集团东方地球物理勘探有限责任公司 | 一种连续介质地区的表层高精度静校正方法 |
CN101625421A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 三维重磁复式采集方法 |
CN101625419A (zh) * | 2008-07-08 | 2010-01-13 | 中国石油集团东方地球物理勘探有限责任公司 | 利用井口时间提高静校正量精度的方法 |
CN101661115A (zh) * | 2008-08-29 | 2010-03-03 | 中国石油集团东方地球物理勘探有限责任公司 | 基于标准格架的快速三维重力、磁力物性反演的方法 |
Non-Patent Citations (2)
Title |
---|
Property Identification of Anomalous Seismic Bodies by GMES Techniques,A Case History;Liu Yunxiang,etc.;《Applied Geophysics》;20050930;第2卷(第3期);全文 * |
自动变速空校的原理及应用;牟中海等;《石油地球物理勘探》;20000630;第35卷(第3期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN102236104A (zh) | 2011-11-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Grasmueck | 3-D ground-penetrating radar applied to fracture imaging in gneiss | |
CN101438176B (zh) | 电磁数据的时滞分析 | |
CN102419452B (zh) | 电性源瞬变电磁法全场区探测方法 | |
Balkaya et al. | Exploration for a cave by magnetic and electrical resistivity surveys: Ayvacık Sinkhole example, Bozdağ, İzmir (western Turkey) | |
CN104656157B (zh) | 一种识别页岩气甜点区的方法及装置 | |
CN108873103A (zh) | 一种结构约束的二维重力梯度和大地电磁联合反演方法 | |
CN104133245A (zh) | 一种地震资料的静校正方法及系统 | |
CN103869368A (zh) | 一种无表层调查资料约束的大炮初至综合建模静校正方法 | |
CN104537714A (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
Araffa | Delineation of groundwater aquifer and subsurface structures on North Cairo, Egypt, using integrated interpretation of magnetic, gravity, geoelectrical and geochemical data | |
Galloway et al. | Analysis and simulation of regional subsidence accompanying groundwater abstraction and compaction of susceptible aquifer systems in the USA | |
CN103149588B (zh) | 一种利用井震标定计算vti各向异性参数的方法及系统 | |
CN104597494A (zh) | 地震地层体分析方法及装置 | |
CN102236104B (zh) | 三维密度变速地震反射空校方法 | |
Li et al. | 3D electrical resistivity inversion using prior spatial shape constraints | |
Papadopoulos et al. | Electrical resistivity tomography for the modelling of cultural deposits and geomophological landscapes at Neolithic sites: a case study from Southeastern Hungary | |
CN102866422B (zh) | 一种深度域地质实体模型生成方法 | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
Šumanovac et al. | Exploration of buried carbonate aquifers by the inverse and forward modelling of the Controlled Source Audio-Magnetotelluric data | |
Tao et al. | Evaluating the joint use of GPR and ERT on mapping shallow subsurface features of karst critical zone in southwest China | |
Tang et al. | Topographic effects on long offset transient electromagnetic response | |
Peng et al. | Joint tomography of multi-cross-hole and borehole-to-surface seismic data for karst detection | |
Lin et al. | Non-invasive characterization of water-bearing strata using a combination of geophysical techniques | |
Nie et al. | Integrated ERT, seismic, and electrical resistivity imaging for geological prospecting on Metro Line R3 in Qingdao, China | |
Zhang et al. | Application of cross-borehole integrated geophysical methods for the detailed investigation of karst in urban metro construction |
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 |