Method for predicting amount of water-sealed gas in high-sulfur water-bearing gas reservoir

A method for predicting an amount of water-sealed gas in a high-sulfur water-bearing gas reservoir, comprising the following steps:

• S100, deriving a material balance equation of the high-sulfur water-bearing gas reservoir in consideration of water-sealed gas and water-soluble gas, which mainly comprises:
• S101, establishing a physical model of the high-sulfur water-bearing gas reservoir in consideration of water-sealed gas and water-soluble gas;
• S102, for a constant-volume gas reservoir in which a sum of a volume of initial gas and a volume of initial formation water is equal to a sum of a volume of remaining gas, an evolution volume of water soluble gas, a volume of irreducible water, a rock expansion volume, and a volume of remaining formation water, establishing the material balance equation;
• and
• S103, substituting each of the volume of remaining gas, the evolution volume of water soluble gas, the volume of irreducible water, the rock expansion volume, and the volume of remaining formation water in S102 into the material balance formula for representation, wherein the material balance equation of the high-sulfur water-bearing gas reservoir in consideration of water-sealed gas and water-soluble gas is given as follows:
• φ
• =
• 1
• -
• R
• -
• δ
• ⁢
• R
• F
• 1
• +
• β
• +
• (
• α
• -
• 1
• )
• ⁢
• R
• F
• ,
• wherein
• ⁢
• ⁢
• δ
• =
• K
• +
• V
• +
• ⁢
• K
• +
• +
• 1
• ,
• α
• =
• r
• s
• ⁢
• B
• g
• ,
• r
• s
• =
• 8.314
• Z
• s
• ⁢
• ρ
• w
• ⁢
• m
• H
• 2
• ⁢
• S
• ⁢
• T
• p
• ,
• β
• =
• B
• g
• ⁢
• S
• wi
• S
• gi
• ⁢
• (
• B
• w
• B
• wi
• ⁢
• r
• s
• -
• r
• si
• )
• -
• (
• C
• f
• +
• C
• w
• ⁢
• S
• wi
• )
• ⁢
• (
• p
• i
• -
• p
• )
• S
• gi
• ,
• and δ denotes a reservoir heterogeneity coefficient;
• φ denotes a dimensionless pseudopressure;
• R denotes a dimensionless recovery percent of a gas reservoir;
• F denotes a dimensionless water invasion constant;
• K+ denotes a dimensionless permeability ratio;
• V+ denotes a dimensionless volume ratio;
• p;
• denotes an initial formation pressure of the gas reservoir, in unit of MPa;
• p denotes a current formation pressure of the gas reservoir, in unit of MPa;
• Bg denotes a dimensionless volume coefficient of natural gas under a pressure p;
• Swi denotes a dimensionless initial formation water saturation;
• Sgi denotes a dimensionless initial gas saturation;
• Bg denotes a dimensionless volume coefficient of formation water under a pressure p;
• Bwi denotes a dimensionless volume coefficient of formation water under a pressure pi;
• rs denotes a solution gas-water ratio of hydrogen sulfide under a pressure p, in unit of m3/m3;
• rsi denotes a solution gas-water ratio of hydrogen sulfide under a pressure pi, in unit of m3/m3, Cf denotes a compressibility coefficient of rock, in unit of MPa−1;
• Cw denotes a compressibility coefficient of formation water, in unit of MPa−1;
• Zs denotes a dimensionless deviation factor of hydrogen sulfide;
• ρw denotes density of water, in unit of kg/m3;
• mH2S denotes solubility of hydrogen sulfide, in unit of mol/kg;
• and T denotes a temperature of the gas reservoir, in unit of K;
• S200, drawing a chart of water-sealed gas in the high-sulfur water-bearing gas reservoir by conducting iterative operation in the material balance equation of the high-sulfur water-bearing gas reservoir in consideration of water-sealed gas and water-soluble gas;
• and
• S300, calculating the amount of the water-sealed gas in the high-sulfur water-bearing gas reservoir, which specifically comprises:
• S301, taking R corresponding to φ=0 as a recovery factor EUR of the high-sulfur water-bearing gas reservoir in consideration of water-sealed gas and water-soluble gas;
• and
• S302, setting GB=(1−EUR)×G;
• wherein G denotes dynamic reserves, in unit of 108 m3, and GB denotes an amount of water-sealed gas, in unit of 108 m3.