页岩气藏数值模拟器输入参数的校正方法及装置与流程

文档序号:12118339阅读:来源:国知局

技术特征:

1.一种页岩气藏数值模拟器输入参数的校正方法,其特征在于,包括:

根据吸附气储气能力和吸附气密度得到岩石物理模型中吸附气体积;

根据岩石物理模型中吸附气体积得到岩石物理模型吸附气孔隙度;

根据岩石物理模型有效孔隙度和岩石物理模型吸附气孔隙度,得到校正的页岩气藏数值模拟器输入参数。

2.根据权利要求1所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,根据岩石物理模型有效孔隙度和岩石物理模型吸附气孔隙度,得到校正的页岩气藏数值模拟器输入参数,包括:

根据岩石物理模型整体孔隙度和岩石物理模型吸附气孔隙度,得到校正的体相孔隙度;

根据校正的体相孔隙度得到校正的束缚水饱和度;

根据校正的束缚水饱和度得到校正的自由气饱和度;

根据岩石物理模型有效孔隙度和校正的体相孔隙度得到不同压力下校正的孔隙度乘子;

根据岩石物理模型有效孔隙度得到校正的渗透率乘子。

3.根据权利要求2所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,根据岩石物理模型整体孔隙度和岩石物理模型吸附气孔隙度,得到校正的体相孔隙度,包括:

根据自由气储气能力和自由气密度得到岩石物理模型中自由气体积;根据岩石物理模型中吸附气体积、岩石物理模型中自由气体积以及岩石物理模型自由气孔隙度,得到岩石物理模型吸附气孔隙度;根据岩石物理模型整体孔隙度和岩石物理模型吸附气孔隙度,得到校正的体相孔隙度;

根据校正的体相孔隙度得到校正的束缚水饱和度,包括:根据校正的体相孔隙度、岩石物理模型整体孔隙度、岩石物理模型束缚水饱和度,得到校正的束缚水饱和度。

4.根据权利要求3所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,

按如下公式得到校正的体相孔隙度:

<mrow> <msubsup> <mi>&phi;</mi> <mi>t</mi> <mo>&prime;</mo> </msubsup> <mo>=</mo> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>-</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>G</mi> <mi>a</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>a</mi> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>

其中,φ′t为校正的体相孔隙度,φt为岩石物理模型整体孔隙度,φgf为岩石物理模型自由气孔隙度,Ga为吸附气储气能力,ρgf为自由气密度,Gf为自由气储气能力,ρga为吸附气密度;

按如下公式得到校正的束缚水饱和度:

<mrow> <msubsup> <mi>S</mi> <mrow> <mi>w</mi> <mi>c</mi> </mrow> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>S</mi> <mrow> <mi>w</mi> <mi>c</mi> </mrow> </msub> </mrow> <msubsup> <mi>&phi;</mi> <mi>t</mi> <mo>&prime;</mo> </msubsup> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>a</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>S</mi> <mrow> <mi>w</mi> <mi>c</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>a</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>-</mo> <msub> <mi>&phi;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>G</mi> <mi>a</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>

其中,S'wc为校正的束缚水饱和度,Swc为岩石物理模型束缚水饱和度;

按如下公式得到校正的自由气饱和度:

<mrow> <msubsup> <mi>S</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mn>1</mn> <mo>-</mo> <msubsup> <mi>S</mi> <mrow> <mi>w</mi> <mi>c</mi> </mrow> <mo>&prime;</mo> </msubsup> <mo>=</mo> <mfrac> <mrow> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>a</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>&CenterDot;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <msub> <mi>S</mi> <mrow> <mi>w</mi> <mi>c</mi> </mrow> </msub> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>&phi;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>G</mi> <mi>a</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> </mrow> <mrow> <msub> <mi>G</mi> <mi>f</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>a</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>&phi;</mi> <mi>t</mi> </msub> <mo>-</mo> <msub> <mi>&phi;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> <mo>&CenterDot;</mo> <msub> <mi>G</mi> <mi>a</mi> </msub> <mo>&CenterDot;</mo> <msub> <mi>&rho;</mi> <mrow> <mi>g</mi> <mi>f</mi> </mrow> </msub> </mrow> </mfrac> <mo>.</mo> </mrow>

5.根据权利要求2所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,根据岩石物理模型有效孔隙度和校正的体相孔隙度得到不同压力下校正的孔隙度乘子,包括:

根据不考虑吸附层的孔隙半径和去除第一吸附层的孔隙半径,得到不同压力下的有效孔隙半径和不同压力下的有效孔隙体积;

根据不同压力下的有效孔隙体积和岩块体积,得到岩石物理模型有效孔隙度;

