CN108343426A - A kind of fracture-pore reservoir oil well output moisture content analogy method - Google Patents
A kind of fracture-pore reservoir oil well output moisture content analogy method Download PDFInfo
- Publication number
- CN108343426A CN108343426A CN201810191769.4A CN201810191769A CN108343426A CN 108343426 A CN108343426 A CN 108343426A CN 201810191769 A CN201810191769 A CN 201810191769A CN 108343426 A CN108343426 A CN 108343426A
- Authority
- CN
- China
- Prior art keywords
- node
- oil
- water
- moisture content
- fracture
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B47/00—Survey of boreholes or wells
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH DRILLING; MINING
- E21B—EARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
Abstract
The invention discloses a kind of fracture-pore reservoir oil well output moisture content analogy methods, include the following steps：S1, the description scheme and characterising parameter for initially setting up more fracture-cavity units oil reservoirs；S2, the fracture-cavity units description scheme and characterising parameter established in S1 are utilized, in conjunction with the flexible drive stage material balance relationship and moisture content model during oil reservoir development, use elastic expansion energy as the virtual injection rate of fracture-cavity units, simulation calculates each fracture-cavity units stage oil production and water yield, obtains the oil production, water yield and cumulative liquid production of oil well；S3, last oil production, water yield and cumulative liquid production according to obtained oil well, calculate the moisture content of producing well with the variation of cumulative liquid production.The inventive method is that the Production development moisture content analysis in fracture-pore reservoir flexible drive stage and understanding seam hole structure contribute new thinking.
Description
Technical field
The present invention relates to oil-gas field development field, especially a kind of fracture-pore reservoir oil well output moisture content analogy method.
Background technology
For fracture and vug carbonate reservoir, reservoir and flow space is dissolution porosity, corrosion hole and crack, fracture hole flow conductivity
By force, solution cavity is main reservoir space, and the grease gravitational differentiation in solution cavity is notable, and crack is main flow channel, and oil well is shown in
Moisture content rises very fast after water, and natural water drive reserves exploitation degree is poor, and recovery ratio is relatively low.It is found according to creation data, oil well life
More broken line forms of the tired production liquid curve of the tired production water-of production, as shown in Figure 1, the water-cut variation of oil well is divided into a, b, c, d tetra- in figure
A stage includes two kinds of general types of moisture content raising and lowering, the case where variation in three steps for water ratio in oil well,
Simulation combination can be carried out with three vug unit models：Assuming that oil well is gone into operation by 1# solution cavities, 2# solution cavities and 3# solution cavities difference
To 1# solution cavity feed flows, 1# solution cavities, 2# solution cavities outlet end be located on original oil-water level, solution cavity is to produce water after first oil-producing, molten
Fluid composition in hole is that water heavy wool is few, and the outlet end of 3# solution cavities is located under original oil-water level, and solution cavity is produced after first producing water
The fluid composition of oil, solution cavity is that the more water of oil are few；Oil well does not produce water in stage a, and 1# cavern waters do not overflow, and 2# cavern waters do not overflow
Go out, 3# solution cavity oil does not overflow, and well oil output is equal to grease overall expansion amount of the oily swell increment plus 2# solution cavities of 1# solution cavities；In rank
It is produced in section b for grease is same, 1# cavern waters overflow, and 2# cavern waters do not overflow, and 3# solution cavity oil does not overflow, and it is molten that stage water yield is equal to 1#
The water swell increment in hole adds the grease overall expansion amount of 3# solution cavities, and the oily swell increment that stage oil production is equal to 1# solution cavities adds 2# solution cavities
Grease overall expansion amount, due to 1# solution cavities production water cause water ratio in oil well rise a step；Grease is the same as production in stage c；1#
Cavern water overflows, and 2# cavern waters overflow, and 3# solution cavity oil does not overflow, and the water swell increment that stage water yield is equal to 1# solution cavities is molten plus 2#
The water swell increment in hole adds the grease overall expansion amount of 3# solution cavities, and the oily swell increment that stage oil production is equal to 1# solution cavities adds 2# solution cavities
Oily swell increment, due to 2# solution cavities production water cause water ratio in oil well rise a step；Grease is the same as production, 1# solution cavities in stage d
Water overflows, and 2# cavern waters overflow, and 3# solution cavity oil spill goes out, and stage water yield is equal to water of the water swell increment plus 2# solution cavities of 1# solution cavities
Swell increment adds the water swell increment of 3# solution cavities, and stage oil production is equal to oily swell increment of the oily swell increment plus 2# solution cavities of 1# solution cavities
In addition the oily swell increment of 3# solution cavities, since 3# solution cavity oil-producings cause water ratio in oil well to decline a step.By to three solution cavities
Cellular manufacture is dynamically analyzed, and the version that can release seam hole structure and moisture content in oil reservoir has certain relationship, so
The present invention proposes a kind of new method to simulate the Production development moisture content in oil reservoir flexible drive stage, to analysis and understanding seam
Hole structure.
Invention content
It is an object of the invention to overcome the prior art, a kind of fracture-pore reservoir oil well generation moisture content mould is provided
Quasi- method contributes new think of for the analysis of Production development moisture content and understanding seam hole structure in fracture-pore reservoir flexible drive stage
Road.
The purpose of the present invention is achieved through the following technical solutions：A kind of fracture-pore reservoir oil well output moisture content simulation side
Method includes the following steps：S1, the description scheme and characterising parameter for initially setting up more fracture-cavity units oil reservoirs；S2, using in step S1
The fracture-cavity units description scheme and characterising parameter of foundation, in conjunction with the flexible drive stage material balance relationship during oil reservoir development
With moisture content model, use elastic expansion amount as the virtual injection rate of fracture-cavity units, simulation calculates the rank of each fracture-cavity units
Section oil production and water yield, obtain the oil production, water yield and cumulative liquid production of oil well；The oil well that S3, last basis obtain
Oil production, water yield and cumulative liquid production calculate the moisture content of producing well with the variation of cumulative liquid production.
Preferably, the description scheme of fracture-cavity units mainly uses tree to indicate in the step S1, the fracture hole list
Solution cavity in member indicates that the crack in fracture-cavity units is indicated with connection path with node, and the tree is closed according to child node
It is built-up to join backward direction root node mode.
Preferably, the characterising parameter of fracture-cavity units includes fracture-cavity units public attribute parameter, solution cavity section in the step S1
Dot element static attribute parameter and solution cavity node dynamic attribute parameter；The fracture-cavity units public attribute parameter includes pressure ginseng
Number, volume factor, the compressed coefficient and fluid viscosity：
Pressure parameter includes strata pressure Pi, current formation pressure P and crude oil saturation pressure Pb, MPa,
Volume factor includes crude oil initial volume coefficient BoiWith water flooding initial volume coefficient Bwi, zero dimension,
The compressed coefficient includes oil compressibility Co, water flooding compressed coefficient CwWith rock compressibility Cf, MPa-1,
Fluid viscosity includes viscosity of crude μoWith water flooding viscosity, muw, mPa.s；
The solution cavity node unit static attribute parameter, including the original fluid reserves of solution cavity and outlet port, wherein j generations
J-th of node of table；
The original fluid reserves of solution cavity include the reserves N of crude oilj, the reserves W of water floodingjAnd the total volume V of solution cavityj,
Middle Vj=Nj+Wj, m3；
Outlet port includes lower spilling coefficient Cd,jWith upper spilling coefficient Cu,j, the lower spilling coefficient Cd,jFor oil-water interfaces
Rise the ratio between water volume and the total volume occurred when water flooding is overflowed, the upper spilling coefficient Cu,jDecline for oil-water interfaces and occurs
The ratio between water volume and total volume when crude oil overflows；
The solution cavity node dynamic attribute parameter includes node oil-water interfaces coefficient Cj, node time step flow into oil volume
Vio,j, node time step flow into water volume Viw,j, node time step effluent oil water volume Vow,jMoisture content is flowed out with the time step of node
fw,j, wherein node oil-water interfaces system CjNumber is the ratio between water volume and total volume under oil-water interfaces.
Preferably, the step S2 is using the more fracture-cavity units tree models established in step S1, using Numerical-Mode
Quasi- method, using a pressure drop step delta p as time step, according to material balance relationship, solution cavity influx is equal to discharge, from tree
The grease that the bottom leaf node of type structure starts calculate node virtually flows into liquid measure and outflow moisture content, and successively scanning calculates, directly
To root node is traced back to, the outflow liquid measure and outflow moisture content, multiple scanning node that obtain oil well reach target until overall presure drop
Value terminates to calculate；The specific method is as follows for it：S21, using a pressure drop step-length as time step, scanned since tree bottom,
The grease of calculate node virtually flows into liquid measure and flows into moisture content, and it includes front end child node that interior joint grease, which virtually flows into liquid measure,
The volume expansion amount of the external influx and node self residual fluid of generation；S22, decision node are flowing into front and back section in the process
The relationship of the oil-water interface position and upper and lower spill point of point calculates sinking in conjunction with the seam hole structure of node in varied situations
Degree and outflow moisture content, obtain the effluent oil water volume of node；The boundary coefficient of S23, more new node, the stage of rear end father node
Influx；S24, it scans up, after executing step S21, step S22 to upper layer father node, step S23 is executed, until tracing back to
Root node obtains the outflow liquid measure and outflow moisture content of oil well；S25, when pressure drop miss the mark value, repeat step S21
To S24, when pressure drop reaches setting value, execution step S26；S26, terminate to calculate.
Preferably, the virtual influx of oil of the step S21 interior joints j：vo,j=vio,j+(1-Cj)Vj(Co+Cf) Δ p,
The virtual influx of water of node j：vw,j=viw,j+CjVj(Cw+Cf) Δ p,
The inflow liquid measure of node j：vow,j=vo,j+vw,j,
The inflow moisture content of node j：fiw,j=vw,j/vow,j, wherein j represents j-th of node.
Preferably, the outflow moisture content of calculate node needs decision node before flowing into liquid measure and flowing into the step S22
The relationship of oil-water interfaces and upper and lower spill point, wherein j represent j-th of node during afterwards：
If (a) crude oil water termination is less than lower spill point and the interface being added after flowing into water is less than lower spill point, i.e. Cj
＜ Cd,jAndThen node oil-producing does not produce water, outflow moisture content fw,j=0, execute step S23；
If (b) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil higher than upper spill point, i.e. Cj
＞ Cu,jAndThen node only produces water not oil-producing, outflow moisture content fw,j=1, execute step S23；
If (c) crude water interface is between upper and lower spill point, i.e. Cd,j≤Cj≤Cu,j, then node grease is the same as production, use
Submergence model calculates outflow moisture content, executes step S23；
If (d) crude oil water termination less than lower spill point and adds the interface after flowing into water to be higher than lower spill point, i.e. Cj＜
Cd,jAndLower spill point is crossed in then node oil-water interfaces rising, and water is not produced in pressure time step leading portion oil-producing,
In pressure time step back segment grease with producing, influx can be decomposed into two parts, and first part is that oil-water interfaces rise to lower spilling
Point, second part enter between upper and lower spill point, and the outflow moisture content after first part's influx inflow of node j is according to (a)
Situation calculates, and the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains egress j
Outflow moisture content after, execute step S23；
Calculate the water influx that first part's oil-water interfaces rise to lower spill point：vw1=Vj(Cj-Cd,j), by flow
Enter the energy expenditure ratio that amount calculates first part：K=vw1/vw,j, so the node input parameter of first part is：
The crude oil of node j virtually flows into volume：vo,j'=kvo,j,
The water flooding of node j virtually flows into volume：vw,j'=kvw,j,
The inflow liquid measure of node j：vow,j'=kvow,j
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
The node input parameter of second part is：
The crude oil of node j virtually flows into volume：vo,j'=(1-k) vo,j,
The water flooding of node j virtually flows into volume：vw,j'=(1-k) vw,j,
The inflow liquid measure of node j：vow,j'=(1-k) vow,j,
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
If (e) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil less than upper spill point, i.e. Cj
＞ Cu,jAndThen node oil-water interfaces are decreased past spill point, only produce water in pressure time step leading portion and do not produce
Oil, in pressure time step back segment grease with producing, influx can be decomposed into two parts, and first part is that oil-water interfaces drop to overflow
Go out a little, second part enters between upper and lower spill point, the first part influx of node j flow into after outflow moisture content according to
(b) situation calculates, and the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains section
After the outflow moisture content of point j, step S23 is executed；
Calculate the oily influx that first part's oil-water interfaces rise to lower spill point：vo1=Vj(Cj-Cu,j), by oil stream
Enter the energy expenditure ratio that amount calculates first part：K=vo1/vo,j, so the node input parameter of first part is：
The crude oil of node j virtually flows into volume：vo,j'=kvo,j,
The water flooding of node j virtually flows into volume：vw,j'=kvw,j,
The inflow liquid measure of node j：vow,j'=kvow,j
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
The node input parameter of second part is：
The crude oil of node j virtually flows into volume：vo,j'=(1-k) vo,j,
The water flooding of node j virtually flows into volume：vw,j'=(1-k) vw,j,
The inflow liquid measure of node j：vow,j'=(1-k) vow,j,
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
Preferably, the specific method of outflow moisture content is calculated in the step S22 with submergence model：
(1) initial value of the pressure drop time step of calculate node j submergences：Sj=(Cj-Cd,j)/(Cu,j-Cd,j)；
(2) submergence model is solved, the submergence of node is obtainedNew value：
(3) with sinking angle valueThe pressure drop time step of new value calculate node flows out moisture content：Wherein
fiw, j is the inflow moisture content of node j, a=μw/μo, b=VA,j/vow,j, VA,jFor node j grease with producing region volume, VA,j=
Vj(Cu,j-Cd,j)。
Preferably, the stage influx of the boundary coefficient and rear end father node of the step S23 more new nodes, wherein j generations
J-th of node of table, k represent the father node of the rear ends node j, and the specific method is as follows：
The original oil volumes of node j：Vo,j=Vj(1-Cj),
The original water volumes of node j：Vw,j=VjCj,
The remaining oil volumes of node j：
The remaining water volumes of node j：
More new node j boundary coefficients：
Update rear end father node k oil influxs：vio,k *=vio,k+vow,j(1-fw,j)。
The present invention has the following advantages：
The seam hole structure that understanding oil reservoir is analyzed by Production development is one of the main target of fracture-pore reservoir dynamic monitoring,
The space structure complex of fracture-pore reservoir is difficult to portray and be described using conventional reservoir modeling technology.The present invention uses
Volume equivalent way describes the size exclusion solution cavity morphology influence of solution cavity Reservoir Body, is described using spill point and oil-water interfaces mode
The variation that solution cavity output is constituted, the connection model between more solution cavities and producing well is described using tree, is closed based on material balance
System quickly calculates the variation of the cumulative liquid production and moisture content under specific pressure drop.The present invention is to identify the moisture content of fracture-pore reservoir
The fracture-cavity units structure and evaluation control that incidence relation, detection oil well between changing pattern and fracture-cavity units tactic pattern are connected to
Reserves provide a kind of analysis approach of quantification.
Description of the drawings
Fig. 1 is the production hydration of more solution cavities into oil well multi-step moisture content schematic diagram；
Fig. 2 is double fracture-cavity units moisture content analog results；
Fig. 3 is the tired production Fluid Dynamics result of double fracture-cavity units；
Fig. 4 is three fracture-cavity units moisture content analog results；
Fig. 5 is the tired production Fluid Dynamics result of three fracture-cavity units.
Specific implementation mode
The present invention will be further described below in conjunction with the accompanying drawings, but protection scope of the present invention is not limited to following institute
It states.
A kind of fracture-pore reservoir oil well output moisture content analogy method, includes the following steps：S1, more fracture hole lists are initially set up
The description scheme and characterising parameter of first oil reservoir；S2, using the fracture-cavity units description scheme and characterising parameter established in S1, in conjunction with oil
The flexible drive stage material balance relationship and moisture content model in development process are hidden, uses elastic energy as the void of fracture-cavity units
Quasi- injection rate, simulation calculate each fracture-cavity units stage oil production and water yield, obtain the oil production of oil well, water yield and tired
Product Liquid output；S3, last oil production, water yield and cumulative liquid production according to obtained oil well, calculate the moisture content of producing well
With the variation of cumulative liquid production.
Preferably, the description scheme of fracture-cavity units mainly uses tree to indicate in the step S1, the fracture hole list
Solution cavity in member indicates that the crack in fracture-cavity units is indicated with connection path with node, and the tree is closed according to child node
It is built-up to join backward direction root node mode.
Preferably, the characterising parameter of fracture-cavity units includes fracture-cavity units public attribute parameter, solution cavity section in the step S1
Dot element static attribute parameter and solution cavity node dynamic attribute parameter；The fracture-cavity units public attribute parameter includes pressure ginseng
Number, volume factor, the compressed coefficient and fluid viscosity：
Pressure parameter includes strata pressure Pi, current formation pressure P and crude oil saturation pressure Pb, MPa,
Volume factor includes crude oil initial volume coefficient BoiWith water flooding initial volume coefficient Bwi, zero dimension,
The compressed coefficient includes oil compressibility Co, water flooding compressed coefficient CwWith rock compressibility Cf, MPa-1,
Fluid viscosity includes viscosity of crude μoWith water flooding viscosity, muw, mPa.s；
The solution cavity node unit static attribute parameter, including the original fluid reserves of solution cavity and outlet port, wherein j generations
J-th of node of table,
The original fluid reserves of solution cavity include the reserves N of crude oilj, the reserves W of water floodingjAnd the total volume V of solution cavityj,
Middle Vj=Nj+Wj, m3,
Outlet port includes lower spilling coefficient Cd,jWith upper spilling coefficient Cu,j, the lower spilling coefficient Cd,jFor oil-water interfaces
Rise the ratio between water volume and the total volume occurred when water flooding is overflowed, the upper spilling coefficient Cu,jDecline for oil-water interfaces and occurs
The ratio between water volume and total volume when crude oil overflows；
The solution cavity node dynamic attribute parameter includes node oil-water interfaces coefficient Cj, node time step flow into oil volume
Vio,j, node time step flow into water volume Viw,j, node time step effluent oil water volume Vow,jMoisture content is flowed out with the time step of node
fw,j, wherein node oil-water interfaces system CjNumber is the ratio between water volume and total volume under oil-water interfaces.
Preferably, the step S2 is using the more fracture-cavity units tree models established in step S1, using Numerical-Mode
Quasi- method, using a pressure drop step-length as time step, according to material balance relationship, solution cavity influx is equal to discharge, from tree-shaped
The grease that the bottom leaf node of structure starts calculate node virtually flows into liquid measure and outflow moisture content, and successively scanning calculates, until
Root node is traced back to, the outflow liquid measure and outflow moisture content, multiple scanning node for obtaining oil well reach desired value until overall presure drop,
Terminate to calculate；The specific method is as follows for it：S21, using a pressure drop step delta p as time step, scanned since tree bottom,
The virtual influx of grease and inflow moisture content of calculate node, the virtual influx of interior joint grease include that front end child node generates
External influx and node self residual fluid volume expansion amount；S22, the decision node process before and after flowing into liquid measure and flowing into
The relationship of the oil-water interfaces of interior joint and upper and lower spill point calculates sinking in conjunction with the seam hole structure of node in varied situations
Degree and outflow moisture content, obtain the effluent oil water volume of node；The boundary coefficient of S23, more new node, the stage of rear end father node
Influx；S24, it scans up, after executing S21, S22 step to upper layer father node, step S23 is executed, until tracing back to root section
Point obtains outflow liquid measure and accumulates numerical value；S25, when pressure drop miss the mark value, repeat step S21, set when pressure drop reaches
Definite value executes S26 steps；S26, terminate to calculate.
Preferably, the virtual influx of oil of the step S21 interior joints j：vo,j=vio,j+(1-Cj)Vj(Co+Cf) Δ p,
The virtual influx of water of node j：vw,j=viw,j+CjVj(Cw+Cf) Δ p,
The inflow liquid measure of node j：vow,j=vo,j+vw,j,
The inflow moisture content of node j：fiw,j=vw,j/vow,j, wherein j represents j-th of node.
Preferably, the outflow moisture content of calculate node needs decision node before flowing into liquid measure and flowing into the step S22
The relationship of oil-water interfaces and upper and lower spill point, wherein j represent j-th of node during afterwards：
If (a) crude oil water termination is less than lower spill point and the interface being added after flowing into water is less than lower spill point, i.e. Cj
＜ Cd,jAndThen node oil-producing does not produce water, outflow moisture content fw,j=0, execute step S23；
If (b) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil higher than upper spill point, i.e. Cj
＞ Cu,jAndThen node only produces water not oil-producing, outflow moisture content fw,j=1, execute step S23；
If (c) crude water interface is between upper and lower spill point, i.e. Cd,j≤Cj≤Cu,j, then node grease is the same as production, use
Submergence model calculates outflow moisture content, executes step S23；
If (d) crude oil water termination less than lower spill point and adds the interface after flowing into water to be higher than lower spill point, i.e. Cj＜
Cd,jAndLower spill point is crossed in then node oil-water interfaces rising, and water is not produced in pressure time step leading portion oil-producing,
In pressure time step back segment grease with producing, influx can be decomposed into two parts, and first part is that oil-water interfaces rise to lower spilling
Point, second part enter between upper and lower spill point, and the outflow moisture content after first part's influx inflow of node j is according to (a)
Situation calculates, and the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains egress j
Outflow moisture content after, execute step S23；
Calculate the water influx that first part's oil-water interfaces rise to lower spill point：vw1=Vj(Cj-Cd,j), by flow
Enter the energy expenditure ratio that amount calculates first part：K=vw1/vw,j, so the node input parameter of first part is：
The crude oil of node j virtually flows into volume：vo,j'=kvo,j,
The water flooding of node j virtually flows into volume：vw,j'=kvw,j,
The inflow liquid measure of node j：vow,j'=kvow,j
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
The node input parameter of second part is：
The crude oil of node j virtually flows into volume：vo,j'=(1-k) vo,j,
The water flooding of node j virtually flows into volume：vw,j'=(1-k) vw,j,
The inflow liquid measure of node j：vow,j'=(1-k) vow,j,
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
If (e) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil less than upper spill point, i.e. Cj
＞ Cu,jAndThen node oil-water interfaces are decreased past spill point, only produce water in pressure time step leading portion and do not produce
Oil, in pressure time step back segment grease with producing, influx can be decomposed into two parts, and first part is that oil-water interfaces drop to overflow
Go out a little, second part enters between upper and lower spill point, the first part influx of node j flow into after outflow moisture content according to
(b) situation calculates, and the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains section
After the outflow moisture content of point j, step S23 is executed；
Calculate the oily influx that first part's oil-water interfaces rise to lower spill point：vo1=Vj(Cj-Cu,j), by oil stream
Enter the energy expenditure ratio that amount calculates first part：K=vo1/vo,j, so the node input parameter of first part is：
The crude oil of node j virtually flows into volume：vo,j'=kvo,j,
The water flooding of node j virtually flows into volume：vw,j'=kvw,j,
The inflow liquid measure of node j：vow,j'=kvow,j
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
The node input parameter of second part is：
The crude oil of node j virtually flows into volume：vo,j'=(1-k) vo,j,
The water flooding of node j virtually flows into volume：vw,j'=(1-k) vw,j,
The inflow liquid measure of node j：vow,j'=(1-k) vow,j,
The inflow moisture content of node j：fiw,j'=vw,j'/vow,j；
Preferably, the specific method of outflow moisture content is calculated in the step S22 with submergence model：
(1) initial value of the pressure drop time step of calculate node j submergences：Sj=(Cj-Cd,j)/(Cu,j-Cd,j)；
(2) submergence model is solved, the submergence of node is obtainedNew value：
(3) with sinking angle valueThe pressure drop time step of new value calculate node flows out moisture content：Wherein
fiw,jFor the inflow moisture content of node j, a=μw/μo, b=VA,j/vow,j, VA,jFor node j grease with producing region volume, by overflowing
Coefficient calculates VA,j=Vj(Cu,j-Cd,j)。
Preferably, the stage influx of the boundary coefficient and rear end father node of the step S23 more new nodes, wherein j generations
J-th of node of table, k represent the father node of the rear ends node j, and the specific method is as follows：
The original oil volumes of node j：Vo,j=Vj(1-Cj),
The original water volumes of node j：Vw,j=VjCj,
The remaining oil volumes of node j：
The remaining water volumes of node j：
More new node j boundary coefficients：
Update rear end father node k oil influxs：vio,k *=vio,k+vow,j(1-fw,j)。
Instance analysis：
Underlying parameter is set：The volume factor 1.05 of crude oil, the volume factor 1.02 of water flooding, oil compressibility
0.0015MPa-1, water flooding compressed coefficient 0.0005MPa-1, rock compressibility 0.0004MPa-1, viscosity of crude 25mPa.s,
Water flooding viscosity 0.25mPa.s, original formation pressure 55MPa, current formation pressure 40MPa.
As shown in Figure 2 and Figure 3, a kind of structure of double fracture-cavity units is shown in Table 1 with the definition of grease reserves, and oil well is connected to fracture hole list
The moisture content of member 1, unit 1 is oil well output moisture content, and fracture-cavity units 2 calculate each fracture-cavity units to 1 feed flow of unit, simulation
Water-cut variation, accumulation production water variation, reflect the stepped ramp type change procedure of fracture-pore reservoir water ratio in oil well.Fracture-cavity units 1
Original oil-water level be less than lower spill point, be water breakthrough after first oil-producing after operation；Original oil-water level in fracture-cavity units 2 is less than
Upper spill point, be higher than lower spill point, go into operation i.e. start produce water, with oil-water interfaces decline moisture content decline, fracture-cavity units 2 it is aqueous
Rate decline causes the moisture content of unit 1 to decline.
1 two fracture-cavity units structures of table are defined with grease reserves
As shown in Figure 4, Figure 5, a kind of structure of three fracture-cavity units is shown in Table 2 with the definition of grease reserves, and oil well is connected to fracture hole
Unit 1, fracture-cavity units 2 and unit 3 are to 1 feed flow of unit, and simulation calculates the water-cut variation of each fracture-cavity units, accumulation production water becomes
Change, reflects the stepped ramp type change procedure of fracture-pore reservoir water ratio in oil well.The original oil-water level of fracture-cavity units 2 and unit 3
It is water breakthrough after first oil-producing less than lower spill point, after operation, because grease reserves and spill point relative position are different, unit 2 and unit 3
Water breakthrough opportunity it is inconsistent；The original oil-water level of fracture-cavity units 1 is higher than lower spill point, goes into operation and starts to produce water, but with unit
2, after the oil injection of unit 3, the moisture content of unit 1 falls to reduced levels, and the later stage leads to unit because of the production water of unit 2, unit 3
The rising of 1 moisture content.
The structure of 2 three fracture-cavity units of table is defined with grease reserves
What has been described above is only a preferred embodiment of the present invention, it is noted that for those of ordinary skill in the art
For, under the premise of different disengaging the invention are conceived, various modifications and improvements can be made, these belong to this hair
Bright protection domain.
Claims (8)
1. a kind of fracture-pore reservoir oil well output moisture content analogy method, it is characterised in that：Include the following steps：
S1, the description scheme and characterising parameter for initially setting up more fracture-cavity units oil reservoirs；
S2, using the fracture-cavity units description scheme and characterising parameter established in step S1, in conjunction with the elasticity during oil reservoir development
Driving stage material balance relationship and moisture content model use elastic expansion amount as the virtual injection rate of fracture-cavity units, simulation meter
Each fracture-cavity units stage oil production and water yield are calculated, the oil production, water yield and cumulative liquid production of oil well are obtained；
S3, last oil production, water yield and cumulative liquid production according to obtained oil well, calculate the moisture content of producing well with accumulation
The variation of Liquid output.
2. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 1, it is characterised in that：It is described
The description scheme of fracture-cavity units mainly uses tree to indicate in step S1, the solution cavity node table in the fracture-cavity units
Show, the crack in fracture-cavity units is indicated with connection path, and the tree is directed toward root node side according to child node association backward
Formula is built-up.
3. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 2, it is characterised in that：It is described
In step S1 the characterising parameter of fracture-cavity units include fracture-cavity units public attribute parameter, solution cavity node unit static attribute parameter and
Solution cavity node dynamic attribute parameter；
The fracture-cavity units public attribute parameter includes pressure parameter, volume factor, the compressed coefficient and fluid viscosity,
Pressure parameter includes strata pressure Pi, current formation pressure P and crude oil saturation pressure Pb, MPa,
Volume factor includes crude oil initial volume coefficient BoiWith water flooding initial volume coefficient Bwi, zero dimension,
The compressed coefficient includes oil compressibility Co, water flooding compressed coefficient CwWith rock compressibility Cf, MPa-1,
Fluid viscosity includes viscosity of crude μoWith water flooding viscosity, muw, mPa.s；
The solution cavity node unit static attribute parameter, including the original fluid reserves of solution cavity and outlet port, wherein j represent
J node,
The original fluid reserves of solution cavity include the reserves N of crude oilj, the reserves W of water floodingjAnd the total volume V of solution cavityj, wherein Vj=
Nj+Wj, m3；
Outlet port includes lower spilling coefficient Cd,jWith upper spilling coefficient Cu,j, the lower spilling coefficient Cd,jRise for oil-water interfaces
The ratio between water volume and the total volume when water flooding is overflowed, the upper spilling coefficient C occursu,jDecline for oil-water interfaces and crude oil occurs
The ratio between water volume and total volume when spilling；
The solution cavity node dynamic attribute parameter includes node oil-water interfaces coefficient Cj, node time step flow into oil volume Vio,j, section
The time step of point flows into water volume Viw,j, node time step effluent oil water volume Vow,jWith the time step outflow moisture content f of nodew,j,
In, node oil-water interfaces system CjNumber is the ratio between water volume and total volume under oil-water interfaces.
4. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 3, it is characterised in that：It is described
Step S2 utilizes the more fracture-cavity units tree models established in step S1, using the method for numerical simulation, a pressure drop
Step delta p is as time step, and according to material balance relationship, solution cavity influx is equal to discharge, from the bottom leaf node of tree
The grease for starting calculate node virtually flows into liquid measure and outflow moisture content, and successively scanning calculates, until tracing back to root node, obtains
The outflow liquid measure and outflow moisture content, multiple scanning node of oil well reach desired value until overall presure drop, terminate to calculate；Its specific side
Method is as follows：
S21, using a pressure drop step delta p as time step, scanned since tree bottom, the grease virtual stream of calculate node
Enter liquid measure and flow into moisture content, it includes the external influx and section that front end child node generates that interior joint grease, which virtually flows into liquid measure,
The volume expansion amount of point self residual fluid；
The relationship of S22, the decision node oil-water interface position of process interior joint and upper and lower spill point before and after inflow, in difference
In the case of, submergence and outflow moisture content are calculated in conjunction with the seam hole structure of node, obtains the effluent oil water volume of node；
The boundary coefficient of S23, more new node, the stage influx of rear end father node；
S24, it scans up, after executing step S21, step S22 to upper layer father node, step S23 is executed, until tracing back to root section
Point obtains the outflow liquid measure and outflow moisture content of oil well；
S25, when pressure drop miss the mark value, repeat step S21 to S24, when pressure drop reaches setting value, execute step S26；
S26, terminate to calculate.
5. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 4, it is characterised in that：It is described
The virtual influx of oil of step S21 interior joints j：vo,j=vio,j+(1-Cj)Vj(Co+Cf) Δ p,
The virtual influx of water of node j：vw,j=viw,j+CjVj(Cw+Cf) Δ p,
The inflow liquid measure of node j：vow,j=vo,j+vw,j,
The inflow moisture content of node j：fiw,j=vw,j/vow,j, wherein j represents j-th of node.
6. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 5, it is characterised in that：It is described
In step S22 the outflow moisture content of calculate node need decision node before and after inflow the oil-water interface position of process interior joint with
The relationship of upper and lower spill point, wherein j represent j-th of node：
If (a) crude oil water termination is less than lower spill point and the interface being added after flowing into water is less than lower spill point, i.e. Cj＜ Cd,j
AndThen node oil-producing does not produce water, outflow moisture content fw,j=0, execute step S23；
If (b) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil higher than upper spill point, i.e. Cj＞ Cu,j
AndThen node only produces water not oil-producing, outflow moisture content fw,j=1, execute step S23；
If (c) crude water interface is between upper and lower spill point, i.e. Cd,j≤Cj≤Cu,j, then node grease is with producing, with sinking
It spends model and calculates outflow moisture content, execute step S23；
If (d) crude oil water termination less than lower spill point and adds the interface after flowing into water to be higher than lower spill point, i.e. Cj＜ Cd,jAndLower spill point is crossed in then node oil-water interfaces rising, water is not produced in pressure time step leading portion oil-producing, in pressure
For time step back segment grease with producing, influx can be decomposed into two parts, and first part is that oil-water interfaces rise to lower spill point, and second
Part enters between upper and lower spill point, and the outflow moisture content after first part's influx inflow of node j is according to (a) situation meter
It calculates, the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains the outflow of egress j
After moisture content, step S23 is executed；
If (e) crude oil water termination higher than upper spill point and subtracts the interface after flowing into oil less than upper spill point, i.e. Cj＞ Cu,j
AndThen node oil-water interfaces are decreased past spill point, only produce water not oil-producing in pressure time step leading portion, are pressing
Power time step back segment grease is with producing, and influx can be decomposed into two parts, and first part is that oil-water interfaces drop to upper spill point, the
Two parts enter between upper and lower spill point, and the outflow moisture content after first part's influx inflow of node j is according to (b) situation
It calculates, the outflow moisture content after then the second part influx of node j flows into is calculated according to (c) situation, obtains the stream of egress j
After going out moisture content, step S23 is executed.
7. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 6, it is characterised in that：It is described
The specific method of outflow moisture content is calculated in step S22 using submergence model：
(1) initial value of the pressure drop time step of calculate node j submergences：Sj=(Cj-Cd,j)/(Cu,j-Cd,j)；
(2) submergence model is solved, the submergence of node is obtainedNew value：
(3) with sinking angle valueThe pressure drop time step of new value calculate node flows out moisture content：Wherein fiw,j
For the inflow moisture content of node j, a=μw/μo, b=VA,j/vow,j, VA,jFor node j grease with producing region volume, VA,j=Vj
(Cu,j-Cd,j)。
8. a kind of fracture-pore reservoir oil well output moisture content analogy method according to claim 7, it is characterised in that：It is described
The boundary coefficient of step S23 more new nodes and the stage influx of rear end father node, wherein j represent j-th of node, and k represents section
The father node of the rear ends point j, the specific method is as follows：
The original oil volumes of node j：Vo,j=Vj(1-Cj),
The original water volumes of node j：Vw,j=VjCj,
The remaining oil volumes of node j：
The remaining water volumes of node j：
More new node j boundary coefficients：
Update rear end father node k oil influxs：vio,k *=vio,k+vow,j(1-fw,j)。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810191769.4A CN108343426B (en) | 2018-03-08 | 2018-03-08 | Simulation method for water content output of fracture-cavity type oil reservoir oil well |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810191769.4A CN108343426B (en) | 2018-03-08 | 2018-03-08 | Simulation method for water content output of fracture-cavity type oil reservoir oil well |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108343426A true CN108343426A (en) | 2018-07-31 |
| CN108343426B CN108343426B (en) | 2021-05-18 |
Family
ID=62956702
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810191769.4A Active CN108343426B (en) | 2018-03-08 | 2018-03-08 | Simulation method for water content output of fracture-cavity type oil reservoir oil well |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108343426B (en) |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102518414A (en) * | 2011-12-28 | 2012-06-27 | 西南石油大学 | Test method for fracture-cavity carbonate condensate gas reservoir water injection substituting gas experiment |
| CN104790926A (en) * | 2015-03-20 | 2015-07-22 | 中国石油大学(北京) | Fracture-vug type oil reservoir water-flooding development effect evaluation method |
| CN105317407A (en) * | 2015-10-15 | 2016-02-10 | 中国石油天然气股份有限公司 | Development method of untabulated reservoir in extra-high water-cut period |
| CN105512458A (en) * | 2015-11-24 | 2016-04-20 | 中国石油天然气股份有限公司 | Comprehensive analyzing method and apparatus for fracture-vug type carbonate reservoir exhaustion exploitation features |
| CN105587298A (en) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | Water content retrieval method for water-driven oil well with multi-stream-tube mode |
| RU2594402C1 (en) * | 2015-08-28 | 2016-08-20 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Method for subsequent flooding of stratified reservoir |
| CN106837297A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | A kind of method for recognizing inter well connectivity and profit dynamic prediction |
-
2018
- 2018-03-08 CN CN201810191769.4A patent/CN108343426B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102518414A (en) * | 2011-12-28 | 2012-06-27 | 西南石油大学 | Test method for fracture-cavity carbonate condensate gas reservoir water injection substituting gas experiment |
| CN104790926A (en) * | 2015-03-20 | 2015-07-22 | 中国石油大学(北京) | Fracture-vug type oil reservoir water-flooding development effect evaluation method |
| RU2594402C1 (en) * | 2015-08-28 | 2016-08-20 | Открытое акционерное общество "Татнефть" им. В.Д. Шашина | Method for subsequent flooding of stratified reservoir |
| CN105317407A (en) * | 2015-10-15 | 2016-02-10 | 中国石油天然气股份有限公司 | Development method of untabulated reservoir in extra-high water-cut period |
| CN105512458A (en) * | 2015-11-24 | 2016-04-20 | 中国石油天然气股份有限公司 | Comprehensive analyzing method and apparatus for fracture-vug type carbonate reservoir exhaustion exploitation features |
| CN105587298A (en) * | 2015-12-17 | 2016-05-18 | 西南石油大学 | Water content retrieval method for water-driven oil well with multi-stream-tube mode |
| CN106837297A (en) * | 2016-12-22 | 2017-06-13 | 中国石油天然气股份有限公司 | A kind of method for recognizing inter well connectivity and profit dynamic prediction |
Non-Patent Citations (1)
| Title |
|---|
| 赵艳艳等: "单缝单洞缝洞型碳酸盐岩底水油藏弹性开采机理研究", 《西安石油大学学报（自然科学版）》 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108343426B (en) | 2021-05-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108229051B (en) | Method for predicting recovery ratio of air foam flooding of oil reservoir | |
| CN106545336A (en) | Consider the Productivity of tight gas reservoir seepage flow mechanism | |
| CN105626006A (en) | CO2 drive technological limit well spacing determination method for low-permeability oil reservoir | |
| CN104975827B (en) | Predict the material balance method of carbon dioxide flooding oil reservoir index | |
| CN104504230A (en) | Estimation method for recovery ratio and limit drainage radius of low-permeability gas well | |
| CN106401535B (en) | A kind of method of determining coal bed gas well mining intensity | |
| CN108518218B (en) | Unconventional oil and gas reservoir multi-section fracturing horizontal well single well dynamic reserve determination method | |
| CN111709847A (en) | Method for predicting recoverable reserves of top gas side water reservoir and evaluating development effect | |
| CN106555565A (en) | A kind of low pressure and easy leaking loses the cementing method that well positive injection is instead squeezed | |
| CN104675371A (en) | Complex oil displacement experimental method by means of alternately injecting gel and polymer solution after polymer flooding and follow-up water flooding | |
| CN110242263A (en) | Two or three combine the recovery factor calculation method under development mode | |
| CN109858177A (en) | A kind of horizontal well with bottom water reservoir water drive method for numerical simulation based on quasi- streamline method | |
| CN106640021A (en) | Calculating method and device for post-fracture blow-off parameters | |
| CN114233270B (en) | Bottom water heavy oil reservoir horizontal well productivity prediction method | |
| CN107437127A (en) | A kind of oil well stop-spraying Formation pressure prediction method | |
| CN108343426A (en) | A kind of fracture-pore reservoir oil well output moisture content analogy method | |
| CN111626001B (en) | Method for improving refined water injection of oil extraction well | |
| CN105840156B (en) | A kind of more slug compound displacement techniques | |
| CN110397436B (en) | Oil reservoir plugging adjustment feasibility analysis method and system | |
| CN112943230A (en) | Residual oil distribution prediction method for common heavy oil reservoir | |
| CN110439532B (en) | Gas well oil testing conclusion dividing method | |
| CN115935843A (en) | Water injection induced dynamic fracture seepage numerical simulation method considering seepage mechanism | |
| CN104392131A (en) | Calculating method of broken rock seepage field in process of driving sand by water | |
| CN114757029A (en) | Method and system for simulating multi-stage emission reduction filling construction of alpha-beta wave of offshore long horizontal well | |
| CN109707352B (en) | Experimental device and experimental method for measuring nitrogen and nitrogen foam assisted gravity oil displacement efficiency |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CB03 | Change of inventor or designer information | ||
| CB03 | Change of inventor or designer information |
Inventor after: Chen Wei Inventor after: Duan Yonggang Inventor after: Fang Quantang Inventor before: Chen Wei Inventor before: Xu Yandong Inventor before: Duan Yonggang Inventor before: Fang Quantang Inventor before: Chen Lijie Inventor before: Zhou Ning Inventor before: Liu Qingshan Inventor before: Wang Qincong Inventor before: Jiang Yingbing |