发明名称 一种基于灰色系统理论的提高采收率潜力预测方法
摘要 一种基于灰色系统理论的提高采收率潜力预测方法,包括以下步骤:A.根据油田实际情况,选取连续时间段的提高采收率数据,并将数据分为模拟数据集和检测数据集;B.利用模拟数据集构建原始序列,获得1-AGO序列;C.利用优化的模糊预测模型背景值公式计算1-AGO序列的紧邻均值生成序列;D.计算发展系数和灰色作用量,求得时间响应式和还原值;E.利用步骤D得到模型计算模拟值做拟合精度检验;F.利用步骤D所得到的还原值计算方法进行提高采收率潜力预测。该方法通过对提高采收率数据分析,消除因影响因素选取不当造成的预测误差。同时,对模型模糊预测模型进行优化改进,模型预测精度更高,使提高采收率预测果更加准确、公正、合理。
申请公布号 CN104715298A 申请公布日期 2015.06.17
申请号 CN201510169962.4 申请日期 2015.04.10
申请人 中国石油大学(华东) 发明人 李克文;王义龙;苏玉亮
分类号 G06Q10/04(2012.01)I;G06Q50/02(2012.01)I 主分类号 G06Q10/04(2012.01)I
代理机构 代理人
主权项 一种基于灰色系统理论的提高采收率潜力预测方法其特征在于,主要包括以下步骤:A.获取数据根据油田实际生产情况,选取连续时间段的原油提高采收率数据,并将数据分为模拟数据集和检测数据集,分别用于提高采收率预测模型的数据模拟和对模型运算精度检验;B.构建预测模型(1)利用模拟数据集,构建原始序列:X<sup>(0)</sup>={x<sup>(0)</sup>(1),x<sup>(0)</sup>(2),…,x<sup>(0)</sup>(n)}其中,x<sup>(0)</sup>(k)≥0,k=1,2,…,n;(2)建立X<sup>(0)</sup>的1‑AGO序列X<sup>(1)</sup>:X<sup>(1)</sup>={x<sup>(1)</sup>(1),x<sup>(1)</sup>(2),…,x<sup>(1)</sup>(n)}其中,<maths num="0001" id="cmaths0001"><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>&Sigma;</mi><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mi>x</mi><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>i</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><mo>,</mo></mrow>]]></math><img file="FDA0000696977770000011.GIF" wi="898" he="157" /></maths>并将x<sup>(1)</sup>(n)作为初始条件值;(3)优化模糊预测模型背景值公式,计算X<sup>(1)</sup>的紧邻均值生成序列Z<sup>(1)</sup>:Z<sup>(1)</sup>={z<sup>(1)</sup>(1),z<sup>(1)</sup>(2),…,z<sup>(1)</sup>(n)},其中,<maths num="0002" id="cmaths0002"><math><![CDATA[<mrow><msup><mi>z</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mfrac><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><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><mi>ln</mi><mo>[</mo><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><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><mo>-</mo><mi>ln</mi><mo>[</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow><mo>]</mo></mrow></mfrac><mo>-</mo><mfrac><mrow><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow><mo>&CenterDot;</mo><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo>)</mo></mrow></mrow><mrow><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><msup><mi>x</mi><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow></mfrac><mo>,</mo><mi>k</mi><mo>=</mo><mn>2,3</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi></mrow>]]></math><img file="FDA0000696977770000012.GIF" wi="1805" he="210" /></maths>(4)计算时间响应函数灰微分方程方程x<sup>(1)</sup>(k)+az<sup>(1)</sup>(k)=b的时间响应函数为:<maths num="0003" id="cmaths0003"><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><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><mfrac><mi>b</mi><mi>a</mi></mfrac><mo>]</mo><msup><mi>e</mi><mrow><mo>-</mo><mi>a</mi><mrow><mo>(</mo><mi>k</mi><mo>-</mo><mi>n</mi><mo>)</mo></mrow></mrow></msup><mo>+</mo><mfrac><mi>b</mi><mi>a</mi></mfrac><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="FDA0000696977770000013.GIF" wi="1121" he="190" /></maths>其中,‑a为发展系数,b为灰色作用量;(5)计算还原值<maths num="0004" id="cmaths0004"><math><![CDATA[<mrow><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>=</mo><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></msup><mrow><mo>(</mo><mi>k</mi><mo>+</mo><mn>1</mn><mo>)</mo></mrow><mo>-</mo><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>1</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><mo>,</mo></mrow>]]></math><img file="FDA0000696977770000021.GIF" wi="1243" he="136" /></maths>是模型的时间响应函数的累减还原值,当k≤n时,称<img file="FDA0000696977770000022.GIF" wi="167" he="105" />为预测模型模拟值;当k&gt;n时,称<img file="FDA0000696977770000023.GIF" wi="169" he="107" />为预测模型预测值;C.拟合精度检验分别对模型进行残差检验和后验差检验,残差检验分别计算:残差:<maths num="0005" id="cmaths0005"><math><![CDATA[<mrow><mi>E</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><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><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>,</mo><mi>k</mi><mo>=</mo><mn>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo></mrow>]]></math><img file="FDA0000696977770000024.GIF" wi="914" he="126" /></maths>相对残差:<maths num="0006" id="cmaths0006"><math><![CDATA[<mrow><mi>e</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>=</mo><mo>[</mo><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><msup><mover><mi>x</mi><mo>^</mo></mover><mrow><mo>(</mo><mn>0</mn><mo>)</mo></mrow></msup><mo>]</mo><mo>/</mo><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>2</mn><mo>,</mo><mo>.</mo><mo>.</mo><mo>.</mo><mo>,</mo><mi>n</mi><mo>;</mo></mrow>]]></math><img file="FDA0000696977770000025.GIF" wi="1162" he="135" /></maths>后验差检验分别计算:x<sup>(0)</sup>的均值:<maths num="0007" id="cmaths0007"><math><![CDATA[<mrow><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</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></mrow>]]></math><img file="FDA0000696977770000026.GIF" wi="448" he="145" /></maths>x<sup>(0)</sup>的方差:<maths num="0008" id="cmaths0008"><math><![CDATA[<mrow><msub><mi>S</mi><mn>1</mn></msub><mo>=</mo><msqrt><mfrac><mn>1</mn><mi>n</mi></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>1</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>[</mo><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><mover><mi>X</mi><mo>&OverBar;</mo></mover><mo>]</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000696977770000027.GIF" wi="650" he="170" /></maths>残差的均值:<maths num="0009" id="cmaths0009"><math><![CDATA[<mrow><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>2</mn></mrow><mi>n</mi></munderover><mi>E</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>;</mo></mrow>]]></math><img file="FDA0000696977770000028.GIF" wi="467" he="148" /></maths>残差的方差:<maths num="0010" id="cmaths0010"><math><![CDATA[<mrow><msub><mi>S</mi><mn>2</mn></msub><mo>=</mo><msqrt><mfrac><mn>1</mn><mrow><mi>n</mi><mo>-</mo><mn>1</mn></mrow></mfrac><munderover><mi>&Sigma;</mi><mrow><mi>k</mi><mo>=</mo><mn>2</mn></mrow><mi>n</mi></munderover><msup><mrow><mo>[</mo><mi>E</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>]</mo></mrow><mn>2</mn></msup></msqrt><mo>;</mo></mrow>]]></math><img file="FDA0000696977770000029.GIF" wi="683" he="170" /></maths>后验差比值:<img file="FDA00006969777700000210.GIF" wi="224" he="156" />小误差概率:<maths num="0011" id="cmaths0011"><math><![CDATA[<mrow><mi>P</mi><mo>=</mo><mi>PE</mi><mo>{</mo><mo>|</mo><mi>E</mi><mrow><mo>(</mo><mi>k</mi><mo>)</mo></mrow><mo>-</mo><mover><mi>E</mi><mo>&OverBar;</mo></mover><mo>|</mo><mo>&lt;</mo><mn>0.6745</mn><msub><mi>S</mi><mn>1</mn></msub><mo>}</mo></mrow>]]></math><img file="FDA00006969777700000211.GIF" wi="771" he="110" /></maths>D.提高采收率预测利用预测模型,通过控制参数k与n的大小,选择要预测的时间点,模型输出即为油田原油提高采收率预测值。
地址 266580 山东省青岛市青岛经济技术开发区长江西路66号