CN111983685B - τ-p域地表非一致性静校正方法 - Google Patents
τ-p域地表非一致性静校正方法 Download PDFInfo
- Publication number
- CN111983685B CN111983685B CN202010701134.1A CN202010701134A CN111983685B CN 111983685 B CN111983685 B CN 111983685B CN 202010701134 A CN202010701134 A CN 202010701134A CN 111983685 B CN111983685 B CN 111983685B
- Authority
- CN
- China
- Prior art keywords
- seismic data
- static correction
- domain
- tau
- time
- 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.)
- Expired - Fee Related
Links
- 238000012937 correction Methods 0.000 title claims abstract description 62
- 230000003068 static effect Effects 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 22
- 230000001131 transforming effect Effects 0.000 claims abstract description 8
- 238000007781 pre-processing Methods 0.000 claims abstract description 4
- 230000009466 transformation Effects 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 239000006185 dispersion Substances 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000012545 processing Methods 0.000 abstract description 3
- 230000007547 defect Effects 0.000 abstract description 2
- 239000010410 layer Substances 0.000 description 16
- 238000004088 simulation Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000012360 testing method Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000005422 blasting Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000002344 surface layer Substances 0.000 description 1
Images
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/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- 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/32—Transforming one recording into another or one representation into another
- G01V1/325—Transforming one representation into another
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
本发明涉及一种τ‑p域地表非一致性长波长静校正方法,属于地球物理勘探的地震资料处理领域。本发明主要包括如下步骤:1)对陆上地震资料进行预处理;2)建立近地表速度模型;3)将预处理后的地震资料从时间‑空间域变换到τ‑p域;4)计算每一个射线参数对应的静校正量;5)对变换后的地震资料进行静校正;6)将静校正后的地震资料从τ‑p域变换到时间‑空间域。应用本发明,可以实现对不同深度界面的反射波沿其在近地表低速层中的实际传播路径进行校正,克服了常用地表一致性静校正方法的不足,得到更加准确的静校正结果。
Description
技术领域
本发明涉及一种τ-p域地表非一致性长波长静校正方法,属于地球物理勘探的地震资料处理领域。
背景技术
目前在陆上地震资料处理过程中广泛采用基于地表一致性假设的静校正方法。地表一致性假设认为近地表低速层的速度远小于其下覆基岩速度,反射波在近地表层中沿垂直方向传播。因此,地表同一位置接收的不同深度界面的反射波在近地表低速层中的传播路径相同,从而传播时间相同。将反射波在近地表低速层中的传播时间消去,就避免了近地表低速层对深部反射波同相轴的影响。这类方法在近地表低速层速度和基岩速度差异大的地区通常可以取得令人满意的结果。
然而,在近地表低速层速度和基岩速度差异小的地区,反射波在近地表低速层中沿倾斜方向传播。反射波在近地表低速层中的传播方向受反射界面的深度、炮检点间距、近地表低速层速度和基岩速度的影响。地表同一位置接收的不同深度界面的反射波在近地表低速层中的传播路径差异较大,具有不同的传播时间。因此,针对不同的反射波需要不同的静校正量。将基于地表一致性假设的静校正方法应用于这些地区采集的地震数据,难以消除近地表低速层对反射波同相轴的影响,降低了后续的成像质量。
发明内容
针对现有技术存在的上述缺陷,本发明提出了一种τ-p域地表非一致性长波长静校正方法,在τ-p域(截距时间-射线参数域)对不同的反射波沿其在近地表低速层中实际的传播路径进行校正,解决了基于地表一致性假设的静校正方法沿垂直方向进行校正的问题。
本发明是采用以下的技术方案实现的:本发明所述一种τ-p域地表非一致性长波长静校正方法,包括如下步骤:
步骤一:对陆上地震资料进行预处理;预处理包括振幅均衡、去除面波、压制随机噪声;
步骤二:利用初至走时数据建立近地表速度模型,近地表速度模型采用正方形网格离散;
步骤三:将预处理后的地震资料从时间-空间域变换到τ-p域;采用下式进行变换:
式中,m是变换后的地震资料,τ是截距时间,p是射线参数,i是道号,n是最大的道号,d是预处理后的地震资料,xi是第i道地震资料对应的偏移距,t是采样时间;
步骤四:计算射线参数对应的静校正量;采用下式计算静校正量:
式中,Δτ是静校正量,k是近地表速度模型离散网格序号,nz是离散网格的最大序号,dzk是第k个网格在垂向上的大小,qk是射线穿过第k个网格时的垂直慢度;
步骤五:对变换后的地震资料进行静校正;采用下式进行静校正:
m′(τ′,p)=m(τ+Δτ,p), (3)
式中,m’是静校正后的地震资料,τ’是静校正后的截距时间,p是射线参数,m是变换后的地震资料,τ是截距时间,Δτ是静校正量;
步骤六:将静校正后的地震资料从τ-p域变换到时间-空间域;采用下式进行变换:
式中,d’是地表非一致性静校正后的地震资料,x是偏移距,t是采样时间,j是射线参数的序号,np是射线参数的最大序号,m’是静校正后的地震资料,τ’是静校正后的截距时间,pj是第j个射线参数。
本发明的有益效果是:将预处理后的地震资料从时间-空间域变换到τ-p域,然后在τ-p域对每一个射线参数,即p值对应的地震资料进行截距时间校正,可以实现对近地表低速层中路径相同的反射波进行相同的校正,最后将地震资料从τ-p域变换到时间-空间域,实现了根据射线路径进行校正而不是根据地表位置对应的垂直路径进行校正。本发明计算简单,易于实现,静校正效果好。
附图说明
图1为本发明的流程图;
图2为真实速度模型图;
图3为正演模拟得到的第23炮地震记录图;
图4为τ-p域的第23炮地震记录图;
图5为第23炮地震记录在τ-p域进行炮点端静校正后的地震记录图;
图6为时间-空间域经过地表非一致性静校正后的第23炮地震记录图;
图7为时间-空间域经过地表一致性静校正后的第23炮地震记录图;
图8为去掉风化层后正演模拟得到的第23炮地震记录图;
具体实施方式
为了使本发明的目的、技术方案更加清楚,下面以理论模型模拟数据为例并结合附图,对本发明作进一步详细说明。
本发明主要包括如下步骤:
步骤一:对陆上地震资料进行预处理;预处理包括振幅均衡、去除面波、压制随机噪声;
步骤二:利用初至走时数据建立近地表速度模型,近地表速度模型采用正方形网格离散;
步骤三:将预处理后的地震资料从时间-空间域变换到τ-p域;采用下式进行变换:
式中,m是变换后的地震资料,τ是截距时间,p是射线参数,i是道号,n是最大的道号,d是预处理后的地震资料,xi是第i道地震资料对应的偏移距,t是采样时间;
步骤四:计算射线参数对应的静校正量;采用下式计算静校正量:
式中,Δτ是静校正量,k是近地表速度模型离散网格序号,nz是离散网格的最大序号,dzk是第k个网格在垂向上的大小,qk是射线穿过第k个网格时的垂直慢度;
步骤五:对变换后的地震资料进行静校正;采用下式进行静校正:
m′(τ′,p)=m(τ+Δτ,p), (3)
式中,m’是静校正后的地震资料,τ’是静校正后的截距时间,p是射线参数,m是变换后的地震资料,τ是截距时间,Δτ是静校正量;
步骤六:将静校正后的地震资料从τ-p域变换到时间-空间域;采用下式进行变换:
式中,d’是地表非一致性静校正后的地震资料,x是偏移距,t是采样时间,j是射线参数的序号,np是射线参数的最大序号,m’是静校正后的地震资料,τ’是静校正后的截距时间,pj是第j个射线参数。
实施例一:
下面结合具体实施方式,对于本发明的理论模型测试进行解释和说明。
为了进一步说明本方法的实现思路及实现过程并证明方法的有效性,用水平层状模型进行测试,并和地表一致性静校正方法结果进行比较。
S1:建立一个真实速度模型,如图2所示。真实速度模型宽度为20km,深度为2km。真实速度模型共有4层。Z方向0.35km至0.8km为近地表低速层,其下方有3个匀速地层。采用正方形网格离散真实速度模型,网格大小为10×10m。
S2:观测系统:采用中间放炮两边接收的观测方式。炮点以100m的间隔均匀分布在6.1-13.8km。每炮有241个检波点接收地震记录。检波点以50m的间隔均匀分布在炮点两侧。最小偏移距为0km,最大偏移距为6km。地震资料的采样时间为3.5s,时间采样间隔为1ms。
S3:由真实速度模型(详见图2)和震源函数为主频15Hz的雷克子波,利用二阶声波方程正演模拟得到地震资料。图3是正演模拟得到的第23炮地震记录。
S4:对正演模拟得到的地震记录进行振幅均衡预处理。
S5:利用公式1,将模拟得到的地震记录变换到τ-p域。图4是τ-p域的第23炮地震记录。
S6:利用公式2,计算每一个射线参数对应的静校正量。
S7:利用公式3,对τ-p域的地震记录进行静校正。图5是第23炮地震记录在τ-p域进行炮点端静校正后的地震记录。和图4相比,图5中同相轴的截距时间变小了。
S8:利用公式4,将静校正后的地震记录变换到时间-空间域,从而实现地表非一致性静校正。图6是时间-空间域经过地表非一致性静校正后的第23炮地震记录。
图7是时间-空间域经过地表一致性静校正后的第23炮地震记录。对比图6和图7,可以发现图6中只有两个双曲线形状的反射同相轴,而图7中有三个双曲线形状的反射同相轴。图8是去掉地表低速层后正演模拟得到的第23炮地震记录。图8中只有两个双曲线形状的反射同相轴。对比图6、图7和图8,可以发现图6和图8在同相轴形态上更加接近,表明经过地表非一致性静校正后的结果更加准确。
以上所述仅为本发明的较佳实施例而己,并不以本发明为限制,凡在本发明的精神和原则之内所作的均等修改、等同替换和改进等,均应包含在本发明的专利涵盖范围内。
Claims (1)
1.τ-p域地表非一致性长波长静校正方法,其特征在于,包括如下步骤:
步骤一:对陆上地震资料进行预处理;预处理包括振幅均衡、去除面波、压制随机噪声;
步骤二:利用初至走时数据建立近地表速度模型,近地表速度模型采用正方形网格离散;
步骤三:将预处理后的地震资料从时间-空间域变换到τ-p域;采用下式进行变换:
式中,m是变换后的地震资料,τ是截距时间,p是射线参数,i是道号,n是最大的道号,d是预处理后的地震资料,xi是第i道地震资料对应的偏移距,t是采样时间;
步骤四:计算射线参数对应的静校正量;采用下式计算静校正量:
式中,Δτ是静校正量,k是近地表速度模型离散网格序号,nz是离散网格的最大序号,dzk是第k个网格在垂向上的大小,qk是射线穿过第k个网格时的垂直慢度;
步骤五:对变换后的地震资料进行静校正;采用下式进行静校正:
式中,m’是静校正后的地震资料,τ’是静校正后的截距时间,p是射线参数,m是变换后的地震资料,τ是截距时间,Δτ是静校正量;
步骤六:将静校正后的地震资料从τ-p域变换到时间-空间域;采用下式进行变换:
式中,d’是地表非一致性静校正后的地震资料,x是偏移距,t是采样时间,j是射线参数的序号,np是射线参数的最大序号,m’是静校正后的地震资料,τ’是静校正后的截距时间,pj是第j个射线参数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010701134.1A CN111983685B (zh) | 2020-07-21 | 2020-07-21 | τ-p域地表非一致性静校正方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010701134.1A CN111983685B (zh) | 2020-07-21 | 2020-07-21 | τ-p域地表非一致性静校正方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN111983685A CN111983685A (zh) | 2020-11-24 |
CN111983685B true CN111983685B (zh) | 2021-11-12 |
Family
ID=73439230
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010701134.1A Expired - Fee Related CN111983685B (zh) | 2020-07-21 | 2020-07-21 | τ-p域地表非一致性静校正方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111983685B (zh) |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4839869A (en) * | 1986-10-06 | 1989-06-13 | Shell Oil Company | Methods for processing converted wave seismic data |
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
WO2015160652A1 (en) * | 2014-04-17 | 2015-10-22 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data |
CN107229073A (zh) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | 地震数据处理方法和装置 |
CN110858003A (zh) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | 基于等效速度的非地表一致性静校正方法 |
CN110879413A (zh) * | 2018-09-05 | 2020-03-13 | 中国石油化工股份有限公司 | 一种射线参数域转换波静校正方法及系统 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP1849026A2 (en) * | 2005-02-18 | 2007-10-31 | BP Corporation North America Inc. | System and method for using time-distance characteristics in acquisition, processing and imaging of t-csfm data |
US7397728B2 (en) * | 2005-08-26 | 2008-07-08 | Westerngeco L.L.C. | Method for processing a record of seismic traces |
CN105093326B (zh) * | 2014-05-23 | 2017-09-19 | 中国石油化工股份有限公司 | 一种地震勘探资料中波长剩余静校正处理方法 |
-
2020
- 2020-07-21 CN CN202010701134.1A patent/CN111983685B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4839869A (en) * | 1986-10-06 | 1989-06-13 | Shell Oil Company | Methods for processing converted wave seismic data |
US5157638A (en) * | 1992-01-13 | 1992-10-20 | Conoco Inc. | Method of deriving statics corrections from common reflection point gathers |
WO2015160652A1 (en) * | 2014-04-17 | 2015-10-22 | Saudi Arabian Oil Company | Generating subterranean imaging data based on vertical seismic profile data |
CN107229073A (zh) * | 2016-03-24 | 2017-10-03 | 中国石油化工股份有限公司 | 地震数据处理方法和装置 |
CN110858003A (zh) * | 2018-08-22 | 2020-03-03 | 中国石油化工股份有限公司 | 基于等效速度的非地表一致性静校正方法 |
CN110879413A (zh) * | 2018-09-05 | 2020-03-13 | 中国石油化工股份有限公司 | 一种射线参数域转换波静校正方法及系统 |
Non-Patent Citations (1)
Title |
---|
"长波长静校正问题的一种解决方法";戴云 等;《石油地球物理勘探》;20000630;第35卷(第3期);第315-325页 * |
Also Published As
Publication number | Publication date |
---|---|
CN111983685A (zh) | 2020-11-24 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102890290B (zh) | 一种起伏地表条件下的叠前深度偏移方法 | |
Yue et al. | Prestack Gaussian beam depth migration under complex surface conditions | |
CN109669212B (zh) | 地震数据处理方法、地层品质因子估算方法与装置 | |
CN108196305B (zh) | 一种山地静校正方法 | |
CN101545986A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
CN104483704B (zh) | 基于avo异常类型约束的剩余相位校正方法 | |
CN105182408A (zh) | 一种合成地震记录的制作方法和装置 | |
CN109765616B (zh) | 一种保幅波场延拓校正方法及系统 | |
CN101021568A (zh) | 基于最大能量旅行时计算的三维积分叠前深度偏移方法 | |
CN116520419B (zh) | 一种热流体裂缝通道识别方法 | |
CN110389377B (zh) | 基于波形互相关系数相乘的微震偏移成像定位方法 | |
CN111983685B (zh) | τ-p域地表非一致性静校正方法 | |
CN111352153B (zh) | 一种基于瞬时相位互相关加权的微地震干涉定位方法 | |
CN111257938A (zh) | 基于小波互相关时移地震虚拟震源波场重构方法和系统 | |
CN112269211A (zh) | τ-p域地表非一致性长波长静校正方法 | |
CN113447981B (zh) | 一种基于共成像点道集的反射全波形反演方法 | |
CN110780346A (zh) | 一种隧道超前探测复杂地震波场的分离方法 | |
CN111929731B (zh) | 地表一致性和非一致性联合静校正方法 | |
CN113406700B (zh) | 一种地表主动源反射波干涉成像方法 | |
CN112099090B (zh) | 地震资料视速度域非一致性长波长静校正方法 | |
CN112269210A (zh) | 地表一致性和非一致性联合静校正方法 | |
CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
CN112379431B (zh) | 复杂地表条件下ps波地震数据偏移成像方法及系统 | |
CN116125535B (zh) | 三维vsp成像的方法及装置 | |
CN112379430B (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 | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20211112 |