CN113049784A - Prediction method suitable for water sensitivity of shale reservoir - Google Patents

Prediction method suitable for water sensitivity of shale reservoir Download PDF

Info

Publication number
CN113049784A
CN113049784A CN202110312332.3A CN202110312332A CN113049784A CN 113049784 A CN113049784 A CN 113049784A CN 202110312332 A CN202110312332 A CN 202110312332A CN 113049784 A CN113049784 A CN 113049784A
Authority
CN
China
Prior art keywords
sample
illite
water
area
sensitive
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.)
Pending
Application number
CN202110312332.3A
Other languages
Chinese (zh)
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.)
Shaanxi Yanchang Petroleum Group Co Ltd
Original Assignee
Shaanxi Yanchang Petroleum Group Co Ltd
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 Shaanxi Yanchang Petroleum Group Co Ltd filed Critical Shaanxi Yanchang Petroleum Group Co Ltd
Priority to CN202110312332.3A priority Critical patent/CN113049784A/en
Publication of CN113049784A publication Critical patent/CN113049784A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N33/00Investigating or analysing materials by specific methods not covered by groups G01N1/00 - G01N31/00
    • G01N33/24Earth materials
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N1/00Sampling; Preparing specimens for investigation
    • G01N1/28Preparing specimens for investigation including physical details of (bio-)chemical methods covered elsewhere, e.g. G01N33/50, C12Q
    • 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
    • G01N15/088Investigating volume, surface area, size or distribution of pores; Porosimetry

