一种基于检波器权重因素的微地震加权定位方法技术

本发明专利技术涉及一种基于检波器权重因素的微地震加权定位方法，具体来说就是通过对地面检波器对定位影响的相关因素分析计算后给予不同的权重，通过赋权提高微地震事件定位精度。本发明专利技术是在分析影响检波器定位权重因子的基础上，通过对震中距和信噪比分配不同的权重，能够在一定程度上提升微地震事件定位精度，并有效提高射孔点周围微地震定位可信度。经试验证明：与现有方法相比，本发明专利技术的方法能够通过加权将射孔事件重定位至其真实值处，并能有效提高射孔点附近微地震事件定位的可信度。

【技术实现步骤摘要】
一种基于检波器权重因素的微地震加权定位方法
本专利技术涉及一种地球物理探测方法，基于检波器权重因素的微地震加权定位方法，具体来说是通过对地面检波器对定位影响的相关因素进行分析计算后给予不同的权重，通过调整权重达到提高微地震事件定位精度的目的。技术背景微地震事件定位技术已经广泛应用于油页岩开采中，而现阶段主要有两个因素在影响着定位：一个是现阶段速度模型不准确，一个是地震波到时不能十分精确的拾取。那么当地震速度模型相对准确的情况下，只考虑检波器本身接收到的地震波信号，分析影响其的相关因素，进而根据各个因素对定位影响的强弱，对不同的因素赋予不同的权重，从而实现定位精度的提升也是值得研究的一个课题。CN201511032278公开了一种《提高地震数据信噪比的方法及装置》，从增强信噪比入手，基于倾斜叠加变换，获取第一炮检距向量片道集数据的属性信息；基于每个第一炮检距向量片道集数据的属性信息和预设权重值，这样通过预设权重，达到提高地震数据信噪比的目的。CN201510672801公开了一种《一种地震属性的优化方法和装置》，提出提取地震数据中待优化地震属性的权重值，将待优化地震属性的权重值进行排序。进而通过对提取的待优化地震属性的权重数据进行定量分析，最后利用数学方法计算得到待优化地震属性的权重值，根据权重值的大小，取出权重值小的地震属性，保留权重值大的地震属性。但上述方法仅考虑到从数据属性提出预设权重，或者根据地震属性剔除相关权重小的属性，没有对与检波器相关的因子进行分析。因此，如何从微地震定位中的检波器影响因素出发，对不同因素赋予不同的权重是一个全新的提升定位精度的方法。
技术实现思路
本专利技术在基于网格搜索的振幅叠加微地震定位方法的基础上，分析了地面检波器的信噪比与震中距这两个因子对微地震定位精度的影响，最后利用模拟退火算法对每个检波器定位权重进行分配，通过调整检波器振幅来提高射孔事件重定位精度。本专利技术的目的是通过以下技术方案实现的：一种基于检波器权重因素的微地震加权定位方法，包括以下步骤：a、以射孔点位置作为中心点定义目标区域；b、给出第i个站的权重wi；c、读取微地震数据；d、对射孔定位并获得能量聚焦值最大的网格中心点坐标；e、以单道信噪比和震中距作为不确定因素，采用模拟退火法对目标层进行调整；f、对射孔事件进行重定位；g、是否满足射孔精度要求，是；h、结束。步骤b所述的第i个站的权重wi为：wi＝aiwdi+biwsnri(1)其中ai为bi第i个采集站的震中距a，i,bi为信噪比所占权系数，wdi为震中距所占权重，wsnri为信噪比所占权重，ai＝(a1,a2,a3,...,an)，bi＝(b1,b2,b3,...,bn)，求解各检波器震中距权值ai,利用目标函数E(a)求得最佳矩阵a，目标函数表达式：其中A为第i个检波器在第j时刻的振幅大小，M为检波器个数，L为时窗长度。步骤e所述的模拟退火法对目标层进行调整；利用模拟退火算法求取模拟退火初始温度T0，首先为初始温度赋一个很小的正值，多次乘以1.2，使之慢慢升温，直到恰好满足对任何矩阵a值的接收概率接近于1为止，降温参数Tk，降温公式为：Tk＝T0exp(-ck1/2N)(3)其中，k为迭代次数，T0为初始退火参数，c为给定常数，令c＝0.5，N为需要调整的检波器个数；求解随机变量x，用来调整各检波器震中距权值a，Vi为调整公式为：其中Ai为各检波器震中距调整最大上限值，由该道检波器所获信号信噪比及检波器距离射孔点远近决定的，产生x值表达式如下：其中sgn为符号函数，μ取值范围在[-1,1]之间，求解替代概率参数P，当E(a′)＜E(a)时，a′替代a作为当前最优解，当E(a′)＞E(a)时以概率替换当前最优解。有益效果：以往加权只考虑属性问题，剔除或者弱化某些属性然后进行归一化操作，没有从检波器方面进行分析。本专利技术是在分析影响检波器定位权重因子的基础上，通过对震中距和信噪比分配不同的权重，能够在一定程度上提升微地震事件定位精度，并有效提高射孔点周围微地震定位可信度。附图说明：图1一种基于检波器权重因素的微地震加权定位方法流程图；图2三维层状起伏结构地层模型；图3建立初始平层结构射线追踪示意图；图4检波器各道正演模拟波形图；图5(a),(b)为采用检波器不加权定位方法结果图；图5(c),(d)为采用检波器加权定位方法结果图；图6各层界面层厚约束范围以及各层旋转角约束范围。具体实施方式：下面结合附图和实施例对本专利技术作进一步的详细说明：一种基于检波器权重因素的微地震加权定位方法，包括以下步骤：a、以射孔点位置作为中心点定义目标区域；b、给出第i个站的权重wi；c、读取微地震数据；d、对射孔定位并获得能量聚焦值最大的网格中心点坐标；e、以单道信噪比和震中距作为不确定因素，采用模拟退火法对目标层进行调整；f、对射孔事件进行重定位；g、是否满足射孔精度要求，是；h、结束。步骤b所述的第i个站的权重wi为：wi＝aiwdi+biwsnri(1)其中ai为bi第i个采集站的震中距a，i,bi为信噪比所占权系数，wdi为震中距所占权重，wsnri为信噪比所占权重，ai＝(a1,a2,a3,...,an)，bi＝(b1,b2,b3,...,bn)，求解各检波器震中距权值ai,利用目标函数E(a)求得最佳矩阵a，目标函数表达式：其中A为第i个检波器在第j时刻的振幅大小，M为检波器个数，L为时窗长度。步骤e所述的模拟退火法对目标层进行调整；利用模拟退火算法求取模拟退火初始温度T0，首先为初始温度赋一个很小的正值，多次乘以1.2，使之慢慢升温，直到恰好满足对任何矩阵a值的接收概率接近于1为止，降温参数Tk，降温公式为：Tk＝T0exp(-ck1/2N)(3)其中，k为迭代次数，T0为初始退火参数，c为给定常数，令c＝0.5；N为需要调整的检波器个数；求解随机变量x，用来调整各检波器震中距权值a，Vi为调整公式为：其中Ai为各检波器震中距调整最大上限值，由该道检波器所获信号信噪比及检波器距离射孔点远近决定的，产生x值表达式如下：其中sgn为符号函数，μ取值范围在[-1,1]之间，求解替代概率参数P，当E(a′)＜E(a)时，a′替代a作为当前最优解，当E(a′)＞E(a)时以概率替换当前最优解。这时就能得到震中距与信噪比在系统中所占权重值矩阵，然后利用振幅叠加进行定位。实施例1：下面结合附图和实施例对本专利技术进行清楚、完整的描述。具体实施步骤如下：如图2所示，利用三维层状起伏结构地层模型原定位方法，实际信号采用有限差分波动方程进行模拟，根据所模拟的地层模型建立初始平层状地层模型，分为4层，其中1～4层的速度值如表1所示，表1各层界面层厚约束范围以及各层旋转角约束范围依次为：1700m/s,2500m/s,3400m/s,5400m/s，模拟射孔位置为(673，194，-1920)，单位为米。a、在射孔附近定义一个三维目标区域，即该区域有可能发生微震事件。在本实施例中，以射孔坐标为基准，选定目标区域。选定的目标区域的三维坐标为X∈[-800,800]，Y∈[-200,200],Z∈[-2100,-1900]，当定位结果发生在该区域内认为定位到事件。微地震地面检波器成星型形本文档来自技高网...