根据岩石物理模型有效孔隙度和校正的体相孔隙度,得到不同压力下校正的孔隙度乘子。

6.根据权利要求5所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,按如下公式得到不同压力下的有效孔隙半径:

<mrow> <msub> <mi>R</mi> <mi>e</mi> </msub> <mo>=</mo> <mroot> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mo>&Pi;</mo> <mo>)</mo> <msubsup> <mi>R</mi> <mn>0</mn> <mn>3</mn> </msubsup> <mo>+</mo> <mo>&Pi;</mo> <msubsup> <mi>R</mi> <mn>1</mn> <mn>3</mn> </msubsup> </mrow> <mn>3</mn> </mroot> <mo>,</mo> </mrow>

其中,Re为不同压力下的有效孔隙半径,R0为不考虑吸附层的孔隙半径,R1为去除第一吸附层的孔隙半径,C为与净热吸附有关的常数,p为气藏压力,Z为气体压缩因子,p0为气体饱和吸附压力,N为吸附层数;

按如下公式得到不同压力下的有效孔隙体积:

<mrow> <msub> <mi>V</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>4</mn> <mi>&pi;</mi> <mo>&lsqb;</mo> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mo>&Pi;</mo> <mo>)</mo> </mrow> <msubsup> <mi>R</mi> <mn>0</mn> <mn>3</mn> </msubsup> <mo>+</mo> <mo>&Pi;</mo> <msubsup> <mi>R</mi> <mn>1</mn> <mn>3</mn> </msubsup> <mo>&rsqb;</mo> </mrow> <mn>3</mn> </mfrac> <mo>=</mo> <msub> <mi>V</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>r</mi> </msub> <msub> <mi>V</mi> <mi>b</mi> </msub> <msub> <mi>V</mi> <mi>L</mi> </msub> <mo>&Pi;</mo> <mo>,</mo> </mrow>

其中,Ve为不同压力下的有效孔隙体积,V0为不考虑吸附气的孔隙体积,ρr为岩石密度,Vb为岩块体积,VL为朗缪尔体积;

按如下公式得到岩石物理模型有效孔隙度:

<mrow> <msub> <mi>&phi;</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <msub> <mi>V</mi> <mi>e</mi> </msub> <msub> <mi>V</mi> <mi>b</mi> </msub> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>V</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>r</mi> </msub> <msub> <mi>V</mi> <mi>b</mi> </msub> <msub> <mi>V</mi> <mi>L</mi> </msub> <mo>&Pi;</mo> </mrow> <msub> <mi>V</mi> <mi>b</mi> </msub> </mfrac> <mo>=</mo> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>r</mi> </msub> <msub> <mi>V</mi> <mi>L</mi> </msub> <mo>&Pi;</mo> <mo>,</mo> </mrow>

其中,φe为岩石物理模型有效孔隙度,φ0为初始条件固有孔隙度,φr为高有效应力下孔隙度,φ为有效应力为零时孔隙度,η为孔隙度应力敏感参数,pi为初始气藏压力;

按如下公式得到不同压力下校正的孔隙度乘子:

<mrow> <msub> <mi>R</mi> <mi>&phi;</mi> </msub> <mo>=</mo> <mfrac> <msub> <mi>&phi;</mi> <mi>e</mi> </msub> <msubsup> <mi>&phi;</mi> <mi>t</mi> <mo>&prime;</mo> </msubsup> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>r</mi> </msub> <mo>+</mo> <mrow> <mo>(</mo> <msub> <mi>&phi;</mi> <mi>&infin;</mi> </msub> <mo>-</mo> <msub> <mi>&phi;</mi> <mi>r</mi> </msub> <mo>)</mo> </mrow> <msup> <mi>e</mi> <mrow> <mo>-</mo> <mi>&eta;</mi> <mrow> <mo>(</mo> <msub> <mi>p</mi> <mi>i</mi> </msub> <mo>-</mo> <mi>p</mi> <mo>)</mo> </mrow> </mrow> </msup> <mo>-</mo> <msub> <mi>&rho;</mi> <mi>r</mi> </msub> <msub> <mi>V</mi> <mi>L</mi> </msub> <mo>&Pi;</mo> </mrow> <msubsup> <mi>&phi;</mi> <mi>t</mi> <mo>&prime;</mo> </msubsup> </mfrac> <mo>,</mo> </mrow>

其中,Rφ为不同压力下校正的孔隙度乘子,φ′t为校正的体相孔隙度。

7.根据权利要求2所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,根据岩石物理模型有效孔隙度得到校正的渗透率乘子,包括:

