一种频率域自适应偏振角计算方法及系统技术方案

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

A Frequency Domain Adaptive Polarization Angle Calculating Method and System

The present invention proposes a frequency domain adaptive polarization angle calculation method and system, which includes: transforming three components of each detector into frequency domain Uij (_); calculating the adaptive threshold range of each frequency; identifying abnormal frequency components by using the adaptive threshold range to determine the frequency set Sj; and calculating the polarization angle by using the frequency set Sj and the correlation objective function. The biggest innovation of the present invention is that the polarization angle can be calculated adaptively in the frequency domain. This method can select high quality frequency components in the frequency domain to participate in the calculation and improve the stability of the calculation of the polarization angle.

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

【技术保护点】
1.一种频率域自适应偏振角计算方法，其特征在于，该方法包括：将每个检波器的三个分量变换到频率域Uij(ω)；计算每个频率的自适应阈值范围；利用自适应阈值范围识别异常频率成份，确定频率集合Sj；利用频率集合Sj，通过相关目标函数计算偏振角。

【技术特征摘要】
1.一种频率域自适应偏振角计算方法，其特征在于，该方法包括：将每个检波器的三个分量变换到频率域Uij(ω)；计算每个频率的自适应阈值范围；利用自适应阈值范围识别异常频率成份，确定频率集合Sj；利用频率集合Sj，通过相关目标函数计算偏振角。2.根据权利要求1所述的频率域自适应偏振角计算方法，其特征在于，将每个检波器的三个分量变换到频率域中，记作Uij(ω)＝|FFT(uij(t))|，其中，i表示三个分量，取值1，2，3，分别对应两个水平分量和垂直分量，j表示检波器的序号1,…,N。3.根据权利要求2所述的频率域自适应偏振角计算方法，其特征在于，设纵波的偏振角为则旋转完成后切向分量uT、径向分量uR和垂向分量uZ为：其中uE、uN为三分量检波器的两个水平分量，uZ为三分量检波器的垂向分量。4.根据权利要求1所述的频率域自适应偏振角计算方法，其特征在于，计算每个频率的自适应阈值范围包括：统计每个频率成份的平均值：统计每个频率成份的方差：5.根据权利要求4所述的频率域自适应偏振角计算方法，其特征在于，每个频率的自适应阈值为：Th(ω)＝E(ω)+σ(ω)(4)。6.根据权利要求5所述的频率域自适应偏振角...

【专利技术属性】
技术研发人员：崔树果，
申请(专利权)人：中国石油化工股份有限公司，中国石油化工股份有限公司石油物探技术研究院，
类型：发明
国别省市：北京,11