主权项 |
一种确定煤层气纵波和横波速度的方法,特点是通过以下步骤实现:1)采集地震勘探和测井数据,处理后得到纵波速度x和横波速度y;2)按照下式分别对纵波速度x和横波速度y进行处理,使它们的平均值分别为零:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfenced open='{' close=''><mtable><mtr><mtd><msub><mi>x</mi><mi>j</mi></msub><mo>=</mo><msub><mi>x</mi><mi>j</mi></msub><mo>-</mo><mrow><mo>(</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>/</mo><mi>L</mi></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>j</mi></msub><mo>=</mo><msub><mi>y</mi><mi>j</mi></msub><mo>-</mo><mrow><mo>(</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>/</mo><mi>L</mi></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000388571180000011.GIF" wi="1430" he="342" /></maths>式中:L是组数,j是纵波速度x和横波速度y的组序列号,x<sub>j</sub>和y<sub>j</sub>分别是纵波速度x和横波速度y的第j个元素;3)按照下式对L组纵波速度x和横波速度y进行归一化处理:<maths num="0002" id="cmaths0002"><math><![CDATA[<mfenced open='{' close='' separators=''><mtable><mtr><mtd><msub><mi>x</mi><mi>j</mi></msub><mo>=</mo><msub><mi>x</mi><mi>j</mi></msub><mo>/</mo><mi>max</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mtd><mtd></mtd></mtr><mtr><mtd><msub><mi>y</mi><mi>j</mi></msub><mo>=</mo><msub><mi>y</mi><mi>j</mi></msub><mo>/</mo><mi>max</mi><mrow><mo>(</mo><mi>y</mi><mo>)</mo></mrow></mtd><mtd><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>L</mi><mo>)</mo></mrow></mtd></mtr></mtable><mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow></mfenced>]]></math><img file="FDA0000388571180000012.GIF" wi="1389" he="214" /></maths>式中,L是组数;j是纵波速度x和横波速度y的组序列号;x<sub>j</sub>和y<sub>j</sub>分别是是纵波速度x和横波速度y的第j个元素;max(x)是求取纵波速度x的最大值;max(y)是求取横波速度y的最大值;4)按下式分别计算纵波速度x和横波速度y最佳变换的初始值φ<sup>0</sup>(x)和θ<sup>0</sup>(y):<maths num="0003" id="cmaths0003"><math><![CDATA[<mfenced open='{' close='' separators=''><mtable><mtr><mtd><msup><mi>φ</mi><mn>0</mn></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><mn>0</mn></mtd><mtd></mtd></mtr><mtr><mtd><msup><mi>θ</mi><mn>0</mn></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>y</mi><mi>j</mi></msub></mtd><mtd><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>L</mi><mo>)</mo></mrow></mtd></mtr></mtable><mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow></mrow></mfenced>]]></math><img file="FDA0000388571180000013.GIF" wi="1448" he="217" /></maths>式中,L是纵波速度x和横波速度y的组数;j是纵波速度x和横波速度y的组序列号;x<sub>j</sub>和y<sub>j</sub>分别是纵波速度x和横波速度y的第j个元素;5)按照下式求取纵波速度x和横波速度y的最佳变换φ<sup>k</sup>(x<sub>j</sub>)和θ<sup>k</sup>(y<sub>j</sub>):<maths num="0004" id="cmaths0004"><math><![CDATA[<mfenced open='{' close='' separators=''><mtable><mtr><mtd><msup><mi>φ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>x</mi></msub><mrow><mo>(</mo><msup><mi>θ</mi><mrow><mi>k</mi><mo>-</mo><mn>1</mn></mrow></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd></mtd></mtr><mtr><mtd><msup><mi>θ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>=</mo><msub><mi>S</mi><mi>y</mi></msub><mo>[</mo><msup><mi>φ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>/</mo><mo>|</mo><mo>|</mo><msub><mi>S</mi><mi>y</mi></msub><mo>[</mo><msup><mi>φ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>]</mo><mo>|</mo><mo>|</mo></mtd><mtd><mrow><mo>(</mo><mi>j</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>·</mo><mo>·</mo><mo>·</mo><mo>,</mo><mi>L</mi><mo>)</mo></mrow></mtd></mtr></mtable><mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow></mfenced>]]></math><img file="FDA0000388571180000021.GIF" wi="1437" he="223" /></maths>式中,k是迭代次数;S是指对数据作平滑,平滑方法采用样条函数法或超级平滑方法;x<sub>j</sub>和y<sub>j</sub>分别是纵波速度x和横波速度y的第j个元素;φ<sup>k</sup>(x<sub>j</sub>)和θ<sup>k</sup>(y<sub>j</sub>)分别是第k次迭代,纵波速度x和横波速度y最佳变换φ<sup>k</sup>(x)和θ<sup>k</sup>(y)的第j个元素;6)按下式计算迭代误差:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mrow><msup><mi>e</mi><mn>2</mn></msup><mo>[</mo><msup><mi>θ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>,</mo><msup><mi>φ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mo>=</mo><mfrac><mn>1</mn><mi>L</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>L</mi></munderover><msup><mrow><mo>[</mo><msup><mi>θ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>-</mo><msup><mi>φ</mi><mi>k</mi></msup><mrow><mo>(</mo><msub><mi>x</mi><mi>j</mi></msub><mo>)</mo></mrow><mo>]</mo></mrow><mn>2</mn></msup><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>]]></math><img file="FDA0000388571180000022.GIF" wi="1635" he="147" /></maths>式中:k是迭代次数;7)重复步骤5)和6),直到相邻两次迭代误差满足下式:|e<sup>2</sup>[θ<sup>k</sup>(y<sub>j</sub>),φ<sup>k</sup>(x<sub>j</sub>)]‑e<sup>2</sup>[θ<sup>k‑1</sup>(y<sub>j</sub>),φ<sup>k‑1</sup>(x<sub>j</sub>)]|≤δ (6)式中,k是迭代次数;δ是误差门槛值,取值范围是0.001到0.01;当相邻两次迭代误差满足公式(6)时,第k次迭代得到θ(y)称为横波速度y的最佳变换函数,得到的φ(x)称为纵波速度x的最佳变换函数;φ<sup>k</sup>(x<sub>j</sub>)和θ<sup>k</sup>(y<sub>j</sub>)分别是第k次迭代,纵波速度x和横波速度y最佳变换φ<sup>k</sup>(x)和θ<sup>k</sup>(y)的第j个元素;8)将步骤7)得到的最佳变换函数θ(y)和φ(x)的函数值,利用线性回归方法进行回归,得到θ(y)和φ(x)的函数关系式:θ(y)=f(φ(x)) (7)式中,f是由线性回归得方法到的θ(y)和φ(x)的函数关系式;9)根据函数关系式θ(y),φ(x)和f,计算得到其他井或地区的纵波速度x或者横波速度y。 |