三分量检波器水平分量的方向检测方法和装置制造方法及图纸

本发明专利技术提供了一种三分量检波器水平分量的方向检测方法和装置，其中，该方法包括：获取三分量地震数据以及三分量地震数据的初至时间；根据三分量地震数据、初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；计算待建立协方差矩阵的主特征向量；计算主特征向量在水平面内的投影与横坐标轴的第一夹角；根据第一夹角、炮检连线的方位角，计算得到第一夹角与方位角之间的第二夹角。在本发明专利技术实施例中，达到了对三分量检波器的方向进行有效分析和质控的目的。

【技术实现步骤摘要】
三分量检波器水平分量的方向检测方法和装置
本专利技术涉及地质勘探
，特别涉及一种三分量检波器水平分量的方向检测方法和装置。
技术介绍
由于施工条件的不确定性，采集多分量地震数据时，难以保证完全按照采集的要求将每个三分量检波器的Z分量都置于垂直方向，将每个三分量检波器的X分量都置于平行于测线的方向。在这种情况下，采集到的地震数据在三分量检波器的每个分量上都会存在能量投影。具体而言，在Z分量上有横波能量，在水平分量上有纵波能量。采集数据的不准确会给波场的分离、成像和反演带来困难。因此，如何有效地监测三分量检波器接收地震数据的状态是多分量地震勘探的关键。目前，三分量检波器的接收状态可以通过以下方式来监测：通过分析三分量地震记录的初至波来估算三分量检波器的方向，并监测三分量检波器的接收状态。对于垂直地震剖面数据、没有记录检波器方向参数的海上多分量数据而言，可以利用多分量地震记录的初至波估算三分量检波器的方向参数，然后对三分量地震数据进行重定向。目前，可以利用初至波矢端曲线来监测检波器重定向前后初至波的变化情况，其中，矢端曲线可以是地震波传播时地下介质中质点振动随时间变化的空间轨迹图形。因此，利用矢端曲线估算三分量检波器的方向时，只能估算到初至波质点的振动方向。
技术实现思路
本专利技术提供了一种三分量检波器水平分量的方向检测方法和装置，以解决现有技术中仅能估算到初至波质点振动方向的问题。本专利技术实施例提供了一种三分量检波器水平分量的方向检测方法，可以包括：获取三分量地震数据以及所述三分量地震数据的初至时间；根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；计算所述待建立协方差矩阵的主特征向量；计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角；根据所述第一夹角、炮检连线的方位角，计算得到所述第一夹角与所述方位角之间的第二夹角，所述方位角是根据三分量检波器的位置以及炮点的位置确定的。在一个实施例中，所述方法还包括：基于所述第一夹角以及所述第二夹角，对所述三分量检波器的水平分量进行监测。在一个实施例中，所述计算所述待建立协方差矩阵的主特征向量包括，可以包括：计算所述待建立协方差矩阵的多个特征值；将所述多个特征值中的最大值作为所述待建立协方差矩阵的最大特征值；将所述最大特征值对应的特征向量作为所述待建立协方差矩阵的主特征向量。在一个实施例中，基于所述第一夹角以及所述第二夹角，对所述三分量检波器的方向进行调整，包括：基于所述三分量检波器的位置以及所述炮点的位置，确定待绘制图形的横坐标范围以及纵坐标范围，并在所述待绘制图形中确定出所述三分量检波器位置以及所述炮点位置；在所述待绘制图形中绘制长度等于所述最大特征值的线段，所述确定出的三分量检波器位置位于所述线段中点，所述线段与所述横坐标轴的夹角等于所述第一夹角；以所述待绘制图形中线段的颜色深浅表示所述第二夹角；基于所述线段与所述横坐标轴的夹角以及所述颜色深浅，对所述三分量检波器的方向进行调整。在一个实施例中，按照以下公式确定所述炮检连线的方位角：上式中，θ表示所述炮检连线的方位角，XR表示所述三分量检波器位置所对应的横坐标，YR表示所述三分量检波器位置所对应的纵坐标，XS表示所述炮点位置所对应的横坐标，YS表示所述炮点位置所对应的纵坐标。在一个实施例中，按照以下公式根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量：上式中，t1＝Tfb(r)/Δt，t2＝(Tfb(r)+Tw)/Δt，n＝t2-t1+1；上式中，表示所述参量，Si(t)和Sj(t)表示所述三分量地震数据，其中，i＝1,2,3，j＝1,2,3，i和j分别表示所述三分量地震数据的分量序号，Tfb(r)表示三分量地震数据中一个分量第r道的初至时间，r＝1,2，...,Ntr，Ntr表示所述三分量地震数据中一个分量的总道数，Tw表示所述预设的时窗长度，Δt表示所述三分量地震数据的采样间隔。在一个实施例中，计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角，包括：对所述主特征向量进行归一化处理；计算归一化处理后的主特征向量在水平面内的投影与横坐标轴的第一夹角。在一个实施例中，按照以下公式计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角：上式中，φ表示所述第一夹角，Wii表示所述主特征向量中的第ii个分量，ii＝1,2。在一个实施例中，按照以下公式计算得到所述第一夹角与所述方位角之间的第二夹角：α＝φ-θ上式中，α表示所述第二夹角，φ表示所述第一夹角，θ表示所述炮检连线的方位角。在一个实施例中，所述预设的时窗长度是根据所述三分量地震数据的地震子波确定的。本专利技术实施例还提供了一种三分量检波器水平分量的方向检测装置，可以包括：数据获取模块，用于获取三分量地震数据以及所述三分量地震数据的初至时间；参量计算模块，用于根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；向量计算模块，用于计算所述待建立协方差矩阵的主特征向量；第一夹角计算模块，用于计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角；第二夹角计算模块，用于根据所述第一夹角、炮检连线的方位角，计算得到所述第一夹角与所述方位角之间的第二夹角，所述方位角是根据三分量检波器的位置以及炮点的位置确定的。本专利技术实施例还提供了一种三分量检波器水平分量的方向检测装置，可以包括处理器以及用于存储处理器可执行指令的存储器，所述处理器执行所述指令时实现：获取三分量地震数据以及所述三分量地震数据的初至时间；根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；计算所述待建立协方差矩阵的主特征向量；计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角；根据所述第一夹角、炮检连线的方位角，计算得到所述第一夹角与所述方位角之间的第二夹角，所述方位角是根据三分量检波器的位置以及炮点的位置确定的。在本专利技术实施例中，基于三分量地震数据以及初至时间计算得到所述待建立协方差矩阵中的参量，再计算所述待建立协方差矩阵的主特征向量在水平面内的投影与横坐标轴的第一夹角，根据所述第一夹角以及炮检连线的方位角，计算得到表征所述三分量检波器水平分量的方向的所述第二夹角。采用上述方法，解决了现有技术中利用矢端曲线监测三分量检波器状态时，只能表示初至波质点振动方向，而不能对初至波质点振动方向偏离炮检连线方向的程度、初至波的能量变化等进行有效监控的问题，达到了对三分量检波器的方向进行有效分析和质控的目的。附图说明为了更清楚地说明本申请实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请中记载的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这些附图获得其他的附图。图1是本申请提供的一种三分量检波器水平分量的方向检测方法流程图；图2是本申请提供的采用三分量检波器水平分量的方向检测方法绘制某三分量检波器偏移量较大地区的待绘制图形的示意图；图3是本申请提供的采用三分量检波器水平分量的方向检测方法绘制某三分量检波器偏移量较小地区的待绘制图形的示意图；图4是本文档来自技高网...