Landscapes

  • Chemical & Material Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Physics & Mathematics (AREA)
  • Pathology (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • Food Science & Technology (AREA)
  • Medicinal Chemistry (AREA)
  • Dispersion Chemistry (AREA)
  • Silicates, Zeolites, And Molecular Sieves (AREA)

Abstract

The invention relates to a prediction method suitable for water sensitivity of a shale reservoir, which comprises the following steps: sample pretreatment; finding an illite and illite mixed layer in the cubic sample in a target view field; and confirming that illite and illite-montmorillonite mixed layers exist in the target vision field through an energy spectrum analysis system; measurement of illite content WI% and content of illite-montmorillonite mixed layer WISPercent; measuring illite mineral area S before water-sensitive injuryI frontIllite-smectite mixed layer mineral area SIS frontAnd pore area Sk front(ii) a Water-sensitive injuries; coating a film on the surface of the cube sample subjected to water-sensitive injury, and then placing the cube sample on the sample seat again at the same position; measuring illite mineral area S after water-sensitive injuryAfter IIllite-smectite mixed layer mineral area SAfter ISAnd pore area Sk is back(ii) a Calculating the illite mineral areaExpansion ratio phiIIllite-smectite mixed layer mineral area expansion rate phiISAnd rate of change of pore area phik(ii) a Then calculating a water sensitivity predicted value Y; the invention solves the problem that the surface porosity cannot be quantitatively analyzed due to the difficulty in die casting of the shale sample casting body slice.

Description

Prediction method suitable for water sensitivity of shale reservoir
Technical Field
The invention relates to the technical field of petroleum, in particular to a prediction method suitable for water sensitivity of a shale reservoir.
Background
During the operation processes of drilling, well completion, fracturing and injection production of oil and gas fields, external fluid enters the stratum and is easy to react with clay minerals in a series of ways, such as clay mineral expansion, migration, dispersion and the like, so that the influence on the permeability and the productivity of reservoir holes is brought to different degrees. Irreversible reservoir damage can also result if not properly constructed. The clay expansion mainly causes the reduction of the radius of a seepage pore or a throat of a reservoir, the migration of particles blocks seepage channel systems such as the pore and the throat, and the permeability of the reservoir is reduced.
The core displacement method is the most common reservoir water sensitivity evaluation method at present, and is also a macroscopic comprehensive water sensitivity evaluation method closest to reservoir conditions. Various fluids associated with formation damage are injected under approximate formation conditions according to darcy's law, and the permeability of the rock sample and its changes are measured to evaluate the extent of reservoir permeability damage. Reference is mainly made to the experimental industry standard (SY/T5358-2010), but the standard is mainly suitable for sensitivity evaluation of clastic rock reservoir rock samples with air permeability of more than 1 millidarcy. No reference standard exists for shale and tight sandstone reservoirs with air permeability less than 1 millidarcy.
Because the shale is very compact, the conventional core displacement experiment method is difficult to complete the water sensitivity test of the shale, and in the displacement process, fluid is difficult to enter a core sample, and other means such as manual seam making, use of an artificial core and the like are usually needed. The natural shale sample contains a large amount of clay minerals, is very easy to swell when meeting water, has very low success rate in the process of artificial seam making and sample drilling, and brings certain difficulty to the preparation of the sample.
Therefore, the water-sensitive damage degree of the shale reservoir stratum is indirectly predicted by combining tests such as a shale expansion instrument test method, a pore throat radius measurement test method, water-sensitive mineral crystal layer spacing measurement and the like.
Shale dilatometer experimental method: grinding a rock core sample to particles with a certain size, soaking the particles in liquid such as formation water, distilled water and the like, and testing the rock sample expansion rate by adopting a shale expansion instrument, wherein the calculation formula is as follows:
Figure BDA0002990325890000011
wherein E is the expansion ratio,%; h istThe height of the clay sample at the time t is mm; h is0Is the initial height, mm, of the clay sample.
Throat radius measurement test: and obtaining the pore structure and pore diameter change of the reservoir by a mercury intrusion method, a BET method and a nuclear magnetic resonance spectrum analysis method.
Testing the interval between crystal layers of the water-sensitive minerals: by using the X-ray diffraction principle, the inter-layer distance between the water-sensitive mineral and the crystal layer before and after expansion can be quantitatively measured according to the Bragg equation (d ═ lambda/2 sin theta).
The method is a water-sensitive injury evaluation method and has the following defects: firstly, mineral types and occurrence state changes which cause reservoir water sensitivity damage cannot be accurately and intuitively obtained; reservoir space structure change and pore area change of reservoir water sensitivity damage cannot be accurately and intuitively obtained; and thirdly, the water-sensitive mineral expansion rate change caused by the water-sensitive damage cannot be intuitively and accurately obtained.
Disclosure of Invention
The invention aims to solve the problems and provides a prediction method suitable for water sensitivity of a shale reservoir.
The technical scheme of the invention is as follows:
a prediction method suitable for water sensitivity of a shale reservoir comprises the following steps:
(1) sample pretreatment: making the core sample into a cubic sample of 10mm multiplied by 10 mm; fixing the cubic sample on a sample holder;
(2) selecting a target view, and finding the illite and the illite-montmorillonite mixed layer in the cubic sample in the target view; and confirming that illite and illite-montmorillonite mixed layers exist in the target vision field through an energy spectrum analysis system;
(3) measuring the content W of illite in the cubic sampleI% content of illite-montmorillonite mixed layer WIS%;
Measurement of illite mineral area S of cubic sample before water-sensitive injuryI frontIllite-smectite mixed layer mineral area SIS frontAnd pore area Sk front
(4) Taking down the cubic sample from the sample seat for water-sensitive damage;
(5) coating a film on the surface of the cube sample subjected to water-sensitive injury, and then placing the cube sample on the sample seat again at the same position; measurement of illite mineral area S of cube sample after water-sensitive injuryAfter IIllite-smectite mixed layer mineral area SAfter ISAnd pore area Sk is back(ii) a Calculating the area expansion rate phi of the illite mineralIIllite-smectite mixed layer mineral area expansion rate phiISAnd rate of change of pore area phik
(6) Calculating a water sensitivity predicted value Y;
Figure BDA0002990325890000021
and (2) when the sample is preprocessed in the step (1), finding a point in the cubic sample, marking the corresponding position of the point and the sample seat, and when the surface of the cubic sample subjected to water-sensitive damage is coated with a film and then placed on the sample seat again at the same position in the step (5), placing the cubic sample at the corresponding position of the point and the sample seat, and determining that the positions of the samples on the sample seat are consistent.
Before the water-sensitive injury in the step (3) and after the water-sensitive injury in the step (5), when the cube sample is measured, the points are fixed to the target vision field, so that the in-situ fixed-point same-multiple observation is realized; and (5) determining that the measurement conditions before and after the water-sensitive injury are consistent.
The process of the water-sensitive injury in the step (4) is as follows: under the condition of simulating the formation pressure and temperature, a cubic sample is pumped out through a vacuumizing saturation device; and adding liquid into the vacuumizing saturation device to soak the cubic sample, and taking out and drying the sample after the sample is fully saturated with the liquid.
The liquid is in-situ layer water, the soaking time is 48 hours, and the drying is carried out at the temperature lower than 100 ℃.
The illite area expansion ratio phiIThe calculation process of (2) is as follows:
Figure BDA0002990325890000031
area expansion ratio phi of illite-montmorillonite mixed layerISThe calculation process of (2) is as follows:
Figure BDA0002990325890000032
rate of change of pore area ΦkThe calculation process of (2) is as follows:
Figure BDA0002990325890000033
the invention has the technical effects that:
the invention determines the water-sensitive mineral types to be illite and illite mixed layers by utilizing a scanning electron microscope and combining an energy spectrum method, and records the illite content WI% and content of illite-montmorillonite mixed layer WISPercent; carrying out water-sensitive damage treatment on the cubic sample by a vacuumizing saturation device; calculating the area expansion rate phi of illite minerals before and after the water-sensitive damage of the cubic sampleIIllite-smectite mixed layer mineral area expansion rate phiISAnd rate of change of pore area phik(ii) a And then calculating to obtain a water sensitivity prediction value. The water-sensitive damage degree is obtained through the method, the microscopic changes of the sensitive mineral types, mineral occurrence states, structures and pore spaces can be obtained, the method is a simple, rapid and visual shale reservoir water-sensitive prediction method, and more real bases can be provided for researching the shale reservoir water-sensitive damage mechanism. Meanwhile, the problem that the surface porosity cannot be quantitatively analyzed due to the fact that the shale sample casting body slice is difficult to die-cast is solved, pollution is reduced, the process is simplified, and meanwhile cost is saved.
Detailed Description
Example 1
A prediction method suitable for water sensitivity of a shale reservoir comprises the following steps:
(1) sample pretreatment: making the core sample into a cubic sample of 10mm multiplied by 10 mm; fixing the cubic sample on a sample holder;
(2) selecting a target view, and finding the illite and the illite-montmorillonite mixed layer in the cubic sample in the target view; and confirming that illite and illite-montmorillonite mixed layers exist in the target vision field through an energy spectrum analysis system;
(3) measuring the content W of illite in the cubic sampleI% content of illite-montmorillonite mixed layer WIS%;
Measurement of illite mineral area S of cubic sample before water-sensitive injuryI frontIllite-smectite mixed layer mineral area SIS frontAnd pore area Sk front
(4) Taking down the cubic sample from the sample seat for water-sensitive damage;
(5) coating a film on the surface of the cube sample subjected to water-sensitive injury, and then placing the cube sample on the sample seat again at the same position; measurement of illite mineral area S of cube sample after water-sensitive injuryAfter IIllite-smectite mixed layer mineral area SAfter ISAnd pore area Sk is back(ii) a Calculating the area expansion rate phi of the illite mineralIIllite-smectite mixed layer mineral area expansion rate phiISAnd rate of change of pore area phik
(6) Calculating a water sensitivity predicted value Y;
Figure BDA0002990325890000034
example 2
On the basis of embodiment 1, the method further comprises the following steps:
and (2) when the sample is preprocessed in the step (1), finding a point in the cubic sample, marking the corresponding position of the point and the sample seat, and when the surface of the cubic sample subjected to water-sensitive damage is coated with a film and then placed on the sample seat again at the same position in the step (5), placing the cubic sample at the corresponding position of the point and the sample seat, and determining that the positions of the samples on the sample seat are consistent.
Before the water-sensitive injury in the step (3) and after the water-sensitive injury in the step (5), when the cube sample is measured, the points are fixed to the target vision field, so that the in-situ fixed-point same-multiple observation is realized; and (5) determining that the measurement conditions before and after the water-sensitive injury are consistent.
The process of the water-sensitive injury in the step (4) is as follows: under the condition of simulating the formation pressure and temperature, a cubic sample is pumped out through a vacuumizing saturation device; and adding liquid into the vacuumizing saturation device to soak the cubic sample, and taking out and drying the sample after the sample is fully saturated with the liquid. Wherein the liquid is in-situ layer water, the soaking time is 48h, and the drying is carried out at the temperature lower than 100 ℃.
Example 3
On the basis of embodiment 2, the method further comprises the following steps:
the illite area expansion ratio phiIThe calculation process of (2) is as follows:
Figure BDA0002990325890000041
area expansion ratio phi of illite-montmorillonite mixed layerISThe calculation process of (2) is as follows:
Figure BDA0002990325890000042
rate of change of pore area ΦkThe calculation process of (2) is as follows:
Figure BDA0002990325890000043
examples of specific applications
Taking the prolonged south shale reservoir of the oil field as an example, the prediction method provided by the invention can obtain the following data:
Figure BDA0002990325890000044
Figure BDA0002990325890000051
finally, the water sensitivity predicted value Y and the water sensitivity index measured value I are respectively carried out on 20 groups of data by utilizing a root mean square error formulaWAnd (4) standard error calculation is carried out, wherein error values of three groups of data are large, account for 15% of the total test data and are lower than 80% of the total test data, the root mean square error (delta) of 20 groups of data is 1.52, and the delta value is smaller than 2, so that the accuracy of the predicted value Y of the water sensitivity of the mathematical model is good, and the predicted value Y can be used as an index for judging the water sensitivity of the shale reservoir.

Claims (6)

1. A prediction method suitable for water sensitivity of a shale reservoir is characterized by comprising the following steps: the method comprises the following steps:
(1) sample pretreatment: making the core sample into a cubic sample of 10mm multiplied by 10 mm; fixing the cubic sample on a sample holder;
(2) selecting a target view, and finding the illite and the illite-montmorillonite mixed layer in the cubic sample in the target view; and confirming that illite and illite-montmorillonite mixed layers exist in the target vision field through an energy spectrum analysis system;
(3) measuring the content W of illite in the cubic sampleI% content of illite-montmorillonite mixed layer WIS%;
Measurement of illite mineral area S of cubic sample before water-sensitive injuryI frontIllite-smectite mixed layer mineral area SIS frontAnd pore area Sk front
(4) Taking down the cubic sample from the sample seat for water-sensitive damage;
(5) coating a film on the surface of the cube sample subjected to water-sensitive injury, and then placing the cube sample on the sample seat again at the same position; measurement of illite mineral area S of cube sample after water-sensitive injuryAfter IIllite-smectite mixed layer mineral area SAfter ISAnd pore area Sk is back(ii) a Calculating the area expansion rate phi of the illite mineralIIllite-smectite mixed layer mineral area expansion rate phiISAnd rate of change of pore area phik
(6) Calculating a water sensitivity predicted value Y;
Figure FDA0002990325880000011
2. the method for predicting the water sensitivity of the shale reservoir according to claim 1, wherein the method comprises the following steps: and (2) when the sample is preprocessed in the step (1), finding a point in the cubic sample, marking the corresponding position of the point and the sample seat, and when the surface of the cubic sample subjected to water-sensitive damage is coated with a film and then placed on the sample seat again at the same position in the step (5), placing the cubic sample at the corresponding position of the point and the sample seat, and determining that the positions of the samples on the sample seat are consistent.
3. The method for predicting the water sensitivity of the shale reservoir according to claim 2, wherein the method comprises the following steps: before the water-sensitive injury in the step (3) and after the water-sensitive injury in the step (5), when the cube sample is measured, the points are fixed to the target vision field, so that the in-situ fixed-point same-multiple observation is realized; and (5) determining that the measurement conditions before and after the water-sensitive injury are consistent.
4. The method for predicting the water sensitivity of the shale reservoir according to claim 3, wherein the method comprises the following steps: the process of the water-sensitive injury in the step (4) is as follows: under the condition of simulating the formation pressure and temperature, a cubic sample is pumped out through a vacuumizing saturation device; and adding liquid into the vacuumizing saturation device to soak the cubic sample, and taking out and drying the sample after the sample is fully saturated with the liquid.
5. The method for predicting the water sensitivity of the shale reservoir according to claim 4, wherein the method comprises the following steps: the liquid is in-situ layer water, the soaking time is 48 hours, and the drying is carried out at the temperature lower than 100 ℃.
6. The method for predicting the water sensitivity of the shale reservoir according to claim 5, wherein the method comprises the following steps: the illite area expansion ratio phiIThe calculation process of (2) is as follows:
Figure FDA0002990325880000012
area expansion ratio phi of illite-montmorillonite mixed layerISThe calculation process of (2) is as follows:
Figure FDA0002990325880000013
rate of change of pore area ΦkThe calculation process of (2) is as follows:
Figure FDA0002990325880000021
CN202110312332.3A 2021-03-24 2021-03-24 Prediction method suitable for water sensitivity of shale reservoir Pending CN113049784A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110312332.3A CN113049784A (en) 2021-03-24 2021-03-24 Prediction method suitable for water sensitivity of shale reservoir

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110312332.3A CN113049784A (en) 2021-03-24 2021-03-24 Prediction method suitable for water sensitivity of shale reservoir

Publications (1)

Publication Number Publication Date
CN113049784A true CN113049784A (en) 2021-06-29

Family

ID=76514683

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110312332.3A Pending CN113049784A (en) 2021-03-24 2021-03-24 Prediction method suitable for water sensitivity of shale reservoir

Country Status (1)

Country Link
CN (1) CN113049784A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113984620A (en) * 2021-10-25 2022-01-28 中国科学院武汉岩土力学研究所 Uranium reservoir acidification and permeation-increasing transformation evaluation method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105067794A (en) * 2015-07-24 2015-11-18 成都理工大学 Method for testing water sensitivity, salt sensitivity and alkali sensitivity of shale reservoir stratum
CN108897906A (en) * 2018-05-24 2018-11-27 西安石油大学 A kind of reservoir sensitivity damage analogy method based on digital cores model
CN111694856A (en) * 2020-06-11 2020-09-22 中国石油大学(北京) Intelligent prediction method and device for reservoir sensitivity
CN111694855A (en) * 2020-06-11 2020-09-22 中国石油大学(北京) Intelligent prediction data processing method and device for reservoir sensitivity

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN105067794A (en) * 2015-07-24 2015-11-18 成都理工大学 Method for testing water sensitivity, salt sensitivity and alkali sensitivity of shale reservoir stratum
CN108897906A (en) * 2018-05-24 2018-11-27 西安石油大学 A kind of reservoir sensitivity damage analogy method based on digital cores model
CN111694856A (en) * 2020-06-11 2020-09-22 中国石油大学(北京) Intelligent prediction method and device for reservoir sensitivity
CN111694855A (en) * 2020-06-11 2020-09-22 中国石油大学(北京) Intelligent prediction data processing method and device for reservoir sensitivity

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴广 等: "储层水敏系统评价方法与水敏机理研究", 《石油天然气学报》 *
蒲春生 等著: "《异常应力构造低渗油藏大段泥页岩井壁稳定》", 31 December 2015, 中国石油大学出版社 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113984620A (en) * 2021-10-25 2022-01-28 中国科学院武汉岩土力学研究所 Uranium reservoir acidification and permeation-increasing transformation evaluation method
CN113984620B (en) * 2021-10-25 2022-11-22 中国科学院武汉岩土力学研究所 Uranium reservoir acidizing permeability-increasing reconstructability evaluation method

Similar Documents

Publication Publication Date Title
US10161891B1 (en) Method for characterizing rock physical characteristics of deeply buried carbonate rocks
Geffen et al. Experimental investigation of factors affecting laboratory relative permeability measurements
CN108956435B (en) Simulation experiment method and device for high-temperature high-pressure reservoir corrosion
CN108627533A (en) Fluid employs the nuclear magnetic resonance experiment method and device of feature in a kind of measurement porous media
CN108590601B (en) Experimental method for optimizing water injection expansion construction parameters
CN110595953B (en) Experimental test device and method for shale mixing wettability
US20240027379A1 (en) Method for quantitative evaluation on sensitivity of shale oil and gas reservoir to injected fluids
CN113514382B (en) Method for evaluating water film thickness after water phase imbibition flowback of gas reservoir rock containing swelling clay
MXPA06005804A (en) Method and apparatus for measuring the wettability of geological formations.
Xin et al. Pore structure evaluation in ultra-deep tight sandstones using NMR measurements and fractal analysis
McPhee et al. Routine core analysis
CN107524437B (en) Method and system for determining opening of reservoir fracture
CN109459795B (en) Method for evaluating sealing property of porous sandstone fracture zone
Ye et al. Laboratory investigation of fluid flow and permeability evolution through shale fractures
Li et al. Dynamic evolution of shale permeability under coupled temperature and effective stress conditions
CN113933148A (en) Method and device for quantitatively analyzing oil content and reservoir space of shale in different occurrence states
CN113218834A (en) Experimental device and method for quantitatively describing seepage damage of dense gas fracturing fluid and reservoir
Liu et al. Effect of wettability on oil and water distribution and production performance in a tight sandstone reservoir
Wang et al. Microscope dynamic characterization of oil charging in tight sandstone using a physical simulation experiment
Yang et al. Permeability evolution characteristics of intact and fractured shale specimens
Yang et al. Modeling water imbibition and penetration in shales: New insights into the retention of fracturing fluids
CN113049784A (en) Prediction method suitable for water sensitivity of shale reservoir
CN104948150B (en) Method and device for determining formation displacement pressure
CN110987755A (en) Fractured reservoir core flow experiment method
Li et al. Quantitative investigation of phase characteristic effects on CO2 breakthrough pressures in unsaturated Low-Permeability sandstone

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210629

RJ01 Rejection of invention patent application after publication