主权项 |
一种抑制稠油热采井温度监测噪声的方法,其特征在于,利用光纤分布式温度传感器实时采集稠油热采井下温度,得到原始随空间‑时间域温度分布信息;在空间域,采用小波多分辨率分析滤除原始空间‑时间域温度分布信息中的高频噪声,得到初步降噪后的空间‑时间域温度分布信息;在时间域,采用移动灰色模型去除该初步降噪后空间‑时间域温度分布信息的非平稳噪声和干扰,并采用最小二乘法递归优化移动灰色模型中各个窗口的灰色模型参数,得到最终降噪后的空间‑时间域温度分布信息;得到的最终降噪后的空间‑时间域温度分布信息,用于恢复稠油热采井内的地质参数及油藏流体信息,指导采油操作;所述在时间域,采用移动灰色模型去除该初步降噪后空间‑时间域温度分布信息的非平稳噪声和干扰,包括:步骤A:设时间域温度信息原始序列的总点数为N,在时间域温度信息原始序列中从起始位置选取长度为W的信息点应用灰色模型GM(1,1)得到预测序列第1个位置的温度信息,其中W和N均为自然数,且W<N;步骤B:将时间域温度信息原始序列中的移动窗口W向后平移1个信息点,利用灰色模型GM(1,1)得到预测序列第2个位置的温度信息;步骤C:重复步骤B,直到预测序列的长度达到N‑W;步骤D:将时间域温度信息原始序列的前W点和预测序列合并组成处理后的时间域温度信息序列;所述灰色模型GM(1,1)是灰色模型GM(n,m)中的一种模型,灰色模型GM(n,m)是一种用于建模灰色系统动态特性的预测方法,由m个变量的n阶差分方程标识;其中,灰色模型GM(1,1)可被表示为:<maths num="0001" id="cmaths0001"><math><![CDATA[<mrow><mfrac><msup><mi>dX</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mi>dt</mi></mfrac><mo>+</mo><msup><mi>aX</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mo>=</mo><mi>b</mi></mrow>]]></math><img file="FDA0000685234100000021.GIF" wi="300" he="121" /></maths> 公式1其中X<sup>(1)</sup>的分量为:<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><munderover><mi>Σ</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>k</mi></munderover><msup><mi>x</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>,</mo><mi>k</mi><mo>=</mo><mn>1,2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>N</mi></mrow>]]></math><img file="FDA0000685234100000022.GIF" wi="718" he="120" /></maths> 公式2x<sup>(0)</sup>=(x<sup>(0)</sup>(1),x<sup>(0)</sup>(2),...,x<sup>(0)</sup>(N))是时间域温度信息原始序列,单位为摄氏度;N是时间域温度信息原始序列总点数,单位为个;a和b是待定系数,a为规则化系数,无单位,b的单位为摄氏度;所述采用最小二乘法递归优化移动灰色模型中各个窗口的灰色模型参数得到最终降噪后的空间‑时间域温度分布信息的步骤中,所述最小二乘法用于求解灰色模型GM(1,1)的待定系数向量U=(a,b)<sup>T</sup>,根据最小二乘理论,U可被表示为U=(B<sup>T</sup>B)<sup>‑1</sup>B<sup>T</sup>Y<sub>l</sub> 公式3其中<maths num="0003" id="cmaths0003"><math><![CDATA[<mrow><mi>B</mi><mo>=</mo><mfenced open='[' close=']'><mtable><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow><mo>+</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>3</mn><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr><mtr><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd><mtd><mo>.</mo><mo>.</mo><mo>.</mo></mtd></mtr><mtr><mtd><mo>-</mo><mfrac><mn>1</mn><mn>2</mn></mfrac><mrow><mo>(</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>N</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>+</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>N</mi><mo>)</mo></mrow><mo>)</mo></mrow></mtd><mtd><mn>1</mn></mtd></mtr></mtable></mfenced></mrow>]]></math><img file="FDA0000685234100000031.GIF" wi="629" he="407" /></maths> 公式4向量Y<sub>l</sub>(单位为摄氏度)是Y<sub>l</sub><sup>=</sup>(x<sup>(0)</sup>(2),x<sup>(0)</sup>(3),...,x<sup>(0)</sup>(N)) 公式5因此,利用公式2、公式3、公式4和公式5能够计算出待定系数向量U=(a,b)<sup>T</sup>。 |