【技术保护点】
1.一种基于检波器权重因素的微地震加权定位方法，其特征在于，包括以下步骤：a、以射孔点位置作为中心点定义目标区域；b、给出第i个站的权重wi；c、读取微地震数据；d、对射孔定位并获得能量聚焦值最大的网格中心点坐标；e、以单道信噪比和震中距作为不确定因素，采用模拟退火法对目标层进行调整；f、对射孔事件进行重定位；g、是否满足射孔精度要求，是；h、结束。

【技术特征摘要】
1.一种基于检波器权重因素的微地震加权定位方法，其特征在于，包括以下步骤：a、以射孔点位置作为中心点定义目标区域；b、给出第i个站的权重wi；c、读取微地震数据；d、对射孔定位并获得能量聚焦值最大的网格中心点坐标；e、以单道信噪比和震中距作为不确定因素，采用模拟退火法对目标层进行调整；f、对射孔事件进行重定位；g、是否满足射孔精度要求，是；h、结束。2.按照权利要求1所述的一种基于检波器权重因素的微地震加权定位方法，其特征在于，步骤b所述的第i个站的权重wi为：wi＝aiwdi+biwsnri(1)其中ai为bi第i个采集站的震中距a，i,bi为信噪比所占权系数，ai＝(a1,a2,a3,...,an)，bi＝(b1,b2,b3,...,bn),wdi为震中距所占权重，wsnri为信噪比所占权重，求解各检波器震中距权值ai利用目标函数E(a)求得最佳矩阵a，目标函数表达式：将其中A为第i个检波器在第j时刻的振幅...

【专利技术属性】
技术研发人员：孙锋，张彦峰，陈祖斌，吕昊，江海宇，唐国峰，原杰，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林,22