CN109268004B - Shale gas reservoir medium structure coupling and seam network state identification method - Google Patents

Shale gas reservoir medium structure coupling and seam network state identification method Download PDF

Info

Publication number
CN109268004B
CN109268004B CN201811110598.4A CN201811110598A CN109268004B CN 109268004 B CN109268004 B CN 109268004B CN 201811110598 A CN201811110598 A CN 201811110598A CN 109268004 B CN109268004 B CN 109268004B
Authority
CN
China
Prior art keywords
stress
shale
seam
net
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.)
Active
Application number
CN201811110598.4A
Other languages
Chinese (zh)
Other versions
CN109268004A (en
Inventor
朱维耀
宋智勇
岳明
刘昀枫
陈震
马东旭
刘文超
亓倩
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology Beijing USTB
Original Assignee
University of Science and Technology Beijing USTB
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Science and Technology Beijing USTB filed Critical University of Science and Technology Beijing USTB
Priority to CN201811110598.4A priority Critical patent/CN109268004B/en
Publication of CN109268004A publication Critical patent/CN109268004A/en
Application granted granted Critical
Publication of CN109268004B publication Critical patent/CN109268004B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing 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
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N15/00Investigating characteristics of particles; Investigating permeability, pore-volume or surface-area of porous materials
    • G01N15/08Investigating permeability, pore-volume, or surface area of porous materials

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Chemical & Material Sciences (AREA)
  • Dispersion Chemistry (AREA)
  • Health & Medical Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Geophysics (AREA)
  • Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)

Abstract

The invention provides a shale gas reservoir medium structure coupling and seam network shape identification method, and belongs to the technical field of oil and natural gas exploitation. The method comprises the steps of firstly, carrying out a fracturing test on a shale core by utilizing an improved Brazilian fracturing test, and classifying the form of a fracture network; secondly, performing shale matrix-fracture structure and stress sensitivity test on rock cores with different fracture forms; thirdly, solving stress sensitivity constant ranges of different seam net forms according to indoor core data; and finally, drawing a seam net shape judging and identifying plate based on an indoor rock core stress sensitivity experiment and in combination with production practice. The chart can be applied to actual fracturing wells, the fracture morphology is directly judged and identified through real-time effective stress and normalized flow, and the chart has profound significance on implementation effect of shale gas reservoir volume fracturing and establishment and development work system.

Description

