公开/公告号CN114994753A
专利类型发明专利
公开/公告日2022-09-02
原文格式PDF
申请/专利权人 中国海洋石油集团有限公司;中海石油（中国）有限公司天津分公司;
申请/专利号CN202210456110.3
申请日2022-04-28
分类号G01V1/30(2006.01);G01V1/28(2006.01);
代理机构天津三元专利商标代理有限责任公司 12203;
代理人高凤荣
地址 100010 北京市东城区朝阳门北大街25号
入库时间 2023-06-19 16:46:06
法律状态公告日
法律状态信息
法律状态
2022-09-20
实质审查的生效 IPC(主分类):G01V 1/30 专利申请号:2022104561103 申请日:20220428
实质审查的生效
2022-09-02
公开
发明专利申请公布
技术领域
本发明属于地震勘探处理领域，尤其涉及井约束转换波深度偏移初始横波速度场建立方法。
背景技术
目前，在转换波地震资料处理过程中，常用的是转换波叠前时间偏移技术，叠前时间偏移技术采用的是转换波时距曲线简化方程。这些简化方程或者公式有二三十种，其中，得到转换波处理业界普遍认可的是李向阳博士所撰写的“系列转换波叠前时间偏移公式”。该公式是从转换波时间域出发，采用四参数对转换波ccp道集进行动校正的一套方法，其结果是时间域的叠前时间偏移结果。但是，当有大面积气云存在的情况下，转换波由于有一部分也是受到气云影响，因此，导致气云区普遍同相轴存在向下弯曲的现象，不利于转换波的成像。
对存在大面积气云的转换波地震数据，一般是进行深度偏移来彻底解决转换波速度问题，从而使转换波偏移归位。而转换波速度又包括：纵波速度和横波速度二种；但是，由于纵波速度可以通过纵波速度场层析反演得到，而横波速度用现有的方法则无法得到。因此，无法得到准确的转换波速度场。
发明内容
本发明目的在于提供一种井约束转换波深度偏移初始横波速度场建立方法,以解决求取转换波深度偏移所需的深度域横波层速度场的技术问题。
为实现上述目的，本发明的井约束转换波深度偏移初始横波速度场建立方法的具体技术方案如下：
一种井约束转换波深度偏移初始横波速度场建立方法，采用以下步骤：
第一步：读取纵波、横波声波测井时差曲线；
第二步：抽取纵波测井速度、横波测井速度平滑曲线，并匹配地求取百分比；
第三步：多井百分比曲线内插外推形成百分比体；
第四步：读取网格层析纵波最终速度场，抽取井位置处速度曲线；
第五步：将纵波测井速度平滑曲线和纵波速度匹配求取百分比；
第六步：纵波单井百分比曲线，在深度上外推至顶底；
第七步：纵波多井百分比曲线内插外推形成百分比体；
第八步：比例到多井约束纵波速度体；
第九步：将横波多井百分比体和多井约束纵波速度体结合，得到多井约束横波速度体。
进一步，所述第一步中：读取纵波、横波声波测井时差曲线，先从测井时差曲线的求取百分比系数公式为：
percent＝v2/v1 (1)
测井时差曲线应用百分比系数公式为：
v2＝percent*v1 (2)
测井时差曲线的内插公式为直线内插公式，即：若 A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上；则(b-b1)/(i-i1)＝(b2-b1)/(i2-i1)＝直线斜率，变换即得所求；
测井时差曲线的外推公式为拟合公式，即：将已有样点代入 y＝ax
测井时差曲线的抽取平滑曲线公式为：拟合公式，即：将已有样点代入y＝ax
进一步，所述第二步中：抽取纵波测井速度、横波测井速度平滑曲线的步骤如下：
1)在纵波声波时差曲线上抽取平滑曲线；抽取平滑曲线需要根据实际需要的精度进行测试；抽取平滑曲线公式为拟合公式，即：将已有样点代入y＝ax
其中，单井百分比曲线在深度上，外推至顶底；外推公式为拟合公式，即将已有样点代入y＝ax
2)在横波声波时差曲线上抽取平滑曲线；抽取平滑曲线需要根据实际需要的精度进行测试；抽取平滑曲线公式为拟合公式，即将已有样点代入y＝ax
其中，单井百分比曲线在深度上，外推至顶底；外推公式为拟合公式，即将已有样点代入y＝ax
3)纵横波匹配求取百分比；
纵波速度和横波速度相同深度的样点值相除求取百分比系数。 percent＝vpw/vsw，其中，vpw为纵波测井速度，vsw为横波测井速度，percent为百分比系数。
进一步，所述第三步中：多井百分比曲线内插外推形成百分比体的步骤如下：
1)外推公式为拟合公式，即：将已有样点代入y＝ax
2)内插公式为直线内插公式，即：若A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上；则(b-b1)/(i-i1)＝ (b2-b1)/(i2-i1)＝直线斜率，变换即得所求；所求百分比体记为 percent1。
进一步，所述第五步中：其纵波测井速度平滑曲线和纵波速度相同深度的样点值相除求取百分比系数为；percent2＝vpw/vp，其中， vpw为：纵波测井平滑速度曲线样点速度值，vp为：网格层析最终纵波速度场提取井位置处速度曲线速度值，percent2为：百分比系数。
进一步，所述第六步中：其外推公式为拟合公式，即：将已有样点代入y＝ax
进一步，所述第七步中：其步骤如下：
1)内插公式为直线内插公式，即：若A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上；则(b-b1)/(i-i1)＝ (b2-b1)/(i2-i1)＝直线斜率，变换即得所求；
2)外推公式为拟合公式，即：将已有样点代入y＝ax
进一步，所述第八步中：其步骤如下：
使用上面百分比体percent2比例纵波网格层析最终速度场到多井约束纵波速度体：
使用公式为vpw＝vp*percent2，其中，vpw为多井约束纵波速度体，vp为纵波网格层析最终速度场,percent2为百分比系数。
进一步，所述第九步中：其所使用公式为：vsw＝vpw*percent1，其中，vsw为：多井约束横波速度体，vpw为：井约束纵波速度体, percent1为第三步所求百分比体。
本发明的井约束转换波深度偏移初始横波速度场建立方法具有以下优点：
1.本发明对存在大面积气云的转换波地震数据，直接采用转换波深度偏移技术。
2.本发明采用射线追踪技术直接用纵波速度追踪转换波的转换点，通过横波速度追踪到检波器，即采用真实的纵波和横波速度进行追踪，计算旅行时，进而开展深度偏移来彻底解决转换波速度问题，使转换波偏移准确归位。
3.本发明的纵波速度，不仅可以通过纵波速度场层析反演得到，而且，横波速度还能够得到较为精确的初始速度。
4.本发明的纵波速度通过转换波层析反演，便可进一步进行优化。
5.本发明通过转换波深度偏移，便能够提供准确的转换波速度场。
附图说明
下面结合附图和实施例对本发明进一步说明。
图1为本发明的井约束求取的横波速度偏移得到的偏移剖面示意图(其为屏幕上的实际图形)；
图2为本发明求取横波速度流程示意图；
图3为本发明纵波、横波声波抽取平滑曲线示意图(其为屏幕上的实际图形)；
图4为本发明网格层析最终纵波速度场示意图(其为屏幕上的实际图形)；
图5为本发明多井约束横波速度体示意图(其为屏幕上的实际图形)。
具体实施方式
为了更好地了解本发明的目的、结构及功能，下面结合附图，对本发明一种井约束转换波深度偏移初始横波速度场建立方法做进一步详细的描述。
如图1，图2所示，本发明在已知最终的纵波深度域层速度场、纵波声波测井曲线、横波声波测井曲线的情况下，通过测井资料约束，从而，求得初始的转换波深度偏移所需的深度域的横波初始速度场。通过进行以下一系列的处理，达到求取转换波深度偏移，所需的深度域横波层速度场的目的。本发明采用以下步骤：
第一步：读取纵波、横波声波测井时差曲线；
读取纵波、横波声波测井时差曲线时，先从测井时差曲线的求取百分比系数公式为：percent＝v2/v1 (1)
测井时差曲线应用百分比系数公式为：
v2＝percent*v1 (2)
测井时差曲线的内插公式为直线内插公式，即：若 A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上；则(b-b1)/(i-i1)＝(b2-b1)/(i2-i1)＝直线斜率，变换即得所求；
测井时差曲线的外推公式为拟合公式，即：将已有样点代入 y＝ax
测井时差曲线的抽取平滑曲线公式为：拟合公式，即：将已有样点代入y＝ax
如图3所示，图3中所显示的为：测井曲线(左边为纵波速度，右边为横波速度，黑色为实测值，左、右中间的竖黑色为平滑后曲线)；
第二步：抽取纵波测井速度、横波测井速度平滑曲线，并匹配地求取百分比；抽取纵波测井速度、横波测井速度平滑曲线的步骤如下：
1)在纵波声波时差曲线上抽取平滑曲线；抽取平滑曲线需要根据实际需要的精度进行测试；抽取平滑曲线公式为拟合公式，即：将已有样点代入y＝ax
其中，单井百分比曲线在深度上，外推至顶底；外推公式为拟合公式，即将已有样点代入y＝ax
2)在横波声波时差曲线上抽取平滑曲线；抽取平滑曲线需要根据实际需要的精度进行测试；抽取平滑曲线公式为拟合公式，即：将已有样点代入y＝ax
其中，单井百分比曲线在深度上，外推至顶底；外推公式为拟合公式，即将已有样点代入y＝ax
3)纵横波匹配求取百分比；
纵波速度和横波速度相同深度的样点值相除求取百分比系数。 percent＝vpw/vsw，其中，vpw为纵波测井速度，vsw为横波测井速度，percent为百分比系数。
第三步：多井百分比曲线内插外推形成百分比体；
多井百分比曲线内插外推形成百分比体的步骤如下：
1)外推公式为拟合公式，即：将已有样点代入y＝ax
2)内插公式为直线内插公式，即若A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上。则(b-b1)/(i-i1)＝ (b2-b1)/(i2-i1)＝直线斜率，变换即得所求。所求百分比体记为 percent1。
第四步：读取网格层析纵波最终速度场，抽取井位置处速度曲线；如图4所示，其为网格层析最终纵波速度场示意展示；
第五步：将纵波测井速度平滑曲线和纵波速度匹配求取百分比；其纵波测井速度平滑曲线和纵波速度相同深度的样点值相除求取百分比系数为：percent2＝vpw/vp，其中，vpw为：纵波测井平滑速度曲线样点速度值，vp为：网格层析最终纵波速度场提取井位置处速度曲线速度值，percent2为：百分比系数。
第六步：纵波单井百分比曲线，在深度上外推至顶底；
其外推公式为拟合公式，即：将已有样点代入y＝ax
第七步：纵波多井百分比曲线内插外推形成百分比体；
其步骤如下：
1)内插公式为直线内插公式，即：若A(i1,b1),B(i2,b2)为两点，则点P(i,b)在上述两点确定的直线上。则(b-b1)/(i-i1)＝ (b2-b1)/(i2-i1)＝直线斜率，变换即得所求。
2)外推公式为拟合公式，即：将已有样点代入y＝ax
第八步：比例到多井约束纵波速度体；
使用上面百分比体percent2比例纵波网格层析最终速度场到多井约束纵波速度体。
使用公式为vpw＝vp*percent2，其中，vpw为多井约束纵波速度体，vp为纵波网格层析最终速度场,percent2为百分比系数。
第九步：将横波多井百分比体和多井约束纵波速度体结合得到多井约束横波速度体；
使用公式为：vsw＝vpw*percent1，其中，vsw为多井约束横波速度体，vpw为井约束纵波速度体,percent1为第三步所求百分比体。如图5所示，其为多井约束横波速度体示意展示；
本发明经过以上各步处理之后，就可以实现了转换波深度偏移井约束横波初始深度域速度场的求取处理。为转换波深度偏移网格层析反演横波精确速度，奠定了深度域横波速度基础。之后，再进行网格层析反演由浅至深迭代，就可以得到精确的转换波深度偏移所需的深度域横波速度场。
上述未作说明的技术为现有技术，故不再赘述。
可以理解，本发明是通过一些实施例进行描述的，本领域技术人员知悉的，在不脱离本发明的精神和范围的情况下，可以对这些特征和实施例进行各种改变或等效替换。另外，在本发明的教导下，可以对这些特征和实施例进行修改以适应具体的情况及材料而不会脱离本发明的精神和范围。因此，本发明不受此处所公开的具体实施例的限制，所有落入本申请的权利要求范围内的实施例都属于本发明所保护的范围内。
机译： 确定模式转换数据深度迁移的横波速度模型的方法
机译： 通过地震，重力和大地电磁数据的同时联合反演建立叠前深度偏移速度模型的方法。
机译： 通过地震，重力和大地电磁数据的同时联合反演建立叠前深度偏移速度模型的方法。
