CN104504472A - Method and device for predicting well productivity of gas storage gas-water interaction area - Google Patents
Method and device for predicting well productivity of gas storage gas-water interaction area Download PDFInfo
- Publication number
- CN104504472A CN104504472A CN201410809184.6A CN201410809184A CN104504472A CN 104504472 A CN104504472 A CN 104504472A CN 201410809184 A CN201410809184 A CN 201410809184A CN 104504472 A CN104504472 A CN 104504472A
- Authority
- CN
- China
- Prior art keywords
- gas
- well
- air water
- gas storage
- reservoir
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 201
- 238000003860 storage Methods 0.000 title claims abstract description 146
- 238000000034 method Methods 0.000 title claims abstract description 103
- 230000003993 interaction Effects 0.000 title abstract 6
- 238000004519 manufacturing process Methods 0.000 claims abstract description 92
- 238000006073 displacement reaction Methods 0.000 claims abstract description 83
- 230000035699 permeability Effects 0.000 claims abstract description 76
- 230000008569 process Effects 0.000 claims abstract description 68
- 238000011160 research Methods 0.000 claims abstract description 47
- 238000002474 experimental method Methods 0.000 claims abstract description 38
- 238000011161 development Methods 0.000 claims abstract description 20
- 238000012360 testing method Methods 0.000 claims abstract description 20
- 238000013461 design Methods 0.000 claims abstract description 18
- 230000008859 change Effects 0.000 claims abstract description 10
- 230000001603 reducing effect Effects 0.000 claims abstract description 6
- 230000002452 interceptive effect Effects 0.000 claims description 47
- 238000005325 percolation Methods 0.000 claims description 35
- 230000004941 influx Effects 0.000 claims description 33
- 239000004576 sand Substances 0.000 claims description 24
- 230000006870 function Effects 0.000 claims description 20
- 238000009933 burial Methods 0.000 claims description 12
- 238000004088 simulation Methods 0.000 claims description 7
- 230000006835 compression Effects 0.000 claims description 6
- 238000007906 compression Methods 0.000 claims description 6
- 238000012407 engineering method Methods 0.000 claims description 6
- JEGUKCSWCFPDGT-UHFFFAOYSA-N h2o hydrate Chemical compound O.O JEGUKCSWCFPDGT-UHFFFAOYSA-N 0.000 claims description 6
- 238000013508 migration Methods 0.000 claims description 6
- 230000005012 migration Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 abstract description 6
- 239000007789 gas Substances 0.000 description 220
- 238000010586 diagram Methods 0.000 description 11
- 238000004590 computer program Methods 0.000 description 7
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 6
- 230000008901 benefit Effects 0.000 description 3
- 238000004364 calculation method Methods 0.000 description 3
- 239000003345 natural gas Substances 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 230000001105 regulatory effect Effects 0.000 description 2
- 230000004888 barrier function Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 238000009533 lab test Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000005303 weighing Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Health & Medical Sciences (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Public Health (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Water Supply & Treatment (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method and a device for predicting the well productivity of a gas storage gas-water interaction area, wherein the method comprises the following steps: determining a functional relation between gas reservoir water intrusion and pressure, and determining a research object well in a gas-water interaction area in the operation process of rebuilding a gas storage according to the functional relation; processing the gas reservoir development data or well testing data of the research object well to obtain a gas well binomial productivity equation, and determining the influence factors of gas-water interaction displacement for reducing productivity; designing a seepage experiment in the gas storage building process according to the influence factors of the gas-water interaction displacement for reducing the productivity, and testing the gas phase permeability change condition of each cycle period of the design life of the gas storage; and correcting a gas well binomial energy production equation by using the seepage experimental result to obtain the capacity of the research object well, which deducts the seepage capability lost due to gas-water interaction displacement in the operation process of the gas storage. The method can improve the accuracy of the well productivity prediction of the gas storage constructed by the flooded gas reservoir.
Description
Technical field
The present invention relates to natural gas storage, gas reservoir engineering, laboratory experiment technical field, particularly relate to Forecasting Methodology and the device of gas storage air water interactive areas well productivity.
Background technology
Gas storage is the powerful measure regulating the rock gas equilibrium of supply and demand, and the important indicator weighing its regulating power is gas production ability, and gas production ability is realized by gas storage well production capacity, and Accurate Prediction well productivity is the technical barrier of puzzlement related scientific research personnel.
At present, during prediction watered-out gas reservoir gas storage reconstruction well productivity, mostly do not consider that air water hands over displacement effect on the impact of production capacity.As the article " the operational factor design of Xiangguo Monastery CARBONIFEROUS GAS RESERVOIR reconstruction underground natural gas storage " of document volume the 2nd February the 32nd in 2012 phase " gas industry ", directly the well productivity in gas reservoir development stage is equal to the well productivity after gas storage reconstruction; The article of document volume the 1st phase February the 32nd in 2011 " Xinjiang oil geology " determination of 21 pieces of underground natural gas storage gas well deliverabilities " Yong'an oilfield forever ", described method is only applicable to the situation of reservoir once water enchroachment (invasion) experience, cannot the physical process of air water repeatedly mutual displacement in accurate description gas storage operational process.
Therefore, there is limitation in the method for calculating productivity of current watered-out gas reservoir gas storage reconstruction, and the method for calculating productivity report being particularly positioned at air water interactive areas reservoir is less, needs relatively accurate computing method to occur.
Summary of the invention
The embodiment of the present invention provides the Forecasting Methodology of a kind of gas storage air water interactive areas well productivity, and in order to improve the accuracy of the well productivity prediction of watered-out gas reservoir gas storage reconstruction, the method comprises:
Determine the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determine according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Process gas reservoir development data or the well test data of described research object well, obtain gas well binomial potential curve and equation, determine that the mutual displacement of air water reduces the influence factor of production capacity;
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
With Seepage Experiment modified result gas well binomial potential curve and equation, obtain deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
In an embodiment, determine the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determine according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process, comprising:
Obtain clean water influx rate W by gas reservoir engineering method or method for numerical simulation, return with reservoir mean pressure P and have functional relation:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then determine that this well is described research object well.
In an embodiment, process gas reservoir development data or the well test data of described research object well, obtain gas well binomial potential curve and equation, determine that the mutual displacement of air water reduces the influence factor of production capacity, comprising:
Select gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to described two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P
e'-P
wf′＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P
e' be the pseudopressure of reservoir pressure, P
wf' be the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Determine that the influence factor that the mutual displacement of air water reduces production capacity is reservoir permeability K.
In an embodiment, reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, and the gas phase permeability situation of change of test gas storage each cycle period designed life, comprising:
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
Until K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
In an embodiment, with Seepage Experiment modified result gas well binomial potential curve and equation, obtain deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water, comprising:
Seepage Experiment result is substituted into the pseudopressure function of real gas, obtains the pseudopressure function expression revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation:
P
en″-P
ewf″＝AQ+BQ
2；
Wherein, P
en" be the modified quasi pressure of reservoir pressure, P
ewf" be the modified quasi pressure of sand face pressure.。
The embodiment of the present invention also provides the prediction unit of a kind of gas storage air water interactive areas well productivity, and in order to improve the accuracy of the well productivity prediction of watered-out gas reservoir gas storage reconstruction, this device comprises:
Research object well determination module, for determining the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determines according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Influence factor determination module, for the treatment of gas reservoir development data or the well test data of described research object well, obtains gas well binomial potential curve and equation, determines that the mutual displacement of air water reduces the influence factor of production capacity;
Seepage Experiment module, for reducing the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
Capability forecasting module, for Seepage Experiment modified result gas well binomial potential curve and equation, obtains deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
In an embodiment, research object well determination module specifically for:
Obtain clean water influx rate W by gas reservoir engineering method or method for numerical simulation, return with reservoir mean pressure P and have functional relation:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then determine that this well is described research object well.
In an embodiment, influence factor determination module specifically for:
Select gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to described two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P
e'-P
wf′＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P
e' be the pseudopressure of reservoir pressure, P
wf' be the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Determine that the influence factor that the mutual displacement of air water reduces production capacity is reservoir permeability K.
In an embodiment, Seepage Experiment module specifically for:
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
Until K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
In an embodiment, capability forecasting module specifically for:
Seepage Experiment result is substituted into the pseudopressure function of real gas, obtains the pseudopressure function expression revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation:
P
en″-P
ewf″＝AQ+BQ
2；
P
en" be the modified quasi pressure of reservoir pressure, P
ewf" be the modified quasi pressure of sand face pressure.
The embodiment of the present invention is when predicting air water interactive areas well productivity, not directly do not get as existing method the constant value that gas reservoir development phase process obtains, but deduct the production capacity of losing due to the mutual displacement of air water, and then give the production history of each cycle period in the gas storage life-span, compared with the conventional method, the result of calculation accuracy of the embodiment of the present invention increases substantially, and along with the prolongation of the gas storage cycle of operation, and technical advantage is entered and highlighted all the more.The prediction scheme of the embodiment of the present invention, provides relatively accurate well productivity technical parameter by for gas storage evaluation in early stage.
Accompanying drawing explanation
In order to be illustrated more clearly in the embodiment of the present invention or technical scheme of the prior art, be briefly described to the accompanying drawing used required in embodiment or description of the prior art below, apparently, accompanying drawing in the following describes is only some embodiments of the present invention, for those of ordinary skill in the art, under the prerequisite not paying creative work, other accompanying drawing can also be obtained according to these accompanying drawings.In the accompanying drawings:
Fig. 1 is the schematic diagram of the Forecasting Methodology of gas storage air water interactive areas well productivity in the embodiment of the present invention;
Fig. 2 is the schematic diagram of position, air water interactive areas in the embodiment of the present invention;
Fig. 3 is the schematic diagram of the prediction unit of gas storage air water interactive areas well productivity in the embodiment of the present invention.
Embodiment
For making the object of the embodiment of the present invention, technical scheme and advantage clearly understand, below in conjunction with accompanying drawing, the embodiment of the present invention is described in further details.At this, schematic description and description of the present invention is for explaining the present invention, but not as a limitation of the invention.
Inventor finds, existing exist circumscribed reason for the method for calculating productivity of air water interactive areas during watered-out gas reservoir gas storage reconstruction and be, not deduct in gas storage operational process the mutual displacement of air water to the reducing effect of production capacity, more can not predict the well productivity of gas storage one by one after cycle deduction air water mutual displacement impact, based on this, for the deficiency of existing air water interactive areas method of calculating productivity, deduct the percolation ability lost due to the mutual displacement of air water in gas storage operational process in embodiments of the present invention, and then obtain the method for a kind of relatively accurate prediction gas storage air water interactive areas well productivity.
Concrete, the embodiment of the present invention is according to the gas-water interface built in the reservoir of storehouse, determine the mutual reservoir area of the air water at gas storage well place, for the feature design physical simulation experiment of air water interactive areas, utilize experimental result to describe the percolation ability lost due to the mutual displacement of air water, finally obtain the consecutive variations Forecasting Methodology of air water interactive areas well productivity.
Fig. 1 is the schematic diagram of the Forecasting Methodology of gas storage air water interactive areas well productivity in the embodiment of the present invention.As shown in Figure 1, in the embodiment of the present invention, the Forecasting Methodology of gas storage air water interactive areas well productivity can comprise:
Step 101, determine the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determine according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Step 102, the gas reservoir development data processing described research object well or well test data, obtain gas well binomial potential curve and equation, determines that the mutual displacement of air water reduces the influence factor of production capacity;
Step 103, reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
Step 104, use Seepage Experiment modified result gas well binomial potential curve and equation, obtain deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
During concrete enforcement, first study the funtcional relationship of the clean water influx rate-pressure of gas reservoir, determine that certain well is in air water interactive areas in gas storage reconstruction operational process.
In embodiment, clean water influx rate W can be obtained by gas reservoir engineering method or method for numerical simulation, returns have functional relation with reservoir mean pressure P:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then this well is suitable for predicting production capacity in the embodiment of the present invention, determines that this well is described research object well.
After determining research object, process gas reservoir development data or well test data, obtain gas well binomial potential curve and equation, specifies the influence factor that the mutual displacement of air water reduces production capacity.In embodiment, can preferred gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to these two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P
e'-P
wf′＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P
e' be the pseudopressure of reservoir pressure, P
wf' be the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Therefore, the parameter affecting production history is mainly reservoir permeability K, needs the Changing Pattern of Study In Reservoir after gas storage reconstruction emphatically.
After determining that the mutual displacement of air water reduces the influence factor of production capacity, the influence factor of production capacity is reduced according to the mutual displacement of described air water, the i.e. feature of the mutual displacement of air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of each cycle period in the test gas storage n designed life cycle.
Fig. 2 is the schematic diagram of position, air water interactive areas in the embodiment of the present invention, as shown in Figure 2, air water interactive areas, reservoir between the gas-water interface that gas storage pressure upper limit is corresponding with threshold pression, namely gas storage fills Qi Shiweiqi district, gas production latter stage is pool, is in the mutual displacement state of air water in gas storage operational process.
According to the feature of the mutual displacement of air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
……
K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
Finally, with Seepage Experiment modified result gas well binomial potential curve and equation, and then obtain deducting the well productivity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
In embodiment, Seepage Experiment result can be substituted into the pseudopressure function of real gas, obtain the pseudopressure function expression revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation:
P
en″-P
ewf″＝AQ+BQ
2；
P
en" be the modified quasi pressure of reservoir pressure, P
ewf" be the modified quasi pressure of sand face pressure.
Can learn, A and B in deliverability equation remains unchanged before and after revising, and then ensure that remaining valid of A, B coefficient that the gas reservoir development stage calculates.
Visible, result by experiment, revise the convenient pseudopressure calculated of deliverability equation left end, A, B coefficient keeping equation right-hand member not easily to calculate is constant, can realize the correction of deliverability equation, obtains the expression formula of well productivity after the mutual displacement of any air water injection-production cycle.
Based on same inventive concept, additionally provide the prediction unit of a kind of gas storage air water interactive areas well productivity in the embodiment of the present invention, as described in the following examples.The principle of dealing with problems due to this device is similar to the Forecasting Methodology of gas storage air water interactive areas well productivity, and therefore the enforcement of this device see the enforcement of the Forecasting Methodology of gas storage air water interactive areas well productivity, can repeat part and repeat no more.
Fig. 3 is the schematic diagram of the prediction unit of gas storage air water interactive areas well productivity in the embodiment of the present invention.As shown in Figure 3, in the embodiment of the present invention, the prediction unit of gas storage air water interactive areas well productivity can comprise:
Research object well determination module 301, for determining the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determines according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Influence factor determination module 302, for the treatment of gas reservoir development data or the well test data of described research object well, obtains gas well binomial potential curve and equation, determines that the mutual displacement of air water reduces the influence factor of production capacity;
Seepage Experiment module 303, for reducing the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
Capability forecasting module 304, for Seepage Experiment modified result gas well binomial potential curve and equation, obtains deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
In embodiment, research object well determination module 301 specifically may be used for:
Obtain clean water influx rate W by gas reservoir engineering method or method for numerical simulation, return with reservoir mean pressure P and have functional relation:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then determine that this well is described research object well.
In embodiment, influence factor determination module 302 specifically may be used for:
Select gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to described two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P
e'-P
wf′＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P
e' be the pseudopressure of reservoir pressure, P
wf' be the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Determine that the influence factor that the mutual displacement of air water reduces production capacity is reservoir permeability K.
In embodiment, Seepage Experiment module 303 specifically may be used for:
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
Until K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
In embodiment, capability forecasting module 304 specifically may be used for:
Seepage Experiment result is substituted into the pseudopressure function of real gas, obtains the pseudopressure function expression revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation:
P
en″-P
ewf″＝AQ+BQ
2；
Wherein, P
en" be the modified quasi pressure of reservoir pressure, P
ewf" be the modified quasi pressure of sand face pressure.
In sum, the embodiment of the present invention is when predicting air water interactive areas well productivity, not directly do not get as existing method the constant value that gas reservoir development phase process obtains, but deduct the production capacity of losing due to the mutual displacement of air water, and then give the production history of each cycle period in the gas storage life-span, compared with the conventional method, the result of calculation accuracy of the embodiment of the present invention increases substantially, and along with the prolongation of the gas storage cycle of operation, technical advantage is entered and is highlighted all the more.The prediction scheme of the embodiment of the present invention, provides relatively accurate well productivity technical parameter by for gas storage evaluation in early stage.
Those skilled in the art should understand, embodiments of the invention can be provided as method, system or computer program.Therefore, the present invention can adopt the form of complete hardware embodiment, completely software implementation or the embodiment in conjunction with software and hardware aspect.And the present invention can adopt in one or more form wherein including the upper computer program implemented of computer-usable storage medium (including but not limited to magnetic disk memory, CD-ROM, optical memory etc.) of computer usable program code.
The present invention describes with reference to according to the process flow diagram of the method for the embodiment of the present invention, equipment (system) and computer program and/or block scheme.Should understand can by the combination of the flow process in each flow process in computer program instructions realization flow figure and/or block scheme and/or square frame and process flow diagram and/or block scheme and/or square frame.These computer program instructions can being provided to the processor of multi-purpose computer, special purpose computer, Embedded Processor or other programmable data processing device to produce a machine, making the instruction performed by the processor of computing machine or other programmable data processing device produce device for realizing the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.
These computer program instructions also can be stored in can in the computer-readable memory that works in a specific way of vectoring computer or other programmable data processing device, the instruction making to be stored in this computer-readable memory produces the manufacture comprising command device, and this command device realizes the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.
These computer program instructions also can be loaded in computing machine or other programmable data processing device, make on computing machine or other programmable devices, to perform sequence of operations step to produce computer implemented process, thus the instruction performed on computing machine or other programmable devices is provided for the step realizing the function of specifying in process flow diagram flow process or multiple flow process and/or block scheme square frame or multiple square frame.
Above-described specific embodiment; object of the present invention, technical scheme and beneficial effect are further described; be understood that; the foregoing is only specific embodiments of the invention; the protection domain be not intended to limit the present invention; within the spirit and principles in the present invention all, any amendment made, equivalent replacement, improvement etc., all should be included within protection scope of the present invention.
Claims (10)
1. a Forecasting Methodology for gas storage air water interactive areas well productivity, is characterized in that, comprising:
Determine the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determine according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Process gas reservoir development data or the well test data of described research object well, obtain gas well binomial potential curve and equation, determine that the mutual displacement of air water reduces the influence factor of production capacity;
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
With Seepage Experiment modified result gas well binomial potential curve and equation, obtain deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
2. the method for claim 1, is characterized in that, determines the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determines the research object well being in air water interactive areas in gas storage reconstruction operational process, comprising according to described funtcional relationship:
Obtain clean water influx rate W by gas reservoir engineering method or method for numerical simulation, return with reservoir mean pressure P and have functional relation:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then determine that this well is described research object well.
3. method as claimed in claim 2, is characterized in that, process gas reservoir development data or the well test data of described research object well, obtain gas well binomial potential curve and equation, determines that the mutual displacement of air water reduces the influence factor of production capacity, comprising:
Select gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to described two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P′
e-P′
wf＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P '
efor the pseudopressure of reservoir pressure, P '
wffor the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Determine that the influence factor that the mutual displacement of air water reduces production capacity is reservoir permeability K.
4. method as claimed in claim 3, it is characterized in that, reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life, comprising:
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
Until K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
5. method as claimed in claim 4, is characterized in that, with Seepage Experiment modified result gas well binomial potential curve and equation, obtains deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water, comprising:
Seepage Experiment result is substituted into the pseudopressure function of real gas, obtains the pseudopressure expression formula revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation:
P″
en-P″
ewf＝AQ+BQ
2；
Wherein, P "
enfor the modified quasi pressure of reservoir pressure, P "
ewffor the modified quasi pressure of sand face pressure.
6. a prediction unit for gas storage air water interactive areas well productivity, is characterized in that, comprising:
Research object well determination module, for determining the funtcional relationship of the clean water influx rate of gas reservoir and pressure, determines according to described funtcional relationship the research object well being in air water interactive areas in gas storage reconstruction operational process;
Influence factor determination module, for the treatment of gas reservoir development data or the well test data of described research object well, obtains gas well binomial potential curve and equation, determines that the mutual displacement of air water reduces the influence factor of production capacity;
Seepage Experiment module, for reducing the influence factor of production capacity according to the mutual displacement of described air water, design gas storage builds the Seepage Experiment of storehouse process, the gas phase permeability situation of change of test gas storage each cycle period designed life;
Capability forecasting module, for Seepage Experiment modified result gas well binomial potential curve and equation, obtains deducting the described research object well capacity losing percolation ability in gas storage operational process due to the mutual displacement of air water.
7. device as claimed in claim 6, is characterized in that, research object well determination module specifically for:
Obtain clean water influx rate W by gas reservoir engineering method or method for numerical simulation, return with reservoir mean pressure P and have functional relation:
W＝f(P)；
Substitute into gas storage pressure upper limit P
max, obtain water influx rate when gas storage moves to pressure upper limit:
W
max＝f(P
max)；
Substitute into gas storage threshold pression P
min, obtain water influx rate when gas storage moves to threshold pression:
W
min＝f(P
min)；
According to the depth of burial residing for research object well, calculate the migration of gas storage gas-water interface to pressure P during this depth of burial
well, and then obtain gas storage and move to pressure P
welltime water influx rate:
W
well=f (P
well);
If W
max<W
well<W
min, then determine that this well is described research object well.
8. device as claimed in claim 7, is characterized in that, influence factor determination module specifically for:
Select gas reservoir development stage two administration measure data points:
Production data point 1 (P
e1, P
wf1, Q
1);
Production data point 2 (P
e2, P
wf2, Q
2);
Wherein, P
e1the reservoir pressure of production data point 1, P
wf1the sand face pressure of production data point 1, Q
1it is the gas flow rate of production data point 1; P
e2the reservoir pressure of production data point 2, P
wf2the sand face pressure of production data point 2, Q
2it is the gas flow rate of production data point 2;
According to described two administration measure data points, process obtains pseudopressure gas well deliverability binomial equation:
P′
e-P′
wf＝AQ+BQ
2；
Wherein, P' is the pseudopressure of real gas:
P '
efor the pseudopressure of reservoir pressure, P '
wffor the pseudopressure of sand face pressure, Q is gas flow rate, and μ is gas viscosity, and Z is Gas Compression Factor, A and B is constant, by described two administration measure data point simultaneous solutions, is the function of permeability K;
Determine that the influence factor that the mutual displacement of air water reduces production capacity is reservoir permeability K.
9. device as claimed in claim 8, is characterized in that, Seepage Experiment module specifically for:
Reduce the influence factor of production capacity according to the mutual displacement of described air water, design gas drive water-water drive gas experiment, describes the physical process of the mutual displacement of air water interactive areas air water, obtains the gas phase relative permeability of the end points of air water interactive areas:
K
1/ K: gas phase relative permeability after the mutual displacement of first time air water, describes gas storage period 1 reservoir remaining relative permeability after a mutual district of air water causes percolation ability to reduce;
K
2/ K: gas phase relative permeability after the mutual displacement of second time air water, describes gas storage reservoir second round remaining relative permeability after the mutual district of secondary air water causes percolation ability to reduce;
K
3/ K: gas phase relative permeability after the mutual displacement of third time air water, describes gas storage period 3 reservoir remaining relative permeability after three mutual districts of air water cause percolation ability to reduce;
Until K
ngas phase relative permeability after/K: the n-th mutual displacement of air water, describes gas storage n-th cycle reservoir remaining relative permeability after n the mutual district of air water causes percolation ability to reduce.
10. device as claimed in claim 9, is characterized in that, capability forecasting module specifically for:
Seepage Experiment result is substituted into the pseudopressure function of real gas, obtains the pseudopressure function expression revised:
Revised pseudopressure expression formula is substituted into gas well binomial potential curve and equation, obtains revising rear binomial potential curve and equation: P "
en-P "
ewf=AQ+BQ
2;
P "
enfor the modified quasi pressure of reservoir pressure, P "
ewffor the modified quasi pressure of sand face pressure.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410809184.6A CN104504472B (en) | 2014-12-23 | 2014-12-23 | Method and device for predicting well productivity of gas storage gas-water interaction area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410809184.6A CN104504472B (en) | 2014-12-23 | 2014-12-23 | Method and device for predicting well productivity of gas storage gas-water interaction area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104504472A true CN104504472A (en) | 2015-04-08 |
| CN104504472B CN104504472B (en) | 2017-09-26 |
Family
ID=52945867
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410809184.6A Active CN104504472B (en) | 2014-12-23 | 2014-12-23 | Method and device for predicting well productivity of gas storage gas-water interaction area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104504472B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104964904A (en) * | 2015-06-16 | 2015-10-07 | 中国石油天然气股份有限公司 | Method and device for testing reservoir permeability of water flooded gas area of gas storage |
| CN107435528A (en) * | 2016-05-20 | 2017-12-05 | 中国石油天然气股份有限公司 | Volcanic gas reservoir gas well production allocation method |
| CN108629471A (en) * | 2017-03-20 | 2018-10-09 | 中国石油天然气股份有限公司 | Gas storage gas well interactive injection and production method |
| CN113250749A (en) * | 2021-06-22 | 2021-08-13 | 中国石油化工股份有限公司西北油田分公司 | Simulation method and system for condensate gas reservoir type gas storage |
| CN114756826A (en) * | 2022-06-16 | 2022-07-15 | 中海油田服务股份有限公司 | Gas well productivity prediction method and device, electronic equipment and storage medium |
| CN115828625A (en) * | 2022-12-21 | 2023-03-21 | 西南石油大学 | Gas storage well productivity calculation method considering reservoir salt deposition influence |
Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187961A (en) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | Output prediction system for emission hole complete well |
| CN101806224A (en) * | 2010-03-09 | 2010-08-18 | 煤炭科学研究总院重庆研究院 | Method for predicting extraction productivity of underground coal-seam gas |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Reservoir production capacity prediction method and device |
| FR2984562A1 (en) * | 2011-12-15 | 2013-06-21 | Terra 3E | METHOD AND SYSTEM FOR DYNAMICALLY MODELING POLYPHASE FLUID FLOW |
| CN103352693A (en) * | 2013-06-20 | 2013-10-16 | 中国石油天然气集团公司 | Productivity prediction system and method based on oil gas content |
| CN103577886A (en) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | Staged fracturing yield prediction method of low-permeability gas reservoir horizontal well |
| CN104134101A (en) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | Low-permeability reservoir natural gas productivity prediction method |
-
2014
- 2014-12-23 CN CN201410809184.6A patent/CN104504472B/en active Active
Patent Citations (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101187961A (en) * | 2007-11-20 | 2008-05-28 | 中国海洋石油总公司 | Output prediction system for emission hole complete well |
| CN101806224A (en) * | 2010-03-09 | 2010-08-18 | 煤炭科学研究总院重庆研究院 | Method for predicting extraction productivity of underground coal-seam gas |
| CN101906966A (en) * | 2010-07-16 | 2010-12-08 | 中国石油天然气股份有限公司 | Reservoir production capacity prediction method and device |
| FR2984562A1 (en) * | 2011-12-15 | 2013-06-21 | Terra 3E | METHOD AND SYSTEM FOR DYNAMICALLY MODELING POLYPHASE FLUID FLOW |
| CN103577886A (en) * | 2012-08-06 | 2014-02-12 | 中国石油化工股份有限公司 | Staged fracturing yield prediction method of low-permeability gas reservoir horizontal well |
| CN103352693A (en) * | 2013-06-20 | 2013-10-16 | 中国石油天然气集团公司 | Productivity prediction system and method based on oil gas content |
| CN104134101A (en) * | 2014-07-23 | 2014-11-05 | 中国石油集团川庆钻探工程有限公司 | Low-permeability reservoir natural gas productivity prediction method |
Non-Patent Citations (2)
| Title |
|---|
| 范新文等: ""涩北气田气井合理配产的综合技术对策"", 《天然气工业》 * |
| 郭平等: ""储气库多次气水互驱实验"", 《大庆石油地质与开发》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104964904A (en) * | 2015-06-16 | 2015-10-07 | 中国石油天然气股份有限公司 | Method and device for testing reservoir permeability of water flooded gas area of gas storage |
| CN107435528A (en) * | 2016-05-20 | 2017-12-05 | 中国石油天然气股份有限公司 | Volcanic gas reservoir gas well production allocation method |
| CN108629471A (en) * | 2017-03-20 | 2018-10-09 | 中国石油天然气股份有限公司 | Gas storage gas well interactive injection and production method |
| CN113250749A (en) * | 2021-06-22 | 2021-08-13 | 中国石油化工股份有限公司西北油田分公司 | Simulation method and system for condensate gas reservoir type gas storage |
| CN113250749B (en) * | 2021-06-22 | 2024-02-27 | 中国石油化工股份有限公司西北油田分公司 | Method and system for simulating condensate gas reservoir type gas storage |
| CN114756826A (en) * | 2022-06-16 | 2022-07-15 | 中海油田服务股份有限公司 | Gas well productivity prediction method and device, electronic equipment and storage medium |
| CN115828625A (en) * | 2022-12-21 | 2023-03-21 | 西南石油大学 | Gas storage well productivity calculation method considering reservoir salt deposition influence |
| CN115828625B (en) * | 2022-12-21 | 2023-11-21 | 西南石油大学 | Gas storage well productivity calculation method considering reservoir salt formation influence |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104504472B (en) | 2017-09-26 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Simons et al. | A model for overland flow and associated processes within the Hydroinformatics Modelling System | |
| CN104504472A (en) | Method and device for predicting well productivity of gas storage gas-water interaction area | |
| Sepehri et al. | Prediction of mining-induced surface subsidence and ground movements at a Canadian diamond mine using an elastoplastic finite element model | |
| CN105401927B (en) | Method and device for determining well pattern and well spacing of oil reservoir | |
| CN104636818A (en) | Method and device for predicting gas storage gas-oil interaction area well productivity | |
| Wu et al. | Consolidation analysis of Nuozhadu high earth-rockfill dam based on the coupling of seepage and stress-deformation physical state | |
| Tal et al. | An improved numerical manifold method for simulation of sequential excavation in fractured rocks | |
| Maihemuti et al. | Numerical simulation analysis of reservoir bank fractured rock-slope deformation and failure processes | |
| Wu et al. | Flow characteristics and escape-route optimization after water inrush in a backward-excavated karst tunnel | |
| CN104392109A (en) | Dynamic discrete fracture simulation method and device for oil reservoir | |
| CN103469780A (en) | Method of calculating deep antiskid stable critical slipping plane of dam foundation of gravity dam | |
| US10650175B2 (en) | Hydrocarbon migration and accumulation methods and systems | |
| CN104346511A (en) | Bran-new method and bran-new device for dynamically monitoring oil deposit | |
| Huang et al. | Simulating groundwater inflow in the underground tunnel with a coupled fracture-matrix model | |
| CN104680244A (en) | Oil and gas field productivity prediction method and device | |
| Wang et al. | Simulated crack and slip plane propagation in soil slopes with embedded discontinuities using XFEM | |
| Beer et al. | Isogeometric boundary element method for the simulation of underground excavations | |
| Deng et al. | 3D finite element modeling of directional hydraulic fracturing based on deformation reinforcement theory | |
| CN103942430A (en) | Building settlement prediction method based on combined model | |
| Fomychov et al. | Determination of parameters of non-uniform fractured rock massif in computing experiment | |
| Han et al. | Analytic solutions of the forces and displacements for multicentre circular arc tunnels | |
| Jiang et al. | Investigation of coordinated development of coal mining and pipeline protection under boundary effect of thick unconsolidated layer | |
| Zhang et al. | Research on the limit operating pressure of a bedded salt gas storage cavern: A case study | |
| Zhang et al. | Distribution for hydraulic head on tunnel structures in water-rich mountainous region considering influences of fault geology using virtual image technique | |
| Liu et al. | Analysis of external water pressure for a tunnel in fractured rocks |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |