CN109655901B - 一种频率域自适应偏振角计算方法及*** - Google Patents

一种频率域自适应偏振角计算方法及*** Download PDF

Info

Publication number
CN109655901B
CN109655901B CN201710942413.5A CN201710942413A CN109655901B CN 109655901 B CN109655901 B CN 109655901B CN 201710942413 A CN201710942413 A CN 201710942413A CN 109655901 B CN109655901 B CN 109655901B
Authority
CN
China
Prior art keywords
frequency
components
frequency domain
polarization angle
component
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
CN201710942413.5A
Other languages
English (en)
Other versions
CN109655901A (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 Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
Original Assignee
China Petroleum and Chemical Corp
Sinopec Geophysical Research Institute
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 China Petroleum and Chemical Corp, Sinopec Geophysical Research Institute filed Critical China Petroleum and Chemical Corp
Priority to CN201710942413.5A priority Critical patent/CN109655901B/zh
Publication of CN109655901A publication Critical patent/CN109655901A/zh
Application granted granted Critical
Publication of CN109655901B publication Critical patent/CN109655901B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/30—Analysis
    • G01V1/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity

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)
  • Measurement Of Mechanical Vibrations Or Ultrasonic Waves (AREA)

Abstract

本发明提出一种频率域自适应偏振角计算方法及***，该方法包括：将每个检波器的三个分量变换到频率域Uij(ω)；计算每个频率的自适应阈值范围；利用自适应阈值范围识别异常频率成份，确定频率集合Sj；利用频率集合Sj，通过相关目标函数计算偏振角。本发明相对于已有技术最大的创新点在于在频率域中自适应计算偏振角，这种方法可以在频率域中选择品质高的频率成份参与计算，提高偏振角计算的稳定性。

Description

一种频率域自适应偏振角计算方法及***
技术领域
本发明属于微地震监测中一种偏振分析方法。具体地，涉及在频率域中自适应选择高品质频率成份计算偏振角。
背景技术
近20年伴随着非常规低渗透油气藏的开采，水力压裂技术已成为油气田增产的一个主要措施。而微地震监测技术在水力压裂过程中起到监测评估的作用。微地震监测技术的地球物理基础是通过观测、分析水力压裂过程中所产生的微小地震事件来评估生产活动效果。井中监测方式是微地震监测方式中的一种重要手段，它是在压裂井中进行水力压裂施工，在临近观测井中布置多级检波器排列进行微地震信号监测。
要实现对微地震事件的定位，偏振角计算是定位过程中必不可少的步骤。计算偏振角有许多种方法，例如矢端曲线能量准则法、协方差法、最大能量准则法。这些方法在信噪比高的情况效果都不错。但是在微地震监测中有监测井距离压裂井较近(约300～600米)，噪音比较强，再加上微地震事件的震级较低，造成了偏振角计算不稳定。
发明内容
为了解决现有技术中存在的问题，本发明给出了一种频率域自适应偏振角计算方法。
根据本发明的一个方面，提供一种频率域自适应偏振角计算方法，该方法包括：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角。
进一步地，将每个检波器的三个分量变换到频率域中，记作Uij(ω)＝|FFT(uij(t))|，
其中，i表示三个分量，取值1，2，3，分别对应两个水平分量和垂直分量，j表示检波器的序号1,…,N。
进一步地，设纵波的偏振角为
Figure BDA0001430982170000021
则旋转完成后切向分量uT、径向分量uR和垂向分量uZ为：
Figure BDA0001430982170000022
其中uE、uN为三分量检波器的两个水平分量，uZ为三分量检波器的垂向分量。
进一步地，计算每个频率的自适应阈值范围包括：
统计每个频率成份的平均值：
Figure BDA0001430982170000023
统计每个频率成份的方差：
Figure BDA0001430982170000024
进一步地，每个频率的自适应阈值为：
Th(ω)＝E(ω)+σ(ω) (4)
进一步地，如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj；反之，则该频率成份属于集合Sj。
进一步地，在频率域中相关目标函数为：
Figure BDA0001430982170000031
式中，Sj为选定频率的集合，j表示检波器序号1,…,N。
根据本发明的另一方面，提供一种频率域自适应偏振角计算***，该***包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角。
进一步地，每个频率的自适应阈值为：
Th(ω)＝E(ω)+σ(ω) (4)
其中每个频率成份的平均值为：
Figure BDA0001430982170000032
每个频率成份的方差为：
Figure BDA0001430982170000033
进一步地，如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj；反之，则该频率成份属于集合Sj；
在频率域中相关目标函数为：
Figure BDA0001430982170000041
式中，Sj为选定频率的集合，j表示检波器序号1,…,N。
本发明相对于已有技术最大的创新点在于在频率域中自适应计算偏振角，这种方法可以在频率域中选择品质高的频率成份参与计算，提高偏振角计算的稳定性。
附图说明
通过结合附图对本公开示例性实施方式进行更详细的描述，本公开的上述以及其它目的、特征和优势将变得更加明显，其中，在本公开示例性实施方式中，相同的参考标号通常代表相同部件。
图1示出了本发明实施例的微地震井中监测示意图。
图2示出了本发明实施例的方法流程图。
具体实施方式
下面将参照附图更详细地描述本公开的优选实施方式。虽然附图中显示了本公开的优选实施方式，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施方式所限制。相反，提供这些实施方式是为了使本公开更加透彻和完整，并且能够将本公开的范围完整地传达给本领域的技术人员。
本发明是一种频率域自适应偏振角计算方法。要实现对微地震事件的定位，偏振角计算是定位过程中必不可少的步骤。计算偏振角有许多种方法，例如矢端曲线能量准则法、协方差法、最大能量准则法。这些方法在信噪比高的情况效果都不错。但是在微地震监测中有监测井距离压裂井较近(约300～600米)，噪音比较强，再加上微地震事件的震级较低，造成了偏振角计算不稳定。因此，本专利给出一种频率域自适应偏振角计算方法。
根据本发明的一个实施方式，提供一种频率域自适应偏振角计算方法，该方法包括：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角。
偏振角分析的理论依据是纵波的偏振方向与波的传播方向一致。设纵波的偏振角为
Figure BDA0001430982170000051
则旋转完成后切向分量uT、径向分量uR和垂向分量uZ为：
Figure BDA0001430982170000052
其中uE、uN为三分量检波器的两个水平分量，uZ为三分量检波器的垂向分量。
在微地震三分量中噪音主要来源于井场噪音，在频率域中表现为近单频的能量极值。因此在频率域中构建目标函数时可以有选择性地选取频率成份参与目标计算，避免噪音干扰目标函数计算。
优选地，频率集合Sj选择过程：
(1)将每个检波器的三个分量变换到频率域中，记作Uij(ω)＝|FFT(uij(t))|，i表示三个分量，取值1，2，3(即E，N，Z)，分别对应两个水平分量和垂直分量，j表示检波器的序号1,…,N。
(2)计算每个频率的自适应阈值范围
统计每个频率成份的平均值：
Figure BDA0001430982170000053
统计每个频率成份的方差：
Figure BDA0001430982170000061
则每个频率的自适应阈值：
Th(ω)＝E(ω)+σ(ω) (4)
(3)利用自适应阈值范围识别异常频率成份
如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj中；反之，则该频率成份属于集合Sj。
经过旋转后，得到切向分量uT(t)垂直于径向分量uR(t)与垂向分量uZ(t)构成平面，其能量最小。而且在检波器的三个分量中微地震纵波到达时间是一致的，因此三个分量的两两相关在零点处取极小值。
因此，在频率域中相关目标函数为：
Figure BDA0001430982170000062
式中，
Figure BDA0001430982170000063
为偏振角，Sj为选定频率的集合，j表示检波器序号，取值1,…,N。
该目标函数(5)的求取可以选用全局寻优计算方法(例如，模拟退火算法、差分进化算法)进行求解。
根据本发明的另一实施方式，提供一种频率域自适应偏振角计算***，该***包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角。
进一步地，每个频率的自适应阈值为：
Th(ω)＝E(ω)+σ(ω) (4)
其中每个频率成份的平均值为：
Figure BDA0001430982170000071
每个频率成份的方差为：
Figure BDA0001430982170000072
进一步地，如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj；反之，则该频率成份属于集合Sj；
在频率域中相关目标函数为：
Figure BDA0001430982170000073
式中，Sj为选定频率的集合，j表示检波器序号1,…,N。
为便于理解本发明实施例的方案及其效果，以下给出一个具体应用示例。本领域技术人员应理解，该示例仅为了便于理解本发明，其任何具体细节并非意在以任何方式限制本发明。
如图2所示，首先输入数据，并将记录数据变换到频率域Uij(ω)。
按照公式(2)-(4)计算平均值、方差和自适应阈值，然后将振幅值小于自适应阈值的频率纳入计算序列，最后按照公式(5)计算偏振角。
本发明的方法可以在频率域中选择品质高的频率成份参与计算，提高偏振角计算的稳定性。
以上已经描述了本公开的各实施例，上述说明是示例性的，并非穷尽性的，并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下，对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。本文中所用术语的选择，旨在最好地解释各实施例的原理、实际应用或对市场中的技术改进，或者使本技术领域的其它普通技术人员能理解本文披露的各实施例。