Shale gas reservoir medium structure coupling and seam network state identification method
Technical Field
The invention relates to the technical field of oil and natural gas exploitation, in particular to a shale gas reservoir medium structure coupling and seam network shape identification method.
Background
In 2017, the dependence of natural gas in China on the outside is as high as 39%, nearly 50% is estimated in 2035, and the country faces a severe oil and gas safety problem. The shale gas reservoir resources in China are rich, and the amount of the resources which can be collected by the shale gas in China is about 36 billion cubic meters and about 20 percent of the whole world in 2017. Compared with foreign shale gas reservoirs, the shale gas reservoir in China has the characteristics of deep burial depth (generally 1500-4000 meters and maximally 6950 meters), low permeability (the permeability is 0.01-1 mD), micro-crack development and complex multi-scale flow mechanism.
Reservoir fracturing modification is a main mode for improving the development effect of shale gas reservoirs, and enables a multi-scale pore network consisting of nano-micron pores, micro fractures and artificial fractures to be formed in the shale reservoirs, so that strong nonlinear flow characteristics exist in the shale gas reservoirs. In the process of exploiting the shale reservoir, the structures of nano-micron pores, micro-cracks and artificial cracks are coupled by a seepage field and a stress field, the mechanism is complex, the nonlinearity is strong, and no technology and method for observing and monitoring the crack morphology exist at present, so that a method for determining the shale gas reservoir multi-medium coupling and the crack network morphology is urgently needed to be established.
Disclosure of Invention
The invention aims to provide a shale gas reservoir medium structure coupling and seam network state identification method.
The method comprises the following steps:
(1) classifying the shapes of the sewing nets: carrying out a fracturing test on the shale core by using an improved Brazilian fracturing test, and dividing the fracture morphology into: a tree-shaped seam net, a feather-shaped seam net, a cluster-shaped seam net and a net-shaped seam net; the improved Brazilian fracturing experiment is a method described in the invention patent of a fracture propagation analysis device and an analysis method of a shale core (CN 201410548672.6);
(2) carrying out shale matrix-fracture structure and stress sensitivity experiment tests: selecting the rock cores with different fracture forms in the step (1), carrying out shale matrix-fracture structure and stress sensitivity experiment tests through a rock core displacement experiment (setting the effective stress range to be 4-22 MPa), and recording rock core permeability and effective stress parameters;
(3) and (3) solving stress sensitivity constants of different seam net forms: fitting experimental data in an exponential function mode according to the permeability and effective stress parameters of the rock cores with different fracture forms under different stress conditions in the step (2), and solving a stress sensitivity constant;
(4) drawing a sewing net shape judgment picture version: and (4) establishing a model of effective stress and production well flow according to the actually measured stress sensitivity constants of the four seam network states obtained in the step (3), and carrying out normalization processing on flow parameters to finally form a seam network state identification chart.
Wherein, the fitting of the experimental data in the form of exponential function in the step (3) is specifically as follows:
according to the method with the effect, the mathematical expression of the change rate of the effective stress and the permeability is as follows:
Figure BDA0001808994320000021
wherein sigma is overburden formation pressure in MPa, p is reservoir pore pressure in MPa, α is effective stress coefficient, and K is0Rock permeability at effective stress zero, in units of 10-3μm2(ii) a K is the permeability under any formation pressure condition and has the unit of 10-3μm2(ii) a b is stress sensitive constant in Mpa-1Wherein the shale reservoir has a microcrack development characteristic, and α is 1.
The stress sensitivity constant b is the basis for distinguishing the shapes of the seam networks, and since K and Q (gas production in the same day) are positively correlated in the actual production process, namely K/K0 is Q/Q0, the exponential relation between the normalized flow and the effective stress can be obtained; therefore, the process of drawing the slot-net shape judging plate in the step (4) comprises the following steps: and (3) drawing stress sensitivity charts of different slit net shapes based on the indoor core stress sensitivity experiment in the step (2), wherein the abscissa is effective stress, and the ordinate is a normalized flow parameter.
Wherein the normalized flow parameters are the daily gas production rate Q and the daily gas production rate peak value QmaxThe ratio of (a) to (b).
The technical scheme of the invention has the following beneficial effects:
in the scheme, the formed chart can be applied to an actual fracturing well, the fracture morphology is directly judged and identified through real-time effective stress and normalized flow, and the method has profound significance on implementation effect of shale gas reservoir volume fracturing and establishment and development work system.
Drawings
FIG. 1 is a flow chart of a shale gas reservoir medium structure coupling and seam network shape identification method of the present invention;
FIG. 2 is a schematic view of different seam net configurations in the present invention, wherein (a) is a feathered seam net, (b) is a tufted seam net, (c) is a net seam net, and (d) is a tree seam net;
fig. 3 is a stress sensitivity characteristic curve of different slot net shapes.
Detailed Description
In order to make the technical problems, technical solutions and advantages of the present invention more apparent, the following detailed description is given with reference to the accompanying drawings and specific embodiments.
The invention provides a shale gas reservoir medium structure coupling and seam network state identification method.
As shown in fig. 1, the method comprises the following steps:
(1) classifying the shapes of the sewing nets: carrying out a fracturing test on the shale core by using an improved Brazilian fracturing test, and dividing the fracture morphology into: a tree-shaped seam net, a feather-shaped seam net, a cluster-shaped seam net and a net-shaped seam net;
(2) carrying out shale matrix-fracture structure and stress sensitivity experiment tests: selecting the rock cores with different fracture forms in the step (1), carrying out shale matrix-fracture structure and stress sensitivity experiment tests through a rock core displacement experiment, and recording rock core permeability and effective stress parameters;
(3) and (3) solving stress sensitivity constants of different seam net forms: fitting experimental data in an exponential function mode according to the permeability and effective stress parameters of the rock cores with different fracture forms under different stress conditions in the step (2), and solving a stress sensitivity constant;
(4) drawing a sewing net shape judgment picture version: and (4) establishing a model of effective stress and production well flow according to the actually measured stress sensitivity constants of the four seam network states obtained in the step (3), and carrying out normalization processing on flow parameters to finally form a seam network state identification chart.
In the specific design, the method is based on the rock core experiment in a certain shale gas reservoir in China, summarizes the shale gas reservoir medium structure coupling and seam net shape identification method, can quantitatively represent the complex seam net formed after the shale gas reservoir is subjected to volume fracturing, and has great significance for on-site actual development and recovery ratio prediction.
The specific implementation mode is as follows:
(1) shale matrix-fracture structure and stress sensitivity experimental test
Matrix shale stress sensitivity test experiments were conducted at room temperature and pressure using nitrogen as the experimental gas to simulate natural gas. Confining pressure equipment uses a high precision plunger displacement pump. The back pressure control system uses a BP-100 air spring back pressure valve produced by American company and uses a high-precision multi-stage plunger displacement pressure pump for control. The adopted experimental method is a differential pressure-flow method, and the effective stresses selected in the experiment are respectively 4MPa, 5MPa, 7MPa, 9MPa, 11MPa, 13MPa, 16MPa, 19MPa and 22 MPa. The experiment is set with a back pressure of 1MPa and an inlet pressure of 3MPa and kept constant.
The experimental steps are as follows:
① drying the core in an oven at a constant temperature of 70 deg.C for 50h, and measuring basic data such as length, diameter, weight, porosity and permeability;
② putting the rock core into a rock core holder, zeroing the initial value of the instrument, adding confining pressure to 4MPa, and adding back pressure to 1MPa to be unchanged;
③ in the whole experiment process, the injection pressure is kept constant at 3MPa, and a gas single-phase seepage experiment is carried out;
④ recording the permeability of rock sample under initial effective stress when the seepage state is stable, adjusting confining pressure according to the pre-designed effective stress value, recording different permeability K in the confining pressure increasing process, determining different confining pressure effective stress values according to the plan, after reaching the maximum effective stress, gradually reducing pressure according to the pressure point drawn up by the boosting experiment, and determining the permeability, and ending the experiment.
The permeability of the manually fractured rock sample is high, and the experimental error is large by adopting a pressure difference-flow method, so that a covering pressure hole permeability instrument KFSY/T08-055 is adopted for testing, the experimental fluid adopts high-purity nitrogen, and the selected rock core is dried in a constant temperature cabinet for 48 hours at 105 ℃. The pressure difference between the inlet and the outlet is 2MPa, and the pressure difference is measured by 4MPa, 7MPa, 14MPa, 21MPa and 27 MPa. The permeability is likewise dimensionless, i.e. expressed as K/K0,K0For initial permeability, see fig. 2 for details.
Experimental results and discussion:
the shale matrix-fracture structure and stress sensitivity experiment result shows that: the shale permeability has larger change amplitude, the initial permeability is lower, the decline amplitude is larger, which is mainly determined by the pore distribution characteristics of the rock sample, the average pore radius of the rock sample with lower permeability is smaller, when the effective stress changes, the small pores are easy to close, thereby leading to the decline of permeability, and the stress sensitivity phenomenon is easier to generate relative to the rock sample with high permeability.
And fitting the experimental curve data, wherein the effective stress and the permeability change rate have a good exponential relationship according to an effective method. The mathematical expression is as follows:
Figure BDA0001808994320000041
wherein sigma is overburden formation pressure in MPa, p is reservoir pore pressure in MPa, α is effective stress coefficient, and K is0Rock permeability at effective stress zero, in units of 10-3μm2(ii) a K is the permeability under any formation pressure condition and has the unit of 10-3μm2(ii) a b is stress sensitive constant in Mpa-1N.r. warlinski and l.w.teufel in 1992 experimentally given the effective stress coefficients of different porous media, consider a reservoir with fracture development, α → 1. shale reservoir has microfracture development characteristics, therefore α is approximately equal to 1.
Under the condition of a gas reservoir, under the influence of rock pressure of an overlying stratum, pores of the reservoir are in a compressed state, after the reservoir is separated from the gas reservoir condition, the pressure borne by a rock framework is relieved, and partial pore channels in the reservoir are opened or enlarged. The condition of a real reservoir stratum cannot be reflected in a low effective stress region in an experiment, and if the reservoir stratum is evaluated, the stress sensitivity of the shale gas reservoir is tested by taking the initial effective pressure of an original stratum as a starting point. The reservoir burial depth of the target area is about 1600m, the overlying formation pressure and the pore pressure of the shale reservoir are respectively about 34MPa and 25MPa according to the information of field formation testing, well logging and the like, and therefore the effective pressure of 9MPa is selected as the effective stress starting point of the reservoir to be evaluated.
All experimental points of the shale rock sample are summed up to be storedThe permeability with the layer effective stress as a starting point is fitted with the effective stress through an exponential function, and an exponential term coefficient in the formula is a stress sensitive constant. The stress sensitivity constant is a parameter reflecting the deformation degree of the porous medium along with stress, and is related to the physical properties of the rock and the fracture network form. The stress sensitivity constant is 0.101MPa calculated based on all experimental points-1~0.322MPa-1(ii) a The stress sensitivity constant ranges for different seam net forms are shown in table 1.
TABLE 1 corresponding table of different seam net forms and stress sensitivity constants
Figure BDA0001808994320000051
The analysis is carried out from four different slot-net shapes, namely feather-shaped slots, cluster-shaped slots, net-shaped slots and tree-shaped slots, and is shown in detail in figure 3. The chart can be applied to actual fracturing wells, and the fracture form can be directly judged and identified through real-time effective stress and normalized flow.
Application example:
aiming at a horizontal well of a shale gas reservoir in Chongqing, the burial depth is 1600m, the initial pressure of a stratum is 35MPa, and the daily gas production at the initial stage of volume fracturing is 20 ten thousand square. At present, the daily gas production is 2 ten thousand square, and the bottom hole flowing pressure is 15 MPa.
The effective stress of the stratum rock is 20MPa according to the conversion of the original stratum pressure and the bottom hole flowing pressure, and the daily gas production is reduced by 90 percent at present, namely Q/Qmax0.1, the stress sensitivity constant is 0.1151. The stress sensitivity constant range table is determined to be a tree-shaped seam network, the conclusion can be directly judged that the fracture network is in the shape of the tree-shaped seam network, and the stress sensitivity constant does not need to be calculated again.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention as defined in the appended claims.

