CN104536046B - 基于地震记录的震源激发信号一致性评价方法 - Google Patents
基于地震记录的震源激发信号一致性评价方法 Download PDFInfo
- Publication number
- CN104536046B CN104536046B CN201410827020.6A CN201410827020A CN104536046B CN 104536046 B CN104536046 B CN 104536046B CN 201410827020 A CN201410827020 A CN 201410827020A CN 104536046 B CN104536046 B CN 104536046B
- Authority
- CN
- China
- Prior art keywords
- focus
- excitation signal
- earthquake
- epicenter excitation
- conformance
- 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
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于地震记录的震源激发信号一致性评价方法,首先提取震源地震记录中近偏移距数据中的直达波,确定不同道数据直达波信号的到时差,然后进行时差校正,将多道直达波信号进行叠加,再通过分析不同记录中的直达波信号的相似系数、相对时延以及能量差异,从而实现对不同震源的激发效果的一致性评价。解决了震源激发信号一致性的检测问题。避免了因各震源点处震源激发信号不一致而引起的地下构造探测结果较大误差的出现。不仅适用于可控震源也适用于炸药震源的一致性评价,不仅能用于单可控震源、组合震源,还可用于相控震源现场质量监控及震源维护,为地震数据处理过程对震源进行一致补偿提供必要的参考信息。
Description
技术领域:
本发明涉及一种信号的一致性评价方法,尤其是基于地震记录的震源激发信号一致性评价方法。
背景技术:
在地震勘探领域,地震数据处理、相控震源和组合震源工作方式,都需要对震源进行一致性评价。例如激发信号一致性差的震源,地震数据在地震数据处理中需要进行震源一致性补偿;激发信号一致性差的震源不能进行相控震源和组合震源方式地震勘探。
为了提高地震勘探深度,地震勘探经常使用组合震源方式,包括炸药震源组合和可控震源组合,而组合震源要求各震源激发信号严格一致。组合震源工作方式为多个震源同时工作,新出现的相控震源工作方式是多个震源协同工作,各震源采用依次延时激发的控制信号,相控震源方式对各震源激发信号的一致性要求更加严格。当前的可控震源地震勘探施工过程,要求对震源质量进行现场质量监控,震源性能达不到质量要求不能施工,这里质量监控的含义即保证震源每一次激发信号一致,并且与震源产品出厂时信号一致。可控震源使用一段时间,设备性能会因零部件老化或使用不当等原因质量下降,因此需要有效的检测方法检测震源激发信号是否与出厂一致,如果不一致要进行可控震源的维护工作。可控震源地震数据处理过程中,一旦发现个别震源点激发信号不一致,必须进行针对性地震源一致性补偿处理,这一过程也要以震源激发信号一致性作为信号补偿的依据。
CN102798889A公布了一种相控震源一致性判定方法,给出了针对相控震源、组合震源方式工作的可控震源激发信号一致性的评价方法,该方法不适用于对炸药震源以及单个可控震源的一致性评价。
发明内容:
本发明的目的就在于针对上述现有技术的不足,提供一种基于地震记录的震源激发信号一致性评价方法。
本发明的主要思想是:
在野外采集的震源地震记录中,近偏移距直达波信号信噪比高,受环境影响小,能够反映震源的激发效果。本方法首先提取震源地震记录中近偏移距数据中的直达波,利用相关检测原理,确定不同道数据直达波信号的到时差,然后进行时差校正,将多道直达波信号进行叠加,得到信噪比较高的直达波信号,再通过分析不同记录中的直达波信号的相似系数、相对时延以及能量差异,从而实现对不同震源的激发效果的一致性评价。
本发明是通过以下技术方案实现的:
基于地震记录的震源激发信号一致性评价方法,包括以下步骤:
a、输入采集系统一致的两个地震记录,记为U1(t,x)、U2(t,x),其中,t为时间,x为正整数,表示地震记录中各道数据的道号。记Ui(t,x)=(ui(t,1),ui(t,2),…,ui(t,n)),i=1,2,ui(t,n)表示第i个地震记录第n道数据;
b、在两炮地震记录中分别截取包含清晰直达波信号的近偏移距道地震记录对应的时间范围为0≤t≤t0,对应的地震道范围为xmin≤x≤xmax,即
c、记设d为中最小偏移道道号,xmin≤d≤xmax,令
d、以Si(t)为参考道,对的所有道数据作互相关运算,公式为k=xmin,xmin+1,…,xmax;
e、记Ri,k最大值对应时刻τi,k,将时移-τi,k得到k=xmin,xmin+1,…,xmax,输出对齐直达波同相轴后的记录
f、叠加的所有道数据,计算得到对应的直达波信号
g、对p1(t)和p2(t)作互相关计算,记Rp的最大值对应的时刻为tp,则两地震记录中震源相对时延Δt=tp;
h、求取p1(t)和p2(t)的相似系数,公式为<·>表示内积运算;
i、定义p1(t)和p2(t)的能量差异系数为ω=(W2-W1)/W1*100%,其中,j表示在一定采样率的情况下得到的离散时间序列;
j、根据地震勘探中对震源的实际要求,确定基于地震记录的震源激发信号一致性判定准则为:|Δt|<0.14ms、γ≥0.9、|ω|≤5%,若以上三个条件同时满足,则判定两个地震记录的震源激发信号一致性强,反之,则两个地震记录的震源激发信号一致性较弱。
有益效果:本发明公开的基于地震记录的震源激发信号一致性评价方法,解决了震源激发信号一致性的检测问题。避免了因各震源点处震源激发信号不一致而引起的地下构造探测结果较大误差的出现。经试验,在不增加施工成本的情况下,不仅适用于可控震源也适用于炸药震源的一致性评价,不仅能用于单可控震源、组合震源,还可用于相控震源现场质量监控,及震源维护,为地震数据处理过程对震源进行一致补偿提供必要的参考信息。
附图说明:
图1两炮地震记录(a)U1(t,x),(b)U2(t,x)
图2提取的直达波信号(a)p1(t),(b)p2(t)
具体实施方式:
下面结合附图和实施例做进一步的详细说明:
在本实施例中以t0=0.16s,xmax=260,xmin=240,采集系统完全一致,道间距为1m的2个地震记录为例进行震源激发信号一致性评价,但震源激发信号一致性评价方法不受实例中给出的参数的限制。
基于地震记录的震源激发信号一致性评价方法,包括以下顺序和步骤:
a、输入采集系统一致两个地震记录,记为U1(t,x)、U2(t,x),记Ui(t,x)=(ui(t,1),ui(t,2),…,ui(t,n)),i=1,2,ui(t,n)表示第i个地震记录第n道数据;
b、在两炮地震记录中分别截取包含清晰直达波信号的近偏移距道地震记录对应的时间范围为0≤t≤0.16s,对应的地震道范围为240≤x≤260,即
c、记两地震记录最小偏移道均为第250道,令
d、以Si(t)为参考道,对的所有道数据作互相关运算,公式为k=240,241,...,260;
e、记Ri,k最大值对应时刻τi,k,将时移-τi,k得到k=240,241,…,260,输出对齐直达波同相轴后的记录
f、叠加的所有道数据,计算得到对应的直达波信号其中i=1,2;
g、对p1(t)和p2(t)作互相关计算,记Rp的最大值对应的时刻tp=0s,则两地震记录中震源相对时延为Δt=tp=0s;
h、p1(t)和p2(t)的相似系数<·>表示内积运算,计算得γ=0.999;
i、定义p1(t)和p2(t)的能量差异系数为ω=(W2-W1)/W1*100%,其中,j表示在一定采样率的情况下得到的离散时间序列,采样间隔为0.2ms,按照公式求取ω=-0.1%;
j、根据震源激发信号一致性判定准则为:|Δt|=0<0.14ms、γ=0.999≥0.9、|ω|=0.1%≤5%,判定两个地震记录的震源激发信号一致性强。
Claims (1)
1.一种基于地震记录的震源激发信号一致性评价的方法,包括以下步骤:
a、输入采集系统一致的两个地震记录,记为U1(t,x)、U2(t,x),其中,t为时间,x为正整数,表示地震记录中各道数据的道号,记Ui(t,x)=(ui(t,1),ui(t,2),…,ui(t,n)),i=1,2,ui(t,n)表示第i个地震记录第n道数据;
b、在两炮地震记录中分别截取包含清晰直达波信号的近偏移距道地震记录对应的时间范围为0≤t≤t0,对应的地震道范围为xmin≤x≤xmax,即
c、记设d为中最小偏移道道号,令
d、以Si(t)为参考道,对的所有道数据作互相关运算,公式为
e、记Ri,k最大值对应时刻为τi,k,将时移-τi,k得到 输出对齐直达波同相轴后的记录
f、叠加的所有道数据,计算得到对应的直达波信号
g、对p1(t)和p2(t)作互相关计算,记Rp的最大值对应的时刻为tp,则两地震记录中震源相对时延Δt=tp;
h、求取p1(t)和p2(t)的相似系数,公式为<·>表示内积运算;
i、定义p1(t)和p2(t)的能量差异系数为ω=(W2-W1)/W1*100%,其中, j表示在一定采样率的情况下得到的离散时间序列;
j、根据地震勘探中对震源的实际要求,确定基于地震记录的震源激发信号一致性判定准则为:|Δt|<0.14ms、γ≥0.9、|ω|≤5%,若以上三个条件同时满足,则判定两个地震记录的震源激发信号一致性强,反之,则两个地震记录的震源激发信号一致性较弱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827020.6A CN104536046B (zh) | 2014-12-26 | 2014-12-26 | 基于地震记录的震源激发信号一致性评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410827020.6A CN104536046B (zh) | 2014-12-26 | 2014-12-26 | 基于地震记录的震源激发信号一致性评价方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104536046A CN104536046A (zh) | 2015-04-22 |
CN104536046B true CN104536046B (zh) | 2017-03-29 |
Family
ID=52851608
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410827020.6A Expired - Fee Related CN104536046B (zh) | 2014-12-26 | 2014-12-26 | 基于地震记录的震源激发信号一致性评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104536046B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107462926B (zh) * | 2016-06-03 | 2019-04-12 | 中国石油化工股份有限公司 | 一种野外试验最佳激发井深的选择方法 |
CN106646603B (zh) * | 2017-01-04 | 2018-11-27 | 中海石油(中国)有限公司 | 一种实际时移地震资料处理差异的可靠性判断方法 |
CN112180469B (zh) * | 2019-07-03 | 2023-10-27 | 中国石油天然气集团有限公司 | 震源有线一致性测试方法 |
CN113109865A (zh) * | 2020-01-09 | 2021-07-13 | 中国石油天然气集团有限公司 | 混源激发地震波的处理方法及装置 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320095A (zh) * | 2008-07-22 | 2008-12-10 | 吉林大学 | 单震源定向照明地震信号合成方法 |
CN101980054A (zh) * | 2010-10-22 | 2011-02-23 | 中国石油化工股份有限公司 | 一种在高密度地震静校正处理中建立近地表速度模型的方法 |
CN102798889A (zh) * | 2012-04-26 | 2012-11-28 | 吉林大学 | 相控震源一致性判定方法 |
CN103344989A (zh) * | 2013-07-11 | 2013-10-09 | 吉林大学 | 可控震源地震记录中脉冲噪声干扰的分析方法 |
-
2014
- 2014-12-26 CN CN201410827020.6A patent/CN104536046B/zh not_active Expired - Fee Related
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101320095A (zh) * | 2008-07-22 | 2008-12-10 | 吉林大学 | 单震源定向照明地震信号合成方法 |
CN101980054A (zh) * | 2010-10-22 | 2011-02-23 | 中国石油化工股份有限公司 | 一种在高密度地震静校正处理中建立近地表速度模型的方法 |
CN102798889A (zh) * | 2012-04-26 | 2012-11-28 | 吉林大学 | 相控震源一致性判定方法 |
CN103344989A (zh) * | 2013-07-11 | 2013-10-09 | 吉林大学 | 可控震源地震记录中脉冲噪声干扰的分析方法 |
Non-Patent Citations (2)
Title |
---|
地震勘探中相控阵震源的方向特性研究;王忠仁等;《地球物理学报》;20060730(第04期);全文 * |
相控震源定向地震波信号分析;姜弢等;《地球物理学报》;20080915(第05期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN104536046A (zh) | 2015-04-22 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US11016214B2 (en) | Dolomite reservoir prediction method and system based on well and seismic combination, and storage medium | |
CN106353792A (zh) | 一种适用于水力压裂微震震源定位的方法 | |
CN101852863B (zh) | 一种利用高精度单道频谱分析技术处理地震数据的方法 | |
Zeng et al. | Surface microseismic monitoring of hydraulic fracturing of a shale‐gas reservoir using short‐period and broadband seismic sensors | |
CN104536046B (zh) | 基于地震记录的震源激发信号一致性评价方法 | |
CN105697002A (zh) | 一种用于识别煤系地层岩性的方法 | |
CN104345345A (zh) | 一种页岩储层总有机碳toc含量预测方法 | |
CN103675907A (zh) | 一种基于岩相约束的avo反演油气检测方法 | |
CN104749624A (zh) | 一种地震岩相识别及其不确定性定量评价同步实现方法 | |
CN110133715A (zh) | 一种基于初至时差和波形叠加的微地震震源定位方法 | |
Sabbione et al. | A robust method for microseismic event detection based on automatic phase pickers | |
CN107678064B (zh) | 一种声波时差实时提取方法 | |
CN104865614A (zh) | 基于可变骨架参数的复杂储层流体识别方法 | |
Sharma et al. | Ground‐motion prediction equations for the Geysers geothermal area based on induced seismicity records | |
Song et al. | Full‐waveform based microseismic event detection and signal enhancement: An application of the subspace approach | |
Maxwell et al. | Tracking microseismic signals from the reservoir to surface | |
CN105093297A (zh) | 一种微地震定位精度的质量控制方法 | |
Iaccarino et al. | Onsite earthquake early warning: predictive models for acceleration response spectra considering site effects | |
Reiter et al. | Source characteristics of seismicity associated with underground wastewater disposal: A case study from the 2008 Dallas-Fort Worth earthquake sequence | |
CN101852864A (zh) | 一种利用地表一致性统计频谱分析技术处理海量地震数据的方法 | |
CN107870359A (zh) | 微地震事件识别方法及装置 | |
Panebianco et al. | Semi-automated template matching and machine-learning based analysis of the August 2020 Castelsaraceno microearthquake sequence (southern Italy) | |
CN103091728B (zh) | 水力压裂次生声波模拟方法和压裂裂缝监测技术评价方法 | |
CN115877464A (zh) | 一种岩性识别方法、装置、计算机设备及存储介质 | |
CN107679614B (zh) | 一种基于粒子群优化的声波时差实时提取方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
C10 | Entry into 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 |
Granted publication date: 20170329 Termination date: 20181226 |
|
CF01 | Termination of patent right due to non-payment of annual fee |