主权项 |
一种地层元素测井仪器解谱分析确定产额的方法,其特征在于:包括以下步骤:1)确定段区数n及各区段左右边界k<sub>n</sub>和k<sub>n+1</sub>;所述k<sub>n</sub>和k<sub>n+1</sub>分别表示第n个段区的左右边界;n为不大于256的整数2)计算标准俘获谱归一化响应矩阵A响应矩阵A的计算方法如下;设a<sub>ij</sub>是响应矩阵A的(i,j)元,即A的第i行第j列的元素;则<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><msub><mi>a</mi><mi>ij</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><msub><mi>k</mi><mi>n</mi></msub></mrow><msub><mi>k</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></munderover><msub><mi>N</mi><mi>kj</mi></msub><mo>;</mo></mrow>]]></math><img file="FDA0000566783850000011.GIF" wi="451" he="156" /></maths>k<sub>n</sub>和k<sub>n+1</sub>分别表示第n个段区的左右边界,N<sub>kj</sub>是准谱归一化之后的第k道的计数;3)计算待解谱在各段区总计数值c<sub>i</sub>设在实际测量中得到的256道γ能谱的第i道计数为c<sub>i</sub>;则<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msub><mi>c</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>k</mi><mo>=</mo><msub><mi>k</mi><mi>n</mi></msub></mrow><msub><mi>k</mi><mrow><mi>n</mi><mo>+</mo><mn>1</mn></mrow></msub></munderover><msub><mi>N</mi><mi>kM</mi></msub><mo>,</mo><mi>n</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><msup><mi>n</mi><mo>,</mo></msup></mrow>]]></math><img file="FDA0000566783850000012.GIF" wi="819" he="140" /></maths>其中,k<sub>n</sub>和k<sub>n+1</sub>分别表示第n个段区的左右边界,N<sub>kM</sub>是实测地层谱归一化后的第k道计数;4)计算对角加权矩阵W设对角加权阵为W,w<sub>ii</sub>是W的(i,i)元,令权系数w<sub>i</sub>=w<sub>ii</sub>,则<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><msub><mi>w</mi><mi>ii</mi></msub><mo>=</mo><msub><mi>w</mi><mi>i</mi></msub><mo>=</mo><mfrac><mn>1</mn><msubsup><mi>c</mi><mi>i</mi><mn>2</mn></msubsup></mfrac></mrow>]]></math><img file="FDA0000566783850000013.GIF" wi="282" he="138" /></maths>则:<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><mi>W</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><msub><mi>w</mi><mn>11</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><msub><mi>w</mi><mn>22</mn></msub></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mo>·</mo><mo>·</mo><mo>·</mo></mtd><mtd><mn>0</mn></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><mn>0</mn></mtd><mtd><msub><mi>w</mi><mi>nn</mi></msub></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000566783850000014.GIF" wi="547" he="312" /></maths>5)计算加权最小二乘法解谱得产额y<sub>j</sub>地层测量谱的第i道的计数可以表示为,地层中所有元素在第i道内产生的γ射线即标准谱的线性叠加;<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><msub><mi>c</mi><mi>i</mi></msub><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mo>+</mo><msub><mi>ϵ</mi><mi>i</mi></msub><mo>,</mo><mi>i</mi><mo>=</mo><mn>1,2,3</mn><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>.</mo><msup><mi>m</mi><mo>,</mo></msup></mrow>]]></math><img file="FDA0000566783850000021.GIF" wi="957" he="131" /></maths>即叠加的系数为该地层中元素的产额;其中a<sub>ij</sub>是由元素的标准谱得到的m×m,的响应矩阵A的(i,j)元,y<sub>j</sub>表示第j种元素的百分含量,即产额,ε<sub>i</sub>表示误差,m为元素总数,m,为总道数,采用加权最小二乘法求解,得:y=[A<sup>T</sup>WA]<sup>‑1</sup>[A<sup>T</sup>WC]其中C=(c<sub>1</sub>,c<sub>2</sub>,....c<sub>n</sub>)<sup>T</sup>,y=(y<sub>1</sub>,y<sub>2</sub>,...y<sub>m</sub>);6)利用残差平方和的平均数Δ<sup>2</sup>判定拟合效果;设残差平方和的平均数为<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><msup><mi>Δ</mi><mn>2</mn></msup><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>(</mo><msub><mi>c</mi><mi>i</mi></msub><mo>-</mo><munderover><mi>Σ</mi><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>a</mi><mi>ij</mi></msub><msub><mi>y</mi><mi>j</mi></msub><mo>)</mo></mrow><mn>2</mn></msup></mrow>]]></math><img file="FDA0000566783850000022.GIF" wi="526" he="149" /></maths>若确定的精度Δ<sup>2</sup><ε,则满足要求,否则返回步骤1)重新进行解谱。 |