本发明涉及油气和矿产资源勘探技术领域，特别涉及一种井周岩石非均匀各向异性的识别方法及系统。
背景技术
地下岩石普遍存在非均匀性。由于岩石圈中岩石的成分、密度、相态和缝洞中流体的组成、温度和压力等不均匀性，实际地层介质弹性参数，如弹性模量、速度和泊松比等，在空间分布上往往是非均匀的。岩石非均匀的各向异性识别对于油气和矿产资源的精细勘探至关重要。
测井尺度的岩石非均匀各向异性可以通过不同方位的时差差异、层析成像等方法进行识别，这些方法利用的主要是声波测井数据中的直达波信息，只能评价井壁附近的岩石非均匀各向异性。目前对于远离井壁的岩石非均匀各向异性的识别方法研究较少。事实上，当远离井壁的岩石存在非均匀各向异性时，不同方位会产生不同强度的散射波，可以利用不同方向声波测井数据中散射波的强弱识别岩石非均匀各向异性。而在声波测井常规数据处理中，散射波和尾波往往被视为噪声而加以压制。现代声波测井仪器已广泛应用于油气勘探，采集了大量正交偶极的阵列声波测井数据。如果能从声波测井数据正确分析和提取岩石的非均匀各向异性信息，将会提供很有价值的地层信息，并为声波测井的处理和解释开辟新的应用前景。
技术实现要素：
针对现有技术的不足，本发明提供一种井周岩石非均匀各向异性的识别方法及系统，计算不同方位偶极声波散射衰减品质因子识别岩石非均匀各向异性，为计算声波测井频率下井周岩石非均匀各向异性识别提供一种实用可行的方法和技术。
为实现上述目的，本发明提供一种井周岩石非均匀各向异性的识别方法，所述方法包括以下步骤：
s1、在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据；
s2、根据步骤s1获取深度区间多个深度位置的正交偶极四分量声波测井数据，并记录每一个深度位置处的仪器方位；
s3、将步骤s2中获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
s4、将方位为
s5、利用希尔伯特变换计算v(t)的能量包络a(t)；
s6、创建目标函数δe，所述目标函数δe定义为理论的能量包络ac(t)和实测的能量包络a(t)的函数：
δe(t,qi,qs)＝|ac(t)-a(t)|2
其中qs为散射衰减品质因子，qi为固有衰减品质因子；
s7、根据步骤s6求解目标函数，得到方位为
s8、重复步骤s3～步骤s7，获得深度区间内岩石非均匀各向异性结果。
在上述技术方案中，步骤s2中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
在上述技术方案中，步骤s4中所述方位
在上述技术方案中，步骤s6中所述理论的能量包络ac(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，id是直达波振幅ad在时间窗(t1,t2)的积分，表示为
在上述技术方案中，步骤s7所述集合
s71、若集合
s72、若集合
本发明还提供一种井周岩石非均匀各向异性的识别系统，其特征在于，所述系统包括采集模块、记录模块、转换模块、滤波模块、计算模块、创建模块、判断模块和循环模块；
采集模块，用于在深度区间内进行陈列声波测井，采集深度区间内正交偶极四分量陈列波形数据；
记录模块，用于在深度区间多个深度位置的正交偶极四分量声波测井数据的获取点，记录深度位置处的仪器方位；
转换模块，用于将获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
滤波模块，用于将方位为
计算模块，用于运用希尔伯特变换计算v(t)的能量包络a(t)；
创建模块，创建目标函数δe，所述目标函数δe定义为理论的能量包络ac(t)和实测的能量包络a(t)的函数：
δe(t,qi,qs)＝|ac(t)-a(t)|2
其中qs为散射衰减品质因子，qi为固有衰减品质因子；
判断模块，用于求解目标函数，得到方位为
循环模块，用于循环重复转换模块～判断模块，获得深度区间内岩石非均匀各向异性结果。
在上述技术方案中，记录模块中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
在上述技术方案中，滤波模块中所述方位
在上述技术方案中，创建模块中所述理论的能量包络ac(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，id是直达波振幅ad在时间窗(t1,t2)的积分，表示为
在上述技术方案中，判断模块中所述集合
相同模块，若集合
相异模块，若集合
本发明一种井周岩石非均匀各向异性的识别方法及系统，具有以下有益效果：本发明采用不同方位散射衰减品质因子的差异来识别岩石非均匀各向异性，与现有技术方案只能判断井壁附近岩石非均匀各向异性相比，本发明可判断远离井壁岩石的非均匀各向异性。岩石非均匀各向异性的结果可用于指导油气开发如井位优化、酸化压裂等，为油气的勘探、开采和生产提供重要的信息。
附图说明
图1为本发明一种井周岩石非均匀各向异性的识别方法流程图；
图2为本发明实施例计算得到的某一深度位置不同方位散射衰减品质因子示意图；
图3为本发明一种井周岩石非均匀各向异性的识别系统模块图；
具体实施方式
下面结合附图对本发明作进一步详细描述，本发明提供一种井周岩石非均匀各向异性的识别方法，所述方法包括以下步骤，如图1所示：
s1、在深度区间内进行陈列声波测井，获取深度区间内正交偶极四分量陈列波形数据；
所述正交偶极四分量阵列数据是用多极子测井仪器的偶极模式采集的数据，正交偶极是相互垂直的两组换能器x和y，与x1和y1发射对应的接收x2和y2，采集的数据为四组数据分别为：x1发x2收、x1发y2收、y1发x2收、y1发y2收。
s2、根据步骤s1获取深度区间多个深度位置的正交偶极四分量声波测井数据，并记录每一个深度位置处的仪器方位；
其中，仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
s3、将步骤s2中获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
其中，通过四组数据组合计算得到不同方位的偶极声波测井数据即为横波数据。
s4、将方位为
其中，所述方位
s5、利用希尔伯特变换计算v(t)的能量包络a(t)；
s6、创建目标函数δe，所述目标函数δe定义为理论的能量包络ac(t)和实测的能量包络a(t)的函数：
δe(t,qi,qs)＝|ac(t)-a(t)|2
其中qs为散射衰减品质因子，qi为固有衰减品质因子；
其中，声波信号的散射衰减品质因子是描述岩石非均匀性的重要参数，不同方位声波测井信号的散射衰减品质因子的差异预示着岩石的非均匀各向异性。
具体的，所述理论的能量包络ac(t)可由三维非均匀模型的能量包络表达式如下计算得到：
其中，f为频率，td为直达波走时，id是直达波振幅ad在时间窗(t1,t2)的积分，表示为
s7、根据步骤s6求解目标函数，得到方位为
具体的，所述集合
s71、若集合
s72、若集合
s8、重复步骤s3～步骤s7，获得深度区间内岩石非均匀各向异性结果。
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明；
如图2所示利用本发明的处理步骤计算得到的某一深度位置不同方位(0度～360度)的散射衰减品质因子，从图中可看出在45度方位的散射衰减品质因子最大(约为435)，150度方位对应的散射衰减品质因子最小(约为145)，利用步骤s7的判断条件，该深度点岩石具有非均匀各向异性。
如图3所示本发明还提供一种井周岩石非均匀各向异性的识别系统，所述系统包括采集模块、记录模块、转换模块、滤波模块、计算模块、创建模块、判断模块和循环模块；
采集模块，用于在深度区间内进行陈列声波测井，采集深度区间内正交偶极四分量陈列波形数据；
记录模块，用于在深度区间多个深度位置的正交偶极四分量声波测井数据的获取点，记录深度位置处的仪器方位；
其中，记录模块中仪器方位为所述任一深度位置的正交偶极四分量声波测井数据的任一个偶极相对于地球北极的方位。
转换模块，用于将获取的任一深度位置偶极四分量测井数据由仪器坐标转换到地球坐标系下，构建不同方位的偶极声波测井数据；
滤波模块，用于将方位为
其中，所述方位
计算模块，用于运用希尔伯特变换计算v(t)的能量包络a(t)；
创建模块，创建目标函数δe，所述目标函数δe定义为理论的能量包络ac(t)和实测的能量包络a(t)的函数：
δe(t,qi,qs)＝|ac(t)-a(t)|2
其中qs为散射衰减品质因子，qi为固有衰减品质因子；
具体的，所述理论的能量包络ac(t)可由三维非均匀模型的能量包络表达式如下计算得到：
判断模块，用于求解目标函数，得到方位为
具体的，所述集合
相同模块，若集合
相异模块，若集合
循环模块，用于循环重复转换模块～判断模块，获得深度区间内岩石非均匀各向异性结果。
说明书中未阐述的部分均为现有技术或公知常识。本实施方式仅用于说明该发明，而不用于限制本发明的范围，本领域技术人员对于本发明所做的等价置换等修改均认为是落入该发明权利要求书所保护范围内。