Claims (3)

1. A shale gas reservoir medium structure coupling and seam network state identification method is characterized by comprising the following steps: the method comprises the following steps:
(1) classifying the shapes of the sewing nets: carrying out a fracturing test on the shale core by using an improved Brazilian fracturing test, and dividing the fracture morphology into: a tree-shaped seam net, a feather-shaped seam net, a cluster-shaped seam net and a net-shaped seam net;
(2) carrying out shale matrix-fracture structure and stress sensitivity experiment tests: selecting the rock cores with different fracture forms in the step (1), carrying out shale matrix-fracture structure and stress sensitivity experiment tests through a rock core displacement experiment, and recording rock core permeability and effective stress parameters;
(3) and (3) solving stress sensitivity constants of different seam net forms: fitting experimental data in an exponential function mode according to the permeability and effective stress parameters of the rock cores with different fracture forms under different stress conditions in the step (2), and solving a stress sensitivity constant;
(4) drawing a sewing net shape judgment picture version: establishing a model of effective stress and production well flow according to the actually measured stress sensitivity constants of the four seam network states obtained in the step (3), and carrying out normalization processing on flow parameters to finally form a seam network state identification chart;
the process of drawing the slot net shape judging plate in the step (4) comprises the following steps: and (3) drawing stress sensitivity charts of different slit net shapes based on the indoor core stress sensitivity experiment in the step (2), wherein the abscissa is effective stress, and the ordinate is a normalized flow parameter.
2. The shale gas reservoir medium structure coupling and seam network form identification method according to claim 1, wherein: the fitting of the experimental data in the form of an exponential function in the step (3) is specifically as follows:
according to the method with the effect, the mathematical expression of the change rate of the effective stress and the permeability is as follows:
Figure FDA0002447568160000011
wherein sigma is overburden formation pressure in MPa, p is reservoir pore pressure in MPa, αIs the effective stress coefficient; k0Rock permeability at effective stress zero, in units of 10-3μm2(ii) a K is the permeability under any formation pressure condition and has the unit of 10-3μm2(ii) a b is stress sensitive constant in Mpa-1Wherein the shale reservoir has a microcrack development characteristic, and α is 1.
3. The shale gas reservoir medium structure coupling and seam network form identification method according to claim 1, wherein: the normalized flow parameter is the ratio of the daily gas production to the peak daily gas production.
CN201811110598.4A 2018-09-21 2018-09-21 Shale gas reservoir medium structure coupling and seam network state identification method Active CN109268004B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811110598.4A CN109268004B (en) 2018-09-21 2018-09-21 Shale gas reservoir medium structure coupling and seam network state identification method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811110598.4A CN109268004B (en) 2018-09-21 2018-09-21 Shale gas reservoir medium structure coupling and seam network state identification method