根据岩石物理模型有效孔隙度、不同压力下的有效孔隙半径和迂曲度,得到有效渗透率;

根据有效渗透率、克努森函数和表面扩散等效渗透率,得到表观渗透率;

根据初始条件固有孔隙度、不考虑吸附层的孔隙半径、表观渗透率和迂曲度,得到校正的渗透率乘子。

8.根据权利要求7所述的页岩气藏数值模拟器输入参数的校正方法,其特征在于,按如下公式得到有效渗透率:

<mrow> <msub> <mi>K</mi> <mi>e</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>e</mi> </msub> <msubsup> <mi>R</mi> <mi>e</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <mn>8</mn> <msub> <mi>&tau;</mi> <mi>m</mi> </msub> </mrow> </mfrac> <mo>,</mo> </mrow>

其中,Ke为有效渗透率,φe为岩石物理模型有效孔隙度,τm为迂曲度,Re为不同压力下的有效孔隙半径,如下:

<mrow> <msub> <mi>R</mi> <mi>e</mi> </msub> <mo>=</mo> <mroot> <mrow> <mo>(</mo> <mn>1</mn> <mo>-</mo> <mo>&Pi;</mo> <mo>)</mo> <msubsup> <mi>R</mi> <mn>0</mn> <mn>3</mn> </msubsup> <mo>+</mo> <mo>&Pi;</mo> <msubsup> <mi>R</mi> <mn>1</mn> <mn>3</mn> </msubsup> </mrow> <mn>3</mn> </mroot> <mo>,</mo> </mrow>

其中,Re为不同压力下的有效孔隙半径,R0为不考虑吸附层的孔隙半径,如下:

<mrow> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>=</mo> <mn>2</mn> <msqrt> <mrow> <mn>2</mn> <msub> <mi>&tau;</mi> <mi>m</mi> </msub> </mrow> </msqrt> <msqrt> <mfrac> <msub> <mi>K</mi> <mrow> <mi>s</mi> <mi>d</mi> </mrow> </msub> <mi>&phi;</mi> </mfrac> </msqrt> <mo>,</mo> </mrow>

其中,Ksd为具有盈利敏感效应的固有渗透率,φ为不同压力下孔隙度;

R1为去除第一吸附层的孔隙半径,如下:

<mrow> <msub> <mi>R</mi> <mn>1</mn> </msub> <mo>=</mo> <msub> <mi>R</mi> <mn>0</mn> </msub> <mo>-</mo> <msub> <mi>&Gamma;d</mi> <mrow> <msub> <mi>CH</mi> <mn>4</mn> </msub> </mrow> </msub> <mo>,</mo> </mrow>

其中,Γ为吸附层校正系数,为甲烷分子半径;

<mrow> <mo>&Pi;</mo> <mo>=</mo> <mfrac> <mrow> <mi>C</mi> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> </mrow> </mfrac> <mo>&lsqb;</mo> <mfrac> <mrow> <mn>1</mn> <mo>-</mo> <mrow> <mo>(</mo> <mi>N</mi> <mo>+</mo> <mn>1</mn> <mo>)</mo> </mrow> <msup> <mrow> <mo>(</mo> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mi>N</mi> </msup> <mo>+</mo> <mi>N</mi> <msup> <mrow> <mo>(</mo> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> </mrow> <mrow> <mn>1</mn> <mo>+</mo> <mrow> <mo>(</mo> <mi>C</mi> <mo>-</mo> <mn>1</mn> <mo>)</mo> </mrow> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> <mo>-</mo> <mi>C</mi> <msup> <mrow> <mo>(</mo> <mi>p</mi> <mo>/</mo> <msub> <mi>Zp</mi> <mn>0</mn> </msub> <mo>)</mo> </mrow> <mrow> <mi>N</mi> <mo>+</mo> <mn>1</mn> </mrow> </msup> </mrow> </mfrac> <mo>&rsqb;</mo> <mo>,</mo> </mrow>

其中,C为与净热吸附有关的常数,p为气藏压力,Z为气体压缩因子,p0为气体饱和吸附压力,N为吸附层数;

按如下公式得到表观渗透率:

<mrow> <msub> <mi>K</mi> <mi>a</mi> </msub> <mo>=</mo> <msub> <mi>K</mi> <mi>e</mi> </msub> <mo>&CenterDot;</mo> <mi>f</mi> <mrow> <mo>(</mo> <mi>K</mi> <mi>n</mi> <mo>)</mo> </mrow> <mo>+</mo> <msub> <mi>K</mi> <mi>s</mi> </msub> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>e</mi> </msub> <msubsup> <mi>R</mi> <mi>e</mi> <mn>2</mn> </msubsup> </mrow> <mrow> <mn>8</mn> <msub> <mi>&tau;</mi> <mi>m</mi> </msub> </mrow> </mfrac> <mo>&lsqb;</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <mi>&alpha;</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>-</mo> <mi>b</mi> <mo>)</mo> </mrow> <msup> <mi>Kn</mi> <mn>2</mn> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mi>K</mi> <mi>n</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>b</mi> <mo>&CenterDot;</mo> <mi>K</mi> <mi>n</mi> </mrow> </mfrac> <mo>&rsqb;</mo> <mo>+</mo> <msub> <mi>K</mi> <mi>s</mi> </msub> <mo>,</mo> </mrow>

其中,Ka为表观渗透率,Ks为表面扩散等效渗透率,f(Kn)为克努森函数,如下:

<mrow> <mi>f</mi> <mrow> <mo>(</mo> <mi>K</mi> <mi>n</mi> <mo>)</mo> </mrow> <mo>=</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <mi>&alpha;</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>-</mo> <mi>b</mi> <mo>)</mo> </mrow> <msup> <mi>Kn</mi> <mn>2</mn> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mi>K</mi> <mi>n</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>b</mi> <mo>&CenterDot;</mo> <mi>K</mi> <mi>n</mi> </mrow> </mfrac> <mo>,</mo> </mrow>

