地震勘探资料采集中精确定位干扰源位置的方法技术

本发明专利技术公开了一种地震勘探资料采集中精确定位干扰源位置的方法，包括如下步骤：在地震勘探野外采集现场，记录外界干扰源产生的干扰波；提取干扰波同一波前面对应的时间；计算三维空间中一个点震源即干扰源到观测点的距离、方位，实现对外界干扰源位置的定位。本发明专利技术通过精确定位干扰源位置，从而实现快速控制采集时外界干扰源产生的噪音，提高资料品质和施工效率。

【技术实现步骤摘要】

本专利技术涉及一种，属于地震勘探资料采集

技术介绍
长期以来，地震勘探资料采集都受到来自外界干扰源产生的环境噪音的影响。定位外界干扰源位置只能采取人工实地寻访的方法，一个20条测线(测线间距500m)、测线总长度200Km的二维地震勘探项目，或勘探面积为IOOKm2三维地震勘探项目，采取人工实地寻访干扰源，需要排出15人次、5车次在崎岖陡峭的山路上寻访，耗时一周左右方能完成； 若工区交通不便完成时间则会相应更长。这种方法不仅效率低、成本高，且存在明显的疏漏，难于发现断续工作或临时性作业新增干扰源，如为灌溉临时启动的抽水机、工矿里间歇运行的机器等干扰。那些疏漏的干扰源势必会影响采集资料的质量，正常采集作业不得不面对临时中断，临时协停干扰，这不仅延长了正常采集作业时间，工期延后，降低了采集施工效率；而且需要增加作业人员的数量、加大了车辆使用频率，造成人力、物力的浪费。在采集过程中，为了准确定位外界干扰源位置，往往要花费大量人力、物力和时间，极大地制约了采集施工效率，同时外界干扰也严重影响了资料品质，但目前没有一种定位外界干扰源位置的技术能解决采集施工过程面临的难题。
技术实现思路
本专利技术的目的在于克服现有地震勘探资料采集时定位外界干扰源位置存在的上述问题，提供一种，本专利技术通过精确定位干扰源位置，从而实现快速控制采集时外界干扰源产生的噪音，提高资料品质和施工效率。为实现上述目的，本专利技术采用的技术方案如下一种，其特征在于，包括如下步骤a、在地震勘探野外采集现场，记录外界干扰源产生的干扰波；b、提取干扰波同一波前面对应的时间；C、计算三维空间中一个点震源即干扰源到观测点的距离、方位，实现对外界干扰源位置的定位。所述a步骤中，采用三维观测系统，布设多条排列进行观测，利用地震仪器记录外界干扰源产生的干扰波。所述b步骤中，波前是指波速分界面上，各质点开始震动的点的面，或者某一时刻介质中各点开始振动形成的曲面为波前；在弹性介质中，已知T时刻的同一波前面上的各个点，看作各个点从该时刻产生子波的新的点震源，子波的包络面就是波Τ+ΔΤ时刻到达的新的波前面。所述b步骤中，同一波前面对应的时间是指干扰源多次震动之中的某一次或同一次震动产生的球面波，经过任何一个Δ T时刻后，被观测和记录。所述b步骤中，外界干扰源产生干扰波，以干扰源为圆心，在弹性介质中以球面波形式向前传播，在不同的地震道拾取干扰波同一波前面相对应的相同时间。所述c步骤中，基于相同时间对应的地震道坐标拟合平面直角坐标系圆曲线，采用最小二乘法拟合圆曲线，求取圆心坐标，求出三维空间中一个点震源到观测点的距离、方位，实现外界干扰源位置的精确定位。所述c步骤中，计算方法为最小二乘法拟合圆曲线R2 = (x-A)2+(y-B)2R2 = x2-2Ax+A2+y2"2By+B2令a = -2Ab = -2Bc = A2+B2-R2得圆曲线方程的另一个形式x2+y2+ax+by+c = 0 (1)求出参数a，b，c，求得圆心半径的参数A- Jh5.=——- 2j《........................................................i = .Ja +b ^ 4c平面直角坐标系中，设圆心坐标为0(xl，yl)，圆上一点P(x2，y2)，已知两点坐标， 计算两点距离及两点方向的方位角；计算公式为S12 = sqr((x2-xl)2+(y2-yl)2) (7)A12 = arcsin((y2-yl)/S12)(8)S12为圆心至圆上一点P的距离；A12为圆心至圆上一点P的坐标方位角；或者A12 = arccos(Ax21/S12)*sgn(Ay21)+360° (9)式中SgnO为取符号函数；由式⑴拟合圆，据式(7)、⑶或(9)求出三维空间中一个点震源到观测点的距离、方位，实现外界干扰源位置的精确定位。采用本专利技术的优点在于一、本专利技术通过精确、高效地定位干扰源位置，从而实现快速控制采集时外界干扰源产生的噪音，保障了资料品质，提高施工效率。二、与现有的人工实地寻方的方法相比，本专利技术投入的人力、物力和时间较少，因而极大地降低了成本。三、本专利技术有利于发现断续工作或临时性作业新增干扰源，探测干扰源不存在任何的疏漏和盲区。四、正常采集作业不会因新增干扰源临时中断，避免了临时协停干扰造成的正常采集作业时间延长，确保了工期。五、避免了临时协停干扰需要增加作业人员，降低了车辆使用频率，提高了人力、 物力的使用能效。具体实施例方式以下对本专利技术作展开说明一种，包括如下步骤a、在地震勘探野外采集现场，采用三维观测系统，布设多条排列进行观测，利用地震仪器记录外界干扰源产生的干扰波；b、提取干扰波同一波前面对应的时间。波前是波速分界面上，各质点开始震动的点的面；或者说，某一时刻介质中各点刚好开始振动，这一曲面叫波前，也叫波阵面。在弹性介质中，已知T时刻的同一波前面上的各个点，可以看作各个点从该时刻产生子波的新的点震源，这些子波的包络面就是波T+Δ T时刻到达的新的波前面。外界干扰源产生干扰波，以干扰源为圆心，在弹性介质中以球面波形式向前传播， 同一波前面对应的时间是指干扰源多次震动之中的某一次(或同一次)震动产生的球面波，经过任何一个Δ T时刻后，被观测并记录下来。依据干扰波的动力学特征，主要是能量特性，一般包括振幅、频率、相位和衰减等特性进行定量分析，正确识别干扰源多次震动之中的同一次震动产生的球面波波前面。依据地震波以球面波形式进行传播的理论，同一波前面对应的时间值具有等时性。在不同的地震道拾取干扰波同一波前面相对应的相同时间。C、计算三维空间中一个点震源(干扰源)到观测点的距离、方位。基于相同时间对应的地震道坐标拟合平面直角坐标系圆曲线。采用最小二乘法拟合圆曲线，求取圆心坐标。最小二乘法(least squares analysis)是一种数学优化技术，它通过最小化误差的平方和找到一组数据的最佳函数匹配。最小二乘法是用最简的方法求得一些绝对不可知的真值，而令误差平方之和为最小。最小二乘法通常用于曲线拟合(least squares fitting)ο 最小二乘法拟合圆曲线R2= (x-A)2+(y-B)2 R2 = x2-2Ax+A2+y2-2By+B2 令a = _2k b = -2B c = A2+B2-R2可得圆曲线方程的另一个形式x2+y2+ax+by+c = 0 (1)只要求出参数a，b，c就可以求得圆心半径的参数即可得至I最小二乘法(least squares analysis)拟合圆曲线的公式推导如下 样本集(Xi^Yi)I e (1，2，3···Ν)中点到圆心的距离为(Ii Cli2 = (Xi-A)2+ (Yi-B)2点(XyYi)到圆边缘的距离的平方与和半径平方的差为 δ j = Cli2-R2 = (Xi-A)2+ (Yi-B) 2-R2 = Xi'+Y^+aXi+bYi+c 令Q(a，b，c)为Si的平方和 Q(a，b，c) =Σ Si2= Σ 2 求参数a，b，c使得Q(a，b，c)的值最小值； 解平方差Q(a，b，c)大于0，因此函数存在大于或等于0的极小值，极大值为无穷大。 F(a, b，c)对a，b，c求偏导，令偏本文档来自技高网...

【技术保护点】
１．一种地震勘探资料采集中精确定位干扰源位置的方法，其特征在于，包括如下步骤：ａ、在地震勘探野外采集现场，记录外界干扰源产生的干扰波；ｂ、提取干扰波同一波前面对应的时间；ｃ、计算三维空间中一个点震源即干扰源到观测点的距离、方位，实现对外界干扰源位置的定位。

【技术特征摘要】

【专利技术属性】
技术研发人员：文中平，张孟，何光明，李正佳，司阳涛，刘泽斌，岳礼，杨智超，管敏，周晓冀，
申请(专利权)人：中国石油天然气集团公司，中国石油集团川庆钻探工程有限公司，
类型：发明
国别省市：11