Publications (2)

Publication Number Publication Date
CN109268004A CN109268004A (en) 2019-01-25
CN109268004B true CN109268004B (en) 2020-06-26

Family

ID=65198289

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811110598.4A Active CN109268004B (en) 2018-09-21 2018-09-21 Shale gas reservoir medium structure coupling and seam network state identification method

Country Status (1)

Country Link
CN (1) CN109268004B (en)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110410052B (en) * 2019-07-22 2020-04-17 四川省科源工程技术测试中心 Device and method for simulating shale gas reservoir fracture drainage gas recovery efficiency
CN113655545B (en) * 2020-08-26 2022-08-02 中国石油大学(北京) Modeling method of stress-sensitive damaged oil-gas layer, damage degree spatial-temporal evolution 4D quantitative and intelligent diagnosis method and system thereof
CN112001134B (en) * 2020-09-04 2021-09-07 北京科技大学 Identification method and device for unconventional gas reservoir flow field structure
CN112417787B (en) 2020-11-26 2021-06-29 北京科技大学 Unconventional oil reservoir two-phase relative permeability curve measuring device and method
CN115749713B (en) * 2022-10-14 2023-06-16 中国矿业大学 Rock stratum variable frequency pulse seam net fracturing method and equipment

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105182423B (en) * 2015-10-22 2016-06-22 中国石油大学(华东) A kind of integrated recognition method in superpressure crack
CN105319603B (en) * 2015-11-06 2018-04-17 中国石油大学(华东) The Forecasting Methodology of tight sandstone reservoir complexity chicken-wire cracking
CN108280288B (en) * 2018-01-22 2021-07-16 中石化石油工程技术服务有限公司 Reservoir permeability model establishing method based on stress sensitivity analysis
CN108518212B (en) * 2018-04-09 2020-10-16 西南石油大学 Method for calculating unsteady state yield of shale gas reservoir complex fracture network