【技术保护点】
一种三分量检波器水平分量的方向检测方法，其特征在于，包括：获取三分量地震数据以及所述三分量地震数据的初至时间；根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；计算所述待建立协方差矩阵的主特征向量；计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角；根据所述第一夹角、炮检连线的方位角，计算得到所述第一夹角与所述方位角之间的第二夹角，所述方位角是根据三分量检波器的位置以及炮点的位置确定的。

【技术特征摘要】
1.一种三分量检波器水平分量的方向检测方法，其特征在于，包括：获取三分量地震数据以及所述三分量地震数据的初至时间；根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量；计算所述待建立协方差矩阵的主特征向量；计算所述主特征向量在水平面内的投影与横坐标轴的第一夹角；根据所述第一夹角、炮检连线的方位角，计算得到所述第一夹角与所述方位角之间的第二夹角，所述方位角是根据三分量检波器的位置以及炮点的位置确定的。2.如权利要求1所述的方法，其特征在于，在计算得到所述第一夹角与所述方位角之间的第二夹角之后，所述方法还包括：基于所述第一夹角以及所述第二夹角，对所述三分量检波器的水平分量进行监测。3.如权利要求1所述的方法，其特征在于，所述计算所述待建立协方差矩阵的主特征向量包括，包括：计算所述待建立协方差矩阵的多个特征值；将所述多个特征值中的最大值作为所述待建立协方差矩阵的最大特征值；将所述最大特征值对应的特征向量作为所述待建立协方差矩阵的主特征向量。4.如权利要求3所述的方法，其特征在于，基于所述第一夹角以及所述第二夹角，对所述三分量检波器的方向进行调整，包括：基于所述三分量检波器的位置以及所述炮点的位置，确定待绘制图形的横坐标范围以及纵坐标范围，并在所述待绘制图形中确定出所述三分量检波器位置以及所述炮点位置；在所述待绘制图形中绘制长度等于所述最大特征值的线段，所述确定出的三分量检波器位置位于所述线段中点，所述线段与所述横坐标轴的夹角等于所述第一夹角；以所述待绘制图形中线段的颜色深浅表示所述第二夹角；基于所述线段与所述横坐标轴的夹角以及所述颜色深浅，对所述三分量检波器的方向进行调整。5.如权利要求1所述的方法，其特征在于，按照以下公式确定所述炮检连线的方位角：上式中，θ表示所述炮检连线的方位角，XR表示所述三分量检波器位置所对应的横坐标，YR表示所述三分量检波器位置所对应的纵坐标，XS表示所述炮点位置所对应的横坐标，YS表示所述炮点位置所对应的纵坐标。6.如权利要求1所述的方法，其特征在于，按照以下公式根据所述三分量地震数据、所述初至时间、预设的时窗长度，计算得到待建立协方差矩阵中的参量：上式中，t1＝Tfb(r)/Δt，t2＝(Tfb(r)+Tw)/Δt，n＝t2-t1+1；上式中，表示所述参量，S...

【专利技术属性】
技术研发人员：张文波，马强，李建峰，聂红梅，孙鹏远，颉冬莲，
申请(专利权)人：中国石油天然气集团公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11