其中,b为滑脱系数,Kn为克努森数;

按如下公式得到校正的渗透率乘子:

<mrow> <msub> <mi>R</mi> <mi>K</mi> </msub> <mo>=</mo> <mfrac> <mrow> <mn>8</mn> <msub> <mi>&tau;</mi> <mi>m</mi> </msub> <msub> <mi>K</mi> <mi>a</mi> </msub> </mrow> <mrow> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <msubsup> <mi>R</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>=</mo> <mfrac> <mrow> <msub> <mi>&phi;</mi> <mi>e</mi> </msub> <msubsup> <mi>R</mi> <mi>e</mi> <mn>2</mn> </msubsup> <mo>&lsqb;</mo> <mn>1</mn> <mo>+</mo> <mfrac> <mrow> <mi>&alpha;</mi> <mrow> <mo>(</mo> <mn>4</mn> <mo>-</mo> <mi>b</mi> <mo>)</mo> </mrow> <msup> <mi>Kn</mi> <mn>2</mn> </msup> <mo>+</mo> <mrow> <mo>(</mo> <mn>4</mn> <mo>+</mo> <mi>&alpha;</mi> <mo>)</mo> </mrow> <mi>K</mi> <mi>n</mi> </mrow> <mrow> <mn>1</mn> <mo>-</mo> <mi>b</mi> <mo>&CenterDot;</mo> <mi>K</mi> <mi>n</mi> </mrow> </mfrac> <mo>&rsqb;</mo> <mo>+</mo> <mn>8</mn> <msub> <mi>&tau;</mi> <mi>m</mi> </msub> <msub> <mi>K</mi> <mi>s</mi> </msub> </mrow> <mrow> <msub> <mi>&phi;</mi> <mn>0</mn> </msub> <msubsup> <mi>R</mi> <mn>0</mn> <mn>2</mn> </msubsup> </mrow> </mfrac> <mo>,</mo> </mrow>

其中,RK为校正的渗透率乘子,φ0为初始条件固有孔隙度,R0为不考虑吸附层的孔隙半径。

9.一种页岩气藏数值模拟器输入参数的校正装置,其特征在于,包括:

吸附气体积模块,用于根据吸附气储气能力和吸附气密度得到岩石物理模型中吸附气体积;

吸附气孔隙度模块,用于根据岩石物理模型中吸附气体积得到岩石物理模型吸附气孔隙度;

输入参数模块,用于根据岩石物理模型有效孔隙度和岩石物理模型吸附气孔隙度,得到校正的页岩气藏数值模拟器输入参数。

10.根据权利要求9所述的页岩气藏数值模拟器输入参数的校正装置,其特征在于,输入参数模块具体用于:

根据岩石物理模型整体孔隙度和岩石物理模型吸附气孔隙度,得到校正的体相孔隙度;

根据校正的体相孔隙度得到校正的束缚水饱和度;

根据校正的束缚水饱和度得到校正的自由气饱和度;

根据岩石物理模型有效孔隙度和校正的体相孔隙度得到不同压力下校正的孔隙度乘子;

根据岩石物理模型有效孔隙度得到校正的渗透率乘子。

当前第2页1 2 3 
网友询问留言 已有0条留言
  • 还没有人留言评论。精彩留言会获得点赞!
1