Claims (10)

1.一种频率域自适应偏振角计算方法，其特征在于，该方法包括：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角；
其中，i表示三个分量，取值1，2，3，分别对应两个水平分量和垂直分量，j表示检波器的序号1,…,N。
2.根据权利要求1所述的频率域自适应偏振角计算方法，其特征在于，将每个检波器的三个分量变换到频率域中，记作Uij(ω)＝|FFT(uij(t))|。
3.根据权利要求2所述的频率域自适应偏振角计算方法，其特征在于，设纵波的偏振角为
Figure FDA0002366894790000011
则旋转完成后切向分量uT、径向分量uR和垂向分量uZ为：
Figure FDA0002366894790000012
其中uE、uN为三分量检波器的两个水平分量，uZ为三分量检波器的垂向分量。
4.根据权利要求1所述的频率域自适应偏振角计算方法，其特征在于，计算每个频率的自适应阈值范围包括：
统计每个频率成份的平均值：
Figure FDA0002366894790000021
统计每个频率成份的方差：
Figure FDA0002366894790000022
5.根据权利要求4所述的频率域自适应偏振角计算方法，其特征在于，每个频率的自适应阈值为：
Th(ω)＝E(ω)+σ(ω) (4)。
6.根据权利要求5所述的频率域自适应偏振角计算方法，其特征在于，如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj；反之，则该频率成份属于集合Sj。
7.根据权利要求5所述的频率域自适应偏振角计算方法，其特征在于，在频率域中相关目标函数为：
Figure FDA0002366894790000023
式中，Sj为选定频率的集合，j表示检波器序号1,…,N，
Figure FDA0002366894790000024
为偏振角。
8.一种频率域自适应偏振角计算***，其特征在于，该***包括：
存储器，存储有计算机可执行指令；
处理器，所述处理器运行所述存储器中的计算机可执行指令，执行以下步骤：
将每个检波器的三个分量变换到频率域Uij(ω)；
计算每个频率的自适应阈值范围；
利用自适应阈值范围识别异常频率成份，确定频率集合Sj；
利用频率集合Sj，通过相关目标函数计算偏振角；
其中，i表示三个分量，取值1，2，3，分别对应两个水平分量和垂直分量，j表示检波器的序号1,…,N。
9.根据权利要求8所述的频率域自适应偏振角计算***，其特征在于，每个频率的自适应阈值为：
Th(ω)＝E(ω)+σ(ω) (4)
其中每个频率成份的平均值为：
Figure FDA0002366894790000031
每个频率成份的方差为：
Figure FDA0002366894790000032
10.根据权利要求9所述的频率域自适应偏振角计算***，其特征在于，如果Uij(ω)大于Th(ω)，则该频率成份不属于集合Sj；反之，则该频率成份属于集合Sj；
在频率域中相关目标函数为：
Figure FDA0002366894790000033
式中，Sj为选定频率的集合，j表示检波器序号1,…,N，
Figure FDA0002366894790000034
为偏振角。
CN201710942413.5A 2017-10-11 2017-10-11 一种频率域自适应偏振角计算方法及*** Active CN109655901B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710942413.5A CN109655901B (zh) 2017-10-11 2017-10-11 一种频率域自适应偏振角计算方法及***

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710942413.5A CN109655901B (zh) 2017-10-11 2017-10-11 一种频率域自适应偏振角计算方法及***

Publications (2)

Publication Number Publication Date
CN109655901A CN109655901A (zh) 2019-04-19
CN109655901B true CN109655901B (zh) 2020-08-25

Family

ID=66109089

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710942413.5A Active CN109655901B (zh) 2017-10-11 2017-10-11 一种频率域自适应偏振角计算方法及***

Country Status (1)

