发明名称 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置
摘要 本发明提供了一种基于旅行时法确定地层裂隙裂缝走向的方法及装置,方法包括:获取三维动校正前的共中心点道集数据、共中心点坐标及纵测线方位角;以共中心点坐标为圆心,纵测线方位角的方位为起始方位,逆时针以预设角度为间隔布设测线;根据方位测线从共中心点道集数据中获取方位共中心点道集数据;根据方位共中心点道集数据选取目标层,确定目标层的多个偏移距方位旅行时;利用统计和加权平均方法,通过对多对方位测线及每对方位测线上不同偏移距方位旅行时差进行最小二乘法线性拟合和加权平均,获得裂隙裂缝的走向。本发明降低现有方位旅行时差法预测裂隙裂缝方向时误差的随机性和不确定性,提高裂隙裂缝预测结果的精度和稳定性。
申请公布号 CN103954996B 申请公布日期 2016.10.26
申请号 CN201410128443.9 申请日期 2014.04.01
申请人 中国石油天然气股份有限公司 发明人 石玉梅;谢韬;马晓宇;李勇根;秦臻
分类号 G01V1/30(2006.01)I 主分类号 G01V1/30(2006.01)I
代理机构 北京三友知识产权代理有限公司 11127 代理人 任默闻
主权项 一种基于旅行时法确定地层裂隙裂缝走向的方法,其特征在于,所述的方法包括:步骤1,获取三维动校正前的共中心点道集数据、共中心点坐标及纵测线方位角,其中,所述纵测线方位角为纵测线与正北方向夹角;步骤2,以所述共中心点坐标为圆心,所述纵测线方位角的方位为起始方位,逆时针以预设角度为间隔布设测线,获得L条方位测线;步骤3,根据所述的方位测线从所述共中心点道集数据中获取方位共中心点道集数据;步骤4,根据所述方位共中心点道集数据选取目标层,确定所述L条方位测线的目标层的多个偏移距方位旅行时;步骤5,将所述的L条方位测线中相互垂直的测线进行两两组合,生成<img file="FDA0000958636960000011.GIF" wi="46" he="119" />对方位测线;步骤6,根据所述L条方位测线的目标层的旅行时和式(1)确定所述<img file="FDA0000958636960000012.GIF" wi="50" he="114" />对方位测线的每对方位测线的第一方位旅行时差,<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>&Delta;t</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msubsup><mi>t</mi><mn>2</mn><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><msubsup><mi>t</mi><mn>1</mn><mi>i</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000958636960000013.GIF" wi="1062" he="71" /></maths>其中,i表示第i对方位测线,x表示偏移距,<img file="FDA0000958636960000014.GIF" wi="97" he="63" />和<img file="FDA0000958636960000015.GIF" wi="99" he="71" />分别为第i对方位测线中第1条方位测线和第2条方位测线上目标层方位旅行时,测线对内测线序号按逆时针编号,Δt<sub>i</sub>(x)为第i对方位测线的第一方位旅行时差;步骤7,对所述每对方位测线的方位旅行时差进行两两组合,按式(2)和式(3)重新计算每组内方位旅行时差,生成<img file="FDA0000958636960000016.GIF" wi="115" he="111" />组第二方位旅行时差,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msubsup><mi>&Delta;t</mi><mi>i</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><msub><mi>&Delta;t</mi><mi>i</mi></msub><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000958636960000017.GIF" wi="1062" he="71" /></maths><maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msubsup><mi>&Delta;t</mi><mi>j</mi><mi>k</mi></msubsup><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow><mo>=</mo><mrow><mo>(</mo><msub><mi>&Delta;t</mi><mi>j</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mo>-</mo><msub><mi>&Delta;t</mi><mi>i</mi></msub><mo>(</mo><mi>x</mi><mo>)</mo><mi>c</mi><mi>o</mi><mi>s</mi><mn>2</mn><msub><mi>&phi;</mi><mn>0</mn></msub><mo>)</mo></mrow><mo>/</mo><mi>sin</mi><mn>2</mn><msub><mi>&phi;</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000958636960000018.GIF" wi="1078" he="87" /></maths>其中,k表示第k组方位旅行时差,由步骤6所述的第i对方位测线的方位旅行时差和第j对方位测线的方位旅行时差组成,Δt<sub>i</sub>(x)和Δt<sub>j</sub>(x)分别表示步骤6所述的第i对方位测线的第一方位旅行时差和第j对方位测线的第一方位旅行时差,<img file="FDA0000958636960000021.GIF" wi="134" he="71" />和<img file="FDA0000958636960000022.GIF" wi="134" he="71" />分别为重新计算后第k组的第二方位旅行时差,φ<sub>0</sub>为进行两两组合的第i对方位测线与第j对方位测线间的夹角;步骤8,分别以各组中的第二方位旅行时差<img file="FDA0000958636960000023.GIF" wi="133" he="69" />和<img file="FDA0000958636960000024.GIF" wi="136" he="71" />为横、纵坐标轴,对每组内多偏移距第二方位旅行时差进行交会,获得离散交会点;步骤9,对所述的每组离散交会点进行线性拟合,根据式(4)进行计算,获得各组裂隙裂缝方位角,<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>=</mo><mrow><mo>(</mo><mi>i</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mi>&theta;</mi><mo>+</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><msub><mi>&phi;</mi><mi>k</mi></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000958636960000025.GIF" wi="934" he="117" /></maths>其中,α<sub>k</sub>为第k组方位旅行时差获得的裂隙裂缝方位角,i为步骤7所述的第i对正交测线,θ为步骤2中布设测线的预设角度,φ<sub>k</sub>为第k组方位旅行时差交会拟合直线与横坐标轴的夹角;步骤10,根据确定的各组的裂隙裂缝方位角、所述纵测线方位角和式(5)确定地层裂隙裂缝走向,<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>&alpha;</mi><mo>=</mo><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>w</mi><mi>k</mi></msub><msub><mi>&alpha;</mi><mi>k</mi></msub><mo>+</mo><msub><mi>&alpha;</mi><mn>0</mn></msub><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000958636960000026.GIF" wi="712" he="135" /></maths>其中,α为裂隙裂缝与正北方向的夹角,N为步骤7所述的第二方位旅行时差组的组数,α<sub>k</sub>为步骤9所述的第k组方位旅行时差获得的裂隙裂缝方位角,α<sub>0</sub>为所述的纵测线方位角,w<sub>k</sub>为加权系数,由所述的离散交会点与拟合直线偏离程度决定,且<img file="FDA0000958636960000027.GIF" wi="219" he="127" />
地址 100007 北京市东城区东直门北大街9号