主权项 |
一种随钻三维反射声波成像测井相控阵声源装置,其特征在于,包括:计算单元,用于计算令井孔中的声源装置在井旁地层中产生的纵波场的主瓣指向目标方向且覆盖所述目标方向上的目标点时所需的控制参数;其中,所述声源装置包括沿钻铤的轴向依次设置的多个相控圆弧阵;所述相控圆弧阵包括按圆周均匀装设于钻铤外表面的多个阵元,所述阵元呈圆弧片状;所述控制参数包括:每个所述相控圆弧阵中各个阵元对应的激励信号的幅度权重,为每个所述相控圆弧阵中相邻阵元施加激励信号的延迟时间,以及为所述声源装置中沿钻铤的轴向的相邻阵元施加激励信号的延迟时间;所述计算单元包括:遍历模块、位置参数分解模块、单层计算模块、邻层计算模块;所述遍历模块,用于遍历每个所述相控圆弧阵,并触发所述位置参数分解模块和所述单层计算模块对当前遍历的所述相控圆弧阵进行计算;所述位置参数分解模块,用于针对当前遍历的所述相控圆弧阵,计算如下各项位置参数:所述目标点在当前遍历的所述相控圆弧阵的中间圆形横截面上的水平投影位置r'和在钻铤的轴向上的垂直投影位置z';所述水平投影位置r'与所述中间圆形横截面的圆心之间的距离r;所述水平投影位置r'在所述中间圆形横截面上所经过的半径与所述中间圆形横截面的基准半径之间的圆心角θ;所述垂直投影位置z'与所述中间圆形横截面的圆心之间的距离z;所述单层计算模块,用于按照如下公式构建关于当前遍历的所述相控圆弧阵中各个阵元对应的激励信号的幅度权重A<sub>1</sub>...A<sub>N</sub>...A<sub>m</sub>和为当前遍历的所述相控圆弧阵中相邻阵元施加激励信号的延迟时间ΔT<sub>1</sub>...ΔT<sub>N</sub>...ΔT<sub>m</sub>的目标函数O(A<sub>1</sub>...A<sub>N</sub>...A<sub>m</sub>,ΔT<sub>1</sub>...ΔT<sub>N</sub>...ΔT<sub>m</sub>):<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mi>O</mi><mrow><mo>(</mo><msub><mi>A</mi><mn>1</mn></msub><mn>...</mn><msub><mi>A</mi><mi>N</mi></msub><mn>...</mn><msub><mi>A</mi><mi>m</mi></msub><mo>,</mo><msub><mi>ΔT</mi><mn>1</mn></msub><mn>...</mn><msub><mi>ΔT</mi><mi>N</mi></msub><mn>...</mn><msub><mi>ΔT</mi><mi>m</mi></msub><mo>)</mo></mrow><mo>=</mo><mfrac><mrow><mi>max</mi><mo>[</mo><msup><mi>Φ</mi><mo>′</mo></msup><mrow><mo>(</mo><mrow><mi>r</mi><mo>,</mo><mi>θ</mi><mo>=</mo><mfrac><mi>π</mi><mn>2</mn></mfrac><mo>~</mo><mfrac><mrow><mn>3</mn><mi>π</mi></mrow><mn>2</mn></mfrac><mo>,</mo><mi>z</mi><mo>,</mo><mi>ω</mi></mrow><mo>)</mo></mrow><mo>]</mo></mrow><mrow><mi>max</mi><mo>[</mo><msup><mi>Φ</mi><mo>′</mo></msup><mrow><mo>(</mo><mrow><mi>r</mi><mo>,</mo><mi>θ</mi><mo>=</mo><mrow><mo>(</mo><mrow><mn>0</mn><mo>~</mo><mfrac><mi>π</mi><mn>2</mn></mfrac><mo>,</mo><mfrac><mrow><mn>3</mn><mi>π</mi></mrow><mn>2</mn></mfrac><mo>~</mo><mn>2</mn><mi>π</mi></mrow><mo>)</mo></mrow><mo>,</mo><mi>z</mi><mo>,</mo><mi>ω</mi></mrow><mo>)</mo></mrow><mo>]</mo></mrow></mfrac></mrow>]]></math><img file="FDA0000962025080000011.GIF" wi="1516" he="295" /></maths>其中,Φ'(r,θ,z,ω)表示当前遍历的所述相控圆弧阵在所述目标点形成的纵波场以位移形式表示的激发谱;且<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>Φ</mi><mo>′</mo></msup><mrow><mo>(</mo><mi>r</mi><mo>,</mo><mi>θ</mi><mo>,</mo><mi>z</mi><mo>,</mo><mi>ω</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>Σ</mo><mrow><mi>n</mi><mo>=</mo><mn>0</mn></mrow><mrow><mo>+</mo><mi>∞</mi></mrow></munderover><munderover><mo>Σ</mo><mrow><mi>N</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>A</mi><mi>N</mi></msub><mo>[</mo><mfrac><mn>1</mn><mrow><mn>2</mn><mi>π</mi></mrow></mfrac><munderover><mo>∫</mo><mrow><mo>-</mo><mi>∞</mi></mrow><mrow><mo>+</mo><mi>∞</mi></mrow></munderover><msubsup><mi>B</mi><mi>n</mi><mrow><mo>(</mo><mi>f</mi><mi>o</mi><mo>)</mo></mrow></msubsup><msub><mi>K</mi><mi>n</mi></msub><mrow><mo>(</mo><mi>q</mi><mi>r</mi><mo>)</mo></mrow><msup><mi>e</mi><mrow><msub><mi>ik</mi><mi>z</mi></msub><mi>z</mi></mrow></msup><msub><mi>dk</mi><mi>z</mi></msub><mo>]</mo><mo>×</mo><mi>exp</mi><mrow><mo>(</mo><mo>-</mo><msub><mi>iωΔT</mi><mi>N</mi></msub><mo>)</mo></mrow><mi>cos</mi><mo>[</mo><mrow><mi>n</mi><mrow><mo>(</mo><mrow><mi>θ</mi><mo>-</mo><mfrac><mrow><mi>N</mi><mi>π</mi><mo>-</mo><mi>π</mi></mrow><mn>4</mn></mfrac></mrow><mo>)</mo></mrow></mrow><mo>]</mo></mrow>]]></math><img file="FDA0000962025080000021.GIF" wi="1789" he="158" /></maths>ω表示所述相控圆弧阵中各个阵元的振动频率;n表示方位阶数;N表示阵元的序号;m表示所述相控圆弧阵中阵元的数量;A<sub>N</sub>表示序号是N的阵元对应的激励信号的幅度权重;ΔT<sub>N</sub>表示为序号是N+1的阵元施加激励信号相比于为序号是N的阵元施加激励信号的延迟时间;<img file="FDA0000962025080000022.GIF" wi="94" he="78" />表示地层中向外传播的纵波的势函数的系数;第二类贝塞尔函数K<sub>n</sub>(qr)与<img file="FDA0000962025080000023.GIF" wi="94" he="71" />共同描述了沿径向由内向外传播的波;q表示纵波的径向波数;i表示虚数单位;k<sub>z</sub>表示液体内纵波的轴向波数;所述单层计算模块,还用于采用非常快速模拟退火算法,计算令所述目标函数O(A<sub>1</sub>...A<sub>N</sub>...A<sub>m</sub>,ΔT<sub>1</sub>...ΔT<sub>N</sub>...ΔT<sub>m</sub>)取极小值时,当前遍历的所述相控圆弧阵中各个所述阵元对应的激励信号的幅度权重A<sub>1</sub>...A<sub>N</sub>...A<sub>m</sub>和为当前遍历的所述相控圆弧阵中相邻阵元施加激励信号的延迟时间ΔT<sub>1</sub>...ΔT<sub>N</sub>...ΔT<sub>m</sub>;所述邻层计算模块,用于计算当井孔中的所述声源装置在井旁地层中产生的纵波场的主瓣指向所述目标方向且覆盖所述目标点时,所述声源装置的主声束偏转角度θ<sub>0</sub>;并按照如下公式计算为所述声源装置中沿钻铤的轴向的相邻阵元施加激励信号的延迟时间τ:<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>τ</mi><mo>=</mo><mfrac><mrow><mi>h</mi><mo>×</mo><msub><mi>sinθ</mi><mn>0</mn></msub></mrow><msub><mi>v</mi><mi>f</mi></msub></mfrac></mrow>]]></math><img file="FDA0000962025080000024.GIF" wi="270" he="142" /></maths>其中,h表示相邻的两个相控圆弧阵的中间圆形横截面之间的距离,v<sub>f</sub>表示井内泥浆的声速;按照所述计算单元针对每个所述相控圆弧阵计算出的幅度权重A<sub>1</sub>...A<sub>N</sub>...A<sub>m</sub>和延迟时间ΔT<sub>1</sub>...ΔT<sub>N</sub>...ΔT<sub>m</sub>激励相应所述相控圆弧阵中的各个阵元振动,并按照所述计算单元计算出的延迟时间τ为所述声源装置中沿钻铤的轴向的相邻阵元施加激励信号。 |