Country Link
CN (1) CN109655901B (zh)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112558147B (zh) * 2019-09-25 2023-10-10 中国石油化工股份有限公司 一种井中微地震数据的偏振分析方法及***
CN113658110A (zh) * 2021-07-22 2021-11-16 西南财经大学 一种基于动态领域自适应学习的医学图像识别方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009007822A3 (en) * 2007-07-06 2009-03-05 Schlumberger Technology Bv Methods and systems for processing microseismic data
CN102879817A (zh) * 2012-07-23 2013-01-16 北京多分量地震技术研究院 基于地面地震数据获取地下裂缝信息的控制方法
CN105785440A (zh) * 2016-02-29 2016-07-20 河南理工大学 一种矿井槽波双分量地震信号频散曲线提取方法
CN106249295A (zh) * 2015-06-15 2016-12-21 中国石油化工股份有限公司 一种井中微地震p、s波联合快速定位方法及***
CN107132575A (zh) * 2017-05-12 2017-09-05 中国地质大学（北京） 基于横波极化分析预测裂缝方位角的方法

Family Cites Families (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2781937B1 (en) * 2013-03-22 2022-09-28 CGG Services SAS Vector-dip filtering of seismic data in the time-frequency domain

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2009007822A3 (en) * 2007-07-06 2009-03-05 Schlumberger Technology Bv Methods and systems for processing microseismic data
CN102879817A (zh) * 2012-07-23 2013-01-16 北京多分量地震技术研究院 基于地面地震数据获取地下裂缝信息的控制方法
CN106249295A (zh) * 2015-06-15 2016-12-21 中国石油化工股份有限公司 一种井中微地震p、s波联合快速定位方法及***
CN105785440A (zh) * 2016-02-29 2016-07-20 河南理工大学 一种矿井槽波双分量地震信号频散曲线提取方法
CN107132575A (zh) * 2017-05-12 2017-09-05 中国地质大学（北京） 基于横波极化分析预测裂缝方位角的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"基于Matlab的偏振分析在微地震三维定位中的应用";夏媛媛 等;《金属矿山》;20100331;第119-121页 *

Also Published As

Publication number Publication date
CN109655901A (zh) 2019-04-19

Similar Documents

Publication Publication Date Title
Ford et al. Network sensitivity solutions for regional moment-tensor inversions
CN102116868B (zh) 一种地震波分解方法
Pioldi et al. Earthquake‐induced structural response output‐only identification by two different Operational Modal Analysis techniques
EP3210391B1 (en) Reverberation estimator
CN102610227A (zh) 声音信号处理设备、声音信号处理方法和程序
CN111399048B (zh) 一种对断溶体计算相关属性及数据加权重构的方法
CN112068195B (zh) 微地震p＆s波匹配事件初至自动拾取方法和计算机存储介质
CN109655901B (zh) 一种频率域自适应偏振角计算方法及***
Li et al. A Denoising Method for Ship‐Radiated Noise Based on Optimized Variational Mode Decomposition with Snake Optimization and Dual‐Threshold Criteria of Correlation Coefficient
CN114609669B (zh) 基于方位弹性阻抗的hti型裂缝储层参数预测方法及***
Xiao et al. Research on fault feature extraction method of rolling bearing based on NMD and wavelet threshold denoising
CN105929444B (zh) 一种基于互相关偏移与最小二乘思想的微地震定位方法
CN101825722B (zh) 一种鲁棒的地震信号瞬时频率的估计方法
CN112068194B (zh) 微地震弱事件p波初至自动拾取方法和计算机存储介质
CN107479091B (zh) 一种提取逆时偏移角道集的方法
Broadhead et al. Performance of some sparseness criterion blind deconvolution methods in the presence of noise
US11322169B2 (en) Target sound enhancement device, noise estimation parameter learning device, target sound enhancement method, noise estimation parameter learning method, and program
CN109425892B (zh) 地震子波的估计方法及***
US10871585B2 (en) System for processing seismic data based upon linear optimization and related methods
CN107561578B (zh) 一种针对宽方位地震资料的面波压制方法
De Martino et al. Very-long-period volcanic tremor at Stromboli, Italy
CN101645701B (zh) 一种基于滤波器组的时延估计方法及其***
US20180038973A1 (en) System for processing seismic data based upon volatility measurement model and related methods
Shumway et al. Mixed signal processing for regional and teleseismic arrays
Zhou et al. A practical methodology for enhancement and detection of transient faults in a gearbox without prior fault feature information

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