Also Published As

Publication number Publication date
CN109268004A (en) 2019-01-25

Similar Documents

Publication Publication Date Title
CN109268004B (en) Shale gas reservoir medium structure coupling and seam network state identification method
WO2020056750A1 (en) Method for identifying medium structure coupling and seam mesh morphology of shale gas reservoirs
AU2020100760A4 (en) Method for identifying medium structure coupling and fracture network morphology of shale gas reservoirs
CN109594968B (en) Shale gas multistage fracturing horizontal well post-fracturing fracture parameter evaluation method and system
CN107506948B (en) Shale oil gas comprehensive yield analysis method based on dynamic drainage volume
US8838427B2 (en) Method for determining the closure pressure of a hydraulic fracture
CN109946436B (en) Fracture-property dense gas layer working fluid damage evaluation method considering both foundation block and fracture
CN106909758B (en) A kind of method of fine and close oily reservoir-level well multistage sub-clustering perforating site optimization design
CN108166963B (en) Method for evaluating fracturing effect of offshore oil and gas well
Hawkins Laboratory study of proppant-pack permeability reduction caused by fracturing fluids concentrated during closure
CN111456709B (en) Horizontal well multistage fracturing segmentation clustering method based on logging curve
CN110656915B (en) Shale gas multi-section fracturing horizontal well multi-working-system productivity prediction method
CN113919240A (en) High-pressure gas reservoir parameter calculation method based on oil and gas well production data
CN113283108B (en) Method and system for quantitatively evaluating fracturing property of shale oil reservoir
CN112211627B (en) Selection method of low-permeability gas reservoir interference well testing test well
CN105930932A (en) Gas index-based shale-gas-layer standardized open-flow capacity obtaining method
CN110595953A (en) Experimental test device and method for shale mixing wettability
CN106640042A (en) Device and method for evaluating single well productivity of gas well
CN113640119A (en) Method for determining stress-related rock dynamic Biot coefficient
CN105678082B (en) Double-pressure drop method for identifying type of oil-gas well acid fracturing communication reservoir
CN206057147U (en) Shale matrix the dynamic damage evaluating apparatus based on fluid pressure pulse
CN105334293A (en) Flow state identification method and device of tight gas in multi mediums
CN115809536A (en) Evaluation method for multi-section fracturing reformation of shale gas well
CN110454152A (en) The non-homogeneous injury skin factor calculation method of fracture-type reservoir
CN111239372B (en) Carbonate rock pore structure classification method based on overburden seepage experiment

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