CN111274528A - Reservoir fracture imbibition quality prediction method and system - Google Patents

Reservoir fracture imbibition quality prediction method and system Download PDF

Info

Publication number
CN111274528A
CN111274528A CN202010135004.6A CN202010135004A CN111274528A CN 111274528 A CN111274528 A CN 111274528A CN 202010135004 A CN202010135004 A CN 202010135004A CN 111274528 A CN111274528 A CN 111274528A
Authority
CN
China
Prior art keywords
fracture
imbibition
height
length
current
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
Application number
CN202010135004.6A
Other languages
Chinese (zh)
Other versions
CN111274528B (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.)
China University of Petroleum Beijing
Original Assignee
China University of Petroleum Beijing
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 China University of Petroleum Beijing filed Critical China University of Petroleum Beijing
Priority to CN202010135004.6A priority Critical patent/CN111274528B/en
Publication of CN111274528A publication Critical patent/CN111274528A/en
Application granted granted Critical
Publication of CN111274528B publication Critical patent/CN111274528B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations
    • 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

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Chemical & Material Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Health & Medical Sciences (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Mathematical Analysis (AREA)
  • Computational Mathematics (AREA)
  • Algebra (AREA)
  • Dispersion Chemistry (AREA)
  • Mathematical Optimization (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Filling Or Discharging Of Gas Storage Vessels (AREA)

Abstract

The invention provides a method and a system for predicting the seepage quality of a reservoir fracture. The method comprises the following steps: determining the fracture imbibition height of each fracture at the current imbibition moment; when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time; otherwise, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time; and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir. The method can quickly predict the accurate seepage quality of the crack.

Description

Reservoir fracture imbibition quality prediction method and system
Technical Field
The invention relates to the technical field of oil reservoir exploration and development, in particular to a method and a system for predicting the seepage and absorption quality of a reservoir fracture.
Background
The development of unconventional oil and gas resources such as dense gas, shale gas and the like has important significance for guaranteeing national energy safety. Different from the conventional oil and gas reservoir, the compact oil and gas reservoir has the characteristics of low pore and low permeability, complex pore throat structure, strong capillary force action and obvious spontaneous imbibition phenomenon. Large-scale hydraulic fracturing is a commonly used technical means for improving the capacity of compact reservoirs, and a large amount of fracturing fluid is injected into the reservoirs to generate artificial fractures. However, in mine applications, more than 50% of the fracturing fluid remains in the underground reservoir, and spontaneous imbibition is one of the main causes of low flowback of the fracturing fluid. Research shows that a large number of natural micro-nano cracks develop in a compact reservoir, and the natural cracks and artificial cracks generated by fracturing are intersected to form a complex crack network structure, so that research on a spontaneous imbibition mechanism of fracturing fluid in the natural/artificial cracks has great significance for development of compact gas and shale gas.
At present, the research on the spontaneous imbibition mechanism of the crack has a weighing method, a volume method, CT scanning, nuclear magnetic resonance, a neutron photography technology and the like, but the methods have the following problems: (1) the porosity of the compact sandstone is generally less than 10%, which means that the mass change of the rock core in the imbibition process is very small, namely the precision requirement on a measuring tool is very high, and the conventional measuring tool cannot achieve satisfactory precision, so that the weighing method and the volume method have large errors although the dynamic change of the imbibition mass of the rock core along with time can be researched. (2) Although the methods such as CT scanning and nuclear magnetic resonance have higher precision, the cost is high, and the dynamic change of the imbibition mass along with the time cannot be researched; although the neutron photography technology can research the dynamic change in the imbibition process, the neutron photography technology has very high requirements on experimental equipment, is very high in experimental cost, is difficult to popularize and apply, and is not beneficial to the development of oil and gas resources.
In addition, sandstone pores (matrix) are often assumed in the prior art as a cluster of capillaries, and the flow in the fracture is assumed to be flow between the plates. However, the spontaneous imbibition rule of the crack and the matrix is obviously different, and errors can be brought by describing the fracture imbibition by using a capillary imbibition model.
Disclosure of Invention
The embodiment of the invention mainly aims to provide a method and a system for predicting the seepage quality of reservoir fractures, so as to rapidly predict the accurate seepage quality of the fractures, reduce the prediction cost and further effectively guide the development process of oil and gas resources.
In order to achieve the above object, an embodiment of the present invention provides a method for predicting a permeability quality of a reservoir fracture, including:
determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture;
creating a first imbibition quality prediction model and a second imbibition quality prediction model;
when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
The embodiment of the invention also provides a system for predicting the seepage quality of the reservoir fracture, which comprises the following steps:
the crack imbibition height unit is used for determining the crack imbibition height of each crack at the current imbibition time according to the obtained current imbibition time, the first physical property parameter and the crack length of each crack;
the imbibition quality prediction model creating unit is used for creating a first imbibition quality prediction model and a second imbibition quality prediction model;
the fracture imbibition quality unit is used for inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the density of the wet-phase fluid into a first imbibition quality prediction model when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, so as to obtain the imbibition quality of the fracture at the current imbibition time; when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and the prediction unit is used for predicting the imbibition quality of the rock core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil and gas resource development scheme of the reservoir.
The embodiment of the invention also provides computer equipment which comprises a memory, a processor and a computer program stored on the memory and operated on the processor, wherein the processor realizes the steps of the reservoir fracture imbibition quality prediction method when executing the computer program.
Embodiments of the present invention further provide a computer-readable storage medium, on which a computer program is stored, where the computer program, when executed by a processor, implements the steps of the method for predicting the imbibition quality of a reservoir fracture.
The reservoir fracture imbibition quality prediction method and system firstly determine the fracture imbibition height of each fracture at the current imbibition time, then determine the imbibition quality of each fracture at the current imbibition time according to the comparison result of the fracture imbibition height and the core height, and finally predict the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust the oil and gas resource development scheme of the reservoir.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments will be briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art that other drawings can be obtained based on these drawings without creative efforts.
FIG. 1 is a flow chart of a reservoir fracture imbibition mass prediction method (without accounting for gravity) in an embodiment of the invention;
FIG. 2 is a schematic illustration of spontaneous imbibition of a wet-phase fluid in a gas-saturated fractured core in an embodiment of the invention;
FIG. 3 is a schematic diagram of a theoretical model of a fractured core comprising a cluster of tortuous slab fractures in an embodiment of the invention;
FIG. 4 is a schematic illustration of spontaneous imbibition of a wet-phase fluid within a single fracture in an embodiment of the invention;
FIG. 5 is a schematic cross-sectional view of a core containing a plurality of fractures according to an embodiment of the present disclosure;
FIG. 6 is a graphical representation of the imbibition mass and imbibition saturation of a core over time for an example of the invention (without regard to gravity);
FIG. 7 is a graph showing the variation of fracture imbibition height over time for both the prior art and the present example (without regard to gravity);
FIG. 8 is a block diagram of a reservoir fracture imbibition quality prediction system (without regard to gravity) in an embodiment of the invention;
fig. 9 is a block diagram of a computer device in the embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
As will be appreciated by one skilled in the art, embodiments of the present invention may be embodied as a system, apparatus, device, method, or computer program product. Accordingly, the present disclosure may be embodied in the form of: entirely hardware, entirely software (including firmware, resident software, micro-code, etc.), or a combination of hardware and software.
In view of the low precision and high cost of the prior art, which are not beneficial to the development of oil and gas resources, the embodiment of the invention provides a reservoir fracture imbibition quality prediction method and a reservoir fracture imbibition quality prediction system, which can rapidly predict the accurate fracture imbibition quality, reduce the prediction cost and further effectively guide the development process of the oil and gas resources. The present invention will be described in detail below with reference to the accompanying drawings.
FIG. 1 is a flow chart of a reservoir fracture imbibition mass prediction method (without regard to gravity) in an embodiment of the invention. As shown in fig. 1, the reservoir fracture imbibition quality prediction method (without considering gravity) includes:
s101: and determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture.
In one embodiment, the first physical property parameter comprises surface tension, contact angle, wet phase fluid viscosity, and average tortuosity. The fracture imbibition height of each fracture (when gravity is not considered) at the current imbibition time can be determined by the following formula:
Figure BDA0002397022590000041
wherein L issThe fracture imbibition height of the fracture (when gravity is not considered) at the current imbibition time (time t) is the fracture length l, n is the ratio of the opening degree of the fracture to the length, l is the fracture length, sigma is surface tension, theta is a contact angle, mu is wet-phase fluid viscosity, tau is average tortuosity, t is imbibition time, and the unit can be s or min.
S102: a first imbibition quality prediction model and a second imbibition quality prediction model are created.
S103: and when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time.
In one embodiment, the first imbibition quality prediction model is as follows:
M(l)=nl2ρLsn(l);
wherein M (L) is a fracture having a fracture length L (L without taking gravity into account)sNot more than H), the unit is kg; ρ is the wet phase fluid density and n (l) is the number of fractures for the fracture length l.
S104: and when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time.
In one embodiment, the second imbibition quality prediction model is as follows:
Meq(l)=ρnl2τHn(l);
wherein M iseq(l) Is a fracture of length L (L without regard to gravity)sWhen is more than H), the unit of the imbibition mass at the current imbibition time is g or kg; h is the core height in cm.
S105: and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
The reservoir fracture imbibition quality prediction method shown in fig. 1 is suitable for application scenarios in which the fracture opening is smaller than 1 micron or the fracture direction is approximately horizontal, the fracture imbibition height is the fracture imbibition height without considering gravity, the imbibition quality is the imbibition quality without considering gravity, and compared with a scheme considering gravity, a relatively accurate prediction result can be obtained more quickly. As can be seen from the flow shown in fig. 1, the reservoir fracture imbibition quality prediction method (without considering gravity) according to the embodiment of the present invention determines the fracture imbibition height of each fracture at the current imbibition time, determines the imbibition quality of each fracture at the current imbibition time according to the comparison result between the fracture imbibition height and the core height, and predicts the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time to create or adjust an oil and gas resource development scheme of the reservoir, so that the accurate fracture imbibition quality can be predicted quickly, the prediction cost is reduced, and the development process of the oil and gas resource is further effectively guided.
Before executing S101, the method further includes: creating a crack number model; acquiring a second physical parameter; and inputting each crack length and the second physical parameter into the crack number model to obtain the crack number of each crack length.
In one embodiment, the second physical parameter includes: core cross section, fractal dimension, maximum fracture length, minimum fracture length, fracture porosity, and fracture length step length. The fracture number model is as follows:
Figure BDA0002397022590000051
wherein n (l) is the number of cracks of the crack length l, l is the crack length, AfIs the cross section of the core, n is the ratio of the opening degree of the crack to the length, DfIs fractal dimension,/maxTo maximum crack length, /)minIs the minimum crack length, phifFor fracture porosity,. DELTA.l is the fracture length step in μm.
Table 1 shows physical property parameters. As shown in Table 1, the physical properties of one example were as follows:
TABLE 1
Figure BDA0002397022590000052
Figure BDA0002397022590000061
Wherein, the cross section of the core can be determined according to the pre-acquired core diameter and the core height. The cross section of the core is
Figure BDA0002397022590000062
In one embodiment, the imbibition saturation of the core at the current imbibition time can be determined according to the imbibition quality of the core at the current imbibition time, and the imbibition saturation can effectively guide the development process of oil and gas resources:
Figure BDA0002397022590000063
wherein, R is the imbibition saturation of the rock core (when gravity is not considered) at the current imbibition moment, and the unit is percent; mTThe imbibition mass of the rock core (without considering gravity) at the current imbibition moment is g or kg; d is the core diameter.
FIG. 2 is a schematic illustration of spontaneous imbibition of a wet phase fluid in a gas saturated fracture-type core in an embodiment of the invention. FIG. 3 is a schematic diagram of a theoretical model of a fractured core comprising a cluster of tortuous slab fractures in an embodiment of the invention. FIG. 4 is a schematic diagram of spontaneous imbibition of a wet phase fluid within a single fracture in an embodiment of the invention, and is also an enlarged view of the box in FIG. 3.
As shown in fig. 2, the fluid flow within the fracture 2 is very similar to the flow between the plates, thus assuming a large number of natural fractures within the core 1 as a cluster of tortuous plate fractures with fractal distribution and non-intersecting with each other, as shown in fig. 3.
As shown in fig. 4, the viscosity and density of the natural gas 4 are negligible compared to the wet phase fluid 3, and the spontaneous imbibition flow in a single fracture can be expressed as:
Figure BDA0002397022590000064
wherein q is the spontaneous imbibition flow in a single fracture in mm3S; a is the crack opening degree in mum; l is the crack length in μm; l isfThe length of the tortuosity from the imbibition front edge to the bottom surface of the rock core is m; l issThe length of a straight line from the imbibition front edge to the bottom surface of the rock core (fracture imbibition height) is m; mu is a wet phase fluidViscosity in mPa · s; sigma is the surface tension of natural gas and wet phase fluid, and the unit is mN/m; (ii) a θ is the contact angle in °; rho is the wet-phase fluid density in kg/m3(ii) a g is the acceleration of gravity in m/s2. Equation (1) can be written as a differential form as follows:
Figure BDA0002397022590000071
the relationship between the tortuosity length of the imbibition front to the core floor and the fracture imbibition height can be expressed as:
Lf=τLs。 (3)
wherein tau is the average tortuosity of the crack and is dimensionless. The formula (3) can be substituted for the formula (2):
Figure BDA0002397022590000072
equation (4) can be written as follows:
Figure BDA0002397022590000073
when gravity is not considered, equation (5) can be simplified as:
Figure BDA0002397022590000074
the solution of equation (6) is:
Figure BDA0002397022590000075
the formula (7) is a crack imbibition height expression of a single crack without considering gravity, and the crack imbibition height and
Figure BDA0002397022590000076
is in direct proportion. Wherein a is nl, n is the ratio of the opening degree of the crack to the length, and is dimensionless. Equation (7) can thus be expressed as:
Figure BDA0002397022590000077
for time t, the imbibition mass m of a single fracture without considering gravity can be expressed as:
Figure BDA0002397022590000078
when the gravity is not considered, the imbibition height of a single crack can finally reach the height of the rock core, and the imbibition mass m at the momenteqCan be expressed as:
meq=ρlaτH。 (9)
wherein m and meqThe unit of (b) is g or kg.
FIG. 5 is a cross-sectional schematic view of a core containing a plurality of fractures according to an embodiment of the present disclosure. As shown in FIG. 5, the core contains a large number of natural fractures, and a representative unit A is selected from the cross section of the coreuTo study, assume representative Unit AuThe internal crack length follows fractal distribution:
Figure BDA0002397022590000081
wherein N (l ≧ zeta) is the number of cracks whose length is greater than zeta, lmaxMaximum crack length in μm; dfFractal dimension, dimensionless. Taking the derivative of both ends of equation (10) with respect to l yields the following equation:
Figure BDA0002397022590000082
thus representative Unit AuTotal pore area A ofpCan be calculated to yield:
Figure BDA0002397022590000083
further representative Unit AuCan be calculated as:
Figure BDA0002397022590000084
wherein phifIs the crack porosity in%; lminIs the minimum crack length in μm. Core Cross section AfThe total number of cracks in (a) can be calculated as:
Figure BDA0002397022590000085
wherein A isfIs the cross section of the core and has the unit of cm2,AuAnd ApAll units of (are cm)2;NfAnd the number of cracks with the length larger than l in the cross section of the core is (more than or equal to l). Derivation of equation (14) for l can be found:
Figure BDA0002397022590000086
the number of fractures n (l) with a fracture length l can be calculated according to equation (14):
Figure BDA0002397022590000091
for the whole core, the fracture imbibition mass with length l is the product of the imbibition mass of a single fracture and the number of the single fracture. When gravity is not considered and imbibition stops:
Figure BDA0002397022590000092
when gravity is not considered and the imbibition is not stopped, the imbibition mass of a fracture of length l at time t is:
Figure BDA0002397022590000093
at a certain moment, the fracture in the core consists of two parts, namely, seepage stopping and seepage not stopping, so that the seepage mass of the core without considering the gravity can be expressed as:
Figure BDA0002397022590000094
wherein lcThe critical crack length, which represents the length of the crack at which the imbibition just stopped at time t, is expressed in μm, and is expressed as follows:
Figure BDA0002397022590000095
the concrete flow of the reservoir fracture imbibition quality prediction method (without considering gravity) is as follows:
1. creating a crack number model; acquiring a second physical parameter; and inputting each crack length and the second physical parameter into the crack number model to obtain the crack number of each crack length.
2. And determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture.
3. A first imbibition quality prediction model and a second imbibition quality prediction model are created.
4. And when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time.
5. And when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time.
6. And predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
7. And determining the imbibition saturation of the core at the current imbibition time according to the imbibition quality of the core at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
FIG. 6 is a graphical representation of the imbibition mass and imbibition saturation of a core as a function of time for an example of the invention (without regard to gravity). As shown in fig. 6, the abscissa in fig. 6 is time in units of min; the ordinate is the imbibition mass and the imbibition saturation in g and% respectively.
FIG. 7 is a graph showing the time-dependent change in fracture imbibition height for both the prior art and the present examples (without considering gravity). As shown in fig. 7, the abscissa in fig. 7 is time in units of s; the ordinate is the fracture imbibition height in mm. As shown in fig. 7, the fracture imbibition height obtained by using the capillary imbibition model in the prior art is obviously different from that of the present application, and the error will gradually expand with time. Compared with the prior art, the result obtained by adopting the tortuosity flat plate crack model is more accurate.
In summary, the reservoir fracture imbibition quality prediction method (without considering gravity) of the embodiment of the invention determines the fracture imbibition height of each fracture at the current imbibition time, determines the imbibition quality of each fracture at the current imbibition time according to the comparison result of the fracture imbibition height and the core height, and predicts the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time to create or adjust an oil and gas resource development scheme of the reservoir, so that the accurate fracture imbibition quality can be rapidly predicted, the prediction cost is reduced, and the development process of the oil and gas resource is further effectively guided.
Based on the same inventive concept, the embodiment of the invention also provides a system for predicting the permeability quality of the reservoir fracture, and as the problem solving principle of the system is similar to that of a method for predicting the permeability quality of the reservoir fracture, the implementation of the system can refer to the implementation of the method, and repeated parts are not repeated.
Fig. 8 is a block diagram of the structure of a reservoir fracture imbibition quality prediction system (without considering gravity) in an embodiment of the invention. As shown in fig. 8, the reservoir fracture imbibition quality prediction system (without considering gravity) includes:
the crack imbibition height unit is used for determining the crack imbibition height of each crack at the current imbibition time according to the obtained current imbibition time, the first physical property parameter and the crack length of each crack;
the imbibition quality prediction model creating unit is used for creating a first imbibition quality prediction model and a second imbibition quality prediction model;
the fracture imbibition quality unit is used for inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the density of the wet-phase fluid into a first imbibition quality prediction model when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, so as to obtain the imbibition quality of the fracture at the current imbibition time; when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and the prediction unit is used for predicting the imbibition quality of the rock core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil and gas resource development scheme of the reservoir.
In one embodiment, the method further comprises the following steps:
a fracture number model creation unit for creating a fracture number model;
an acquisition unit configured to acquire a second physical parameter;
and the crack number unit is used for inputting each crack length and the second physical parameter into the crack number model to obtain the crack number of each crack length.
In one embodiment, the first physical property parameter comprises surface tension, contact angle, wet phase fluid viscosity, and average tortuosity;
the second physical parameter includes: core cross section, fractal dimension, maximum fracture length, minimum fracture length, fracture porosity, and fracture length step length.
In one embodiment, the method further comprises the following steps:
and the cross section determining unit is used for determining the cross section of the core according to the pre-acquired diameter and height of the core.
In summary, the reservoir fracture imbibition quality prediction method (without considering gravity) of the embodiment of the invention determines the fracture imbibition height of each fracture at the current imbibition time, determines the imbibition quality of each fracture at the current imbibition time according to the comparison result of the fracture imbibition height and the core height, and predicts the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time to create or adjust an oil and gas resource development scheme of the reservoir, so that the accurate fracture imbibition quality can be rapidly predicted, the prediction cost is reduced, and the development process of the oil and gas resource is further effectively guided.
The embodiment of the invention also provides a specific implementation mode of computer equipment capable of realizing all the steps in the reservoir fracture imbibition quality prediction method in the embodiment. Fig. 9 is a block diagram of a computer device in an embodiment of the present invention, and referring to fig. 9, the computer device specifically includes the following:
a processor (processor)901 and a memory (memory) 902.
The processor 901 is configured to call a computer program in the memory 902, and the processor executes the computer program to implement all the steps of the method for predicting the permeability quality of the reservoir fracture in the above embodiment, for example, the processor executes the computer program to implement the following steps:
determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture;
creating a first imbibition quality prediction model and a second imbibition quality prediction model;
when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
To sum up, the computer device of the embodiment of the invention firstly determines the fracture imbibition height of each fracture at the current imbibition time, then determines the imbibition quality of each fracture at the current imbibition time according to the comparison result of the fracture imbibition height and the core height, and finally predicts the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust the oil and gas resource development scheme of the reservoir, so that the accurate fracture imbibition quality can be rapidly predicted, the prediction cost is reduced, and the development process of the oil and gas resource is further effectively guided.
An embodiment of the present invention further provides a computer-readable storage medium capable of implementing all the steps in the method for predicting the imbibition quality of a reservoir fracture in the foregoing embodiment, where the computer-readable storage medium stores a computer program, and the computer program, when executed by a processor, implements all the steps in the method for predicting the imbibition quality of a reservoir fracture in the foregoing embodiment, for example, the processor implements the following steps when executing the computer program:
determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture;
creating a first imbibition quality prediction model and a second imbibition quality prediction model;
when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the wet-phase fluid density into a first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into a second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
To sum up, the computer-readable storage medium of the embodiment of the invention determines the fracture imbibition height of each fracture at the current imbibition time, determines the imbibition quality of each fracture at the current imbibition time according to the comparison result of the fracture imbibition height and the core height, and predicts the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time to create or adjust the oil and gas resource development scheme of the reservoir, so that the accurate fracture imbibition quality can be rapidly predicted, the prediction cost is reduced, and the development process of the oil and gas resource is further effectively guided.
The above-mentioned embodiments are intended to illustrate the objects, technical solutions and advantages of the present invention in further detail, and it should be understood that the above-mentioned embodiments are only exemplary embodiments of the present invention, and are not intended to limit the scope of the present invention, and any modifications, equivalent substitutions, improvements and the like made within the spirit and principle of the present invention should be included in the scope of the present invention.
Those of skill in the art will further appreciate that the various illustrative logical blocks, units, and steps described in connection with the embodiments disclosed herein may be implemented as electronic hardware, computer software, or combinations of both. To clearly illustrate the interchangeability of hardware and software, various illustrative components, elements, and steps have been described above generally in terms of their functionality. Whether such functionality is implemented as hardware or software depends upon the particular application and design requirements of the overall system. Skilled artisans may implement the described functionality in varying ways for each particular application, but such implementation decisions should not be interpreted as causing a departure from the scope of the present embodiments.
The various illustrative logical blocks, or elements, or devices described in connection with the embodiments disclosed herein may be implemented or performed with a general purpose processor, a digital signal processor, an Application Specific Integrated Circuit (ASIC), a field programmable gate array or other programmable logic device, discrete gate or transistor logic, discrete hardware components, or any combination thereof designed to perform the functions described herein. A general-purpose processor may be a microprocessor, but in the alternative, the processor may be any conventional processor, controller, microcontroller, or state machine. A processor may also be implemented as a combination of computing devices, e.g., a digital signal processor and a microprocessor, a plurality of microprocessors, one or more microprocessors in conjunction with a digital signal processor core, or any other similar configuration.
The steps of a method or algorithm described in connection with the embodiments disclosed herein may be embodied directly in hardware, in a software module executed by a processor, or in a combination of the two. A software module may be stored in RAM memory, flash memory, ROM memory, EPROM memory, EEPROM memory, registers, hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art. For example, a storage medium may be coupled to the processor such the processor can read information from, and write information to, the storage medium. In the alternative, the storage medium may be integral to the processor. The processor and the storage medium may reside in an ASIC, which may be located in a user terminal. In the alternative, the processor and the storage medium may reside in different components in a user terminal.
In one or more exemplary designs, the functions described above in connection with the embodiments of the invention may be implemented in hardware, software, firmware, or any combination of the three. If implemented in software, the functions may be stored on or transmitted over as one or more instructions or code on a computer-readable medium. Computer-readable media includes both computer storage media and communication media that facilitate transfer of a computer program from one place to another. Storage media may be any available media that can be accessed by a general purpose or special purpose computer. For example, such computer-readable media can include, but is not limited to, RAM, ROM, EEPROM, CD-ROM or other optical disk storage, magnetic disk storage or other magnetic storage devices, or any other medium which can be used to carry or store program code in the form of instructions or data structures and which can be read by a general-purpose or special-purpose computer, or a general-purpose or special-purpose processor. Additionally, any connection is properly termed a computer-readable medium, and, thus, is included if the software is transmitted from a website, server, or other remote source via a coaxial cable, fiber optic cable, twisted pair, Digital Subscriber Line (DSL), or wirelessly, e.g., infrared, radio, and microwave. Such discs (disk) and disks (disc) include compact disks, laser disks, optical disks, DVDs, floppy disks and blu-ray disks where disks usually reproduce data magnetically, while disks usually reproduce data optically with lasers. Combinations of the above may also be included in the computer-readable medium.

Claims (10)

1. A reservoir fracture imbibition quality prediction method is characterized by comprising the following steps:
determining the fracture imbibition height of each fracture at the current imbibition moment according to the obtained current imbibition moment, the first physical property parameter and the fracture length of each fracture;
creating a first imbibition quality prediction model and a second imbibition quality prediction model;
when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the density of the wet-phase fluid into the first imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into the second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and predicting the imbibition quality of the core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil-gas resource development scheme of the reservoir.
2. The method for predicting the seepage quality of the reservoir fractures according to claim 1, further comprising:
creating a crack number model;
acquiring a second physical parameter;
and inputting each crack length and the second physical parameter into the crack number model to obtain the crack number of each crack length.
3. The reservoir fracture imbibition quality prediction method of claim 2, wherein:
the first physical property parameters include surface tension, contact angle, wet phase fluid viscosity, and average tortuosity;
the second physical parameter includes: core cross section, fractal dimension, maximum fracture length, minimum fracture length, fracture porosity, and fracture length step length.
4. The method for predicting the seepage quality of the reservoir fractures according to claim 3, further comprising:
and determining the cross section of the core according to the pre-acquired diameter of the core and the height of the core.
5. A reservoir fracture imbibition quality prediction system, comprising:
the crack imbibition height unit is used for determining the crack imbibition height of each crack at the current imbibition time according to the obtained current imbibition time, the first physical property parameter and the crack length of each crack;
the imbibition quality prediction model creating unit is used for creating a first imbibition quality prediction model and a second imbibition quality prediction model;
the fracture imbibition quality unit is used for inputting the fracture length of the fracture, the fracture imbibition height of the fracture at the current imbibition time, the fracture number of the fracture length and the density of the wet-phase fluid into the first imbibition quality prediction model when the fracture imbibition height of the fracture at the current imbibition time is less than or equal to the core height, so as to obtain the imbibition quality of the fracture at the current imbibition time; when the fracture imbibition height of the fracture at the current imbibition time is greater than the core height, inputting the fracture length, the core height, the fracture number of the fracture length, the wet-phase fluid density and the average tortuosity of the fracture into the second imbibition quality prediction model to obtain the imbibition quality of the fracture at the current imbibition time;
and the prediction unit is used for predicting the imbibition quality of the rock core at the current imbibition time according to the imbibition quality of each fracture at the current imbibition time so as to create or adjust an oil and gas resource development scheme of the reservoir.
6. The reservoir fracture imbibition quality prediction system of claim 5, further comprising:
a fracture number model creation unit for creating a fracture number model;
an acquisition unit configured to acquire a second physical parameter;
and the crack number unit is used for inputting each crack length and the second physical parameter into the crack number model to obtain the crack number of each crack length.
7. The reservoir fracture imbibition quality prediction system of claim 6, wherein:
the first physical property parameters include surface tension, contact angle, wet phase fluid viscosity, and average tortuosity;
the second physical parameter includes: core cross section, fractal dimension, maximum fracture length, minimum fracture length, fracture porosity, and fracture length step length.
8. The reservoir fracture imbibition quality prediction system of claim 7, further comprising:
and the cross section determining unit is used for determining the cross section of the core according to the pre-acquired diameter of the core and the height of the core.
9. A computer device comprising a memory, a processor and a computer program stored on the memory and executed on the processor, wherein the processor when executing the computer program implements the steps of the method of reservoir fracture imbibition quality prediction of any one of claims 1 to 4.
10. A computer-readable storage medium, having stored thereon a computer program, when being executed by a processor, for carrying out the steps of the method for reservoir fracture imbibition quality prediction according to any one of claims 1 to 4.
CN202010135004.6A 2020-03-02 2020-03-02 Reservoir fracture imbibition quality prediction method and system Active CN111274528B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010135004.6A CN111274528B (en) 2020-03-02 2020-03-02 Reservoir fracture imbibition quality prediction method and system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010135004.6A CN111274528B (en) 2020-03-02 2020-03-02 Reservoir fracture imbibition quality prediction method and system

Publications (2)

Publication Number Publication Date
CN111274528A true CN111274528A (en) 2020-06-12
CN111274528B CN111274528B (en) 2021-09-17

Family

ID=70999278

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010135004.6A Active CN111274528B (en) 2020-03-02 2020-03-02 Reservoir fracture imbibition quality prediction method and system

Country Status (1)

Country Link
CN (1) CN111274528B (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111305805A (en) * 2020-03-02 2020-06-19 中国石油大学(北京) Reservoir fracture imbibition quality prediction method and system
CN112163360A (en) * 2020-10-13 2021-01-01 中国石油大学(北京) Crude oil recovery rate prediction method and system based on surfactant spontaneous imbibition
CN112329975A (en) * 2020-09-04 2021-02-05 中国石油天然气股份有限公司 Oil-gas resource prediction method and device, computer equipment and storage medium
CN113554613A (en) * 2021-07-21 2021-10-26 中国电子科技集团公司信息科学研究院 Image processing method and device based on fractal theory
CN113836767A (en) * 2021-09-18 2021-12-24 中国石油大学(华东) Optimization method for shut-in time after lamination of shale oil reservoir
CN114112771A (en) * 2020-08-25 2022-03-01 中国石油天然气股份有限公司 Rock core imbibition quality analysis method and device for oil and gas reservoir
CN114165205A (en) * 2021-12-06 2022-03-11 西安石油大学 Method for calculating cross-well flux of fracturing fluid considering imbibition

Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020016702A1 (en) * 2000-05-26 2002-02-07 Emmanuel Manceau Method for modelling flows in a fractured medium crossed by large fractures
GB2461389A (en) * 2008-07-01 2010-01-06 Logined Bv Modelling nonlinear hysteresis response of reservoir media
US20100250216A1 (en) * 2009-03-24 2010-09-30 Chevron U.S.A. Inc. System and method for characterizing fractures in a subsurface reservoir
US20120109603A1 (en) * 2009-06-22 2012-05-03 Ning Li Quantitative calculation method for oil (gas) saturation of fractured reservoir
US20120136578A1 (en) * 2010-10-20 2012-05-31 The Petroleum Institute Rise in core wettability characterization method
CN105676309A (en) * 2016-01-28 2016-06-15 中国石油大学(华东) Fractured reservoir reverse imbibition recovery prediction method
US20160341652A1 (en) * 2015-05-22 2016-11-24 Aramco Services Company Method for determining unconventional liquid imbibition in low-permeability materials
CN106370582A (en) * 2016-10-28 2017-02-01 陕西延长石油(集团)有限责任公司研究院 Fractured extra-low-permeability reservoir dynamic imbibition stimulation experiment device and application thereof
US20180106708A1 (en) * 2015-05-20 2018-04-19 Schlumberger Technology Corporation Hydraulic fracturability index using high resolution core measurements
CN108153933A (en) * 2017-11-30 2018-06-12 中国石油天然气股份有限公司 A kind of compact reservoir production forecast method and device
CN108548745A (en) * 2018-03-27 2018-09-18 中国石油大学(北京) A kind of spontaneous imbibition test method and its device of rock core
CN108827986A (en) * 2018-06-22 2018-11-16 数岩科技(厦门)股份有限公司 A kind of rock core self-priming experimental method based on CT scan
CN109682850A (en) * 2018-12-24 2019-04-26 西南石油大学 A kind of online imbibition experiment nuclear magnetic resonance test device and experimental method
CA3037513A1 (en) * 2019-03-08 2019-05-23 Research Institute Of Shaanxi Yanchang Petroleum Group, Ltd. Systems, apparatuses and methods for studying moderate waterflooding in fractured tight oil reservoirs
CN109884269A (en) * 2019-02-02 2019-06-14 中国石油大学(北京) The spontaneous imbibition Forecasting method for enhancing oil recovery of rock core and system
CN109902918A (en) * 2019-01-16 2019-06-18 西南石油大学 A kind of shale pressure break self-supporting crack seam net calculation of permeability
WO2019155184A1 (en) * 2018-02-07 2019-08-15 Independence Oilfield Chemicals Llc Fluids and methods
CN110159251A (en) * 2019-06-26 2019-08-23 西南石油大学 A method of test prediction shale fracturing effect
CN110173251A (en) * 2019-06-26 2019-08-27 中国石油集团渤海钻探工程有限公司 Compact oil reservoir CO2Auxiliary energy-storage imbibition fracturing process
CN110306960A (en) * 2019-05-22 2019-10-08 西安石油大学 The double fracturing fluid imbibition volume increase analogy methods for seeping medium reservoir of diplopore and device
CN110441206A (en) * 2019-07-26 2019-11-12 中国石油大学(北京) Integrated shale imbibition device is cut in imbibition and imbibition efficiency parameters determine method
CN110472372A (en) * 2019-09-10 2019-11-19 中国石油大学(北京) Permeability Prediction method and system based on dual media

Patent Citations (22)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20020016702A1 (en) * 2000-05-26 2002-02-07 Emmanuel Manceau Method for modelling flows in a fractured medium crossed by large fractures
GB2461389A (en) * 2008-07-01 2010-01-06 Logined Bv Modelling nonlinear hysteresis response of reservoir media
US20100250216A1 (en) * 2009-03-24 2010-09-30 Chevron U.S.A. Inc. System and method for characterizing fractures in a subsurface reservoir
US20120109603A1 (en) * 2009-06-22 2012-05-03 Ning Li Quantitative calculation method for oil (gas) saturation of fractured reservoir
US20120136578A1 (en) * 2010-10-20 2012-05-31 The Petroleum Institute Rise in core wettability characterization method
US20180106708A1 (en) * 2015-05-20 2018-04-19 Schlumberger Technology Corporation Hydraulic fracturability index using high resolution core measurements
US20160341652A1 (en) * 2015-05-22 2016-11-24 Aramco Services Company Method for determining unconventional liquid imbibition in low-permeability materials
CN105676309A (en) * 2016-01-28 2016-06-15 中国石油大学(华东) Fractured reservoir reverse imbibition recovery prediction method
CN106370582A (en) * 2016-10-28 2017-02-01 陕西延长石油(集团)有限责任公司研究院 Fractured extra-low-permeability reservoir dynamic imbibition stimulation experiment device and application thereof
CN108153933A (en) * 2017-11-30 2018-06-12 中国石油天然气股份有限公司 A kind of compact reservoir production forecast method and device
WO2019155184A1 (en) * 2018-02-07 2019-08-15 Independence Oilfield Chemicals Llc Fluids and methods
CN108548745A (en) * 2018-03-27 2018-09-18 中国石油大学(北京) A kind of spontaneous imbibition test method and its device of rock core
CN108827986A (en) * 2018-06-22 2018-11-16 数岩科技(厦门)股份有限公司 A kind of rock core self-priming experimental method based on CT scan
CN109682850A (en) * 2018-12-24 2019-04-26 西南石油大学 A kind of online imbibition experiment nuclear magnetic resonance test device and experimental method
CN109902918A (en) * 2019-01-16 2019-06-18 西南石油大学 A kind of shale pressure break self-supporting crack seam net calculation of permeability
CN109884269A (en) * 2019-02-02 2019-06-14 中国石油大学(北京) The spontaneous imbibition Forecasting method for enhancing oil recovery of rock core and system
CA3037513A1 (en) * 2019-03-08 2019-05-23 Research Institute Of Shaanxi Yanchang Petroleum Group, Ltd. Systems, apparatuses and methods for studying moderate waterflooding in fractured tight oil reservoirs
CN110306960A (en) * 2019-05-22 2019-10-08 西安石油大学 The double fracturing fluid imbibition volume increase analogy methods for seeping medium reservoir of diplopore and device
CN110159251A (en) * 2019-06-26 2019-08-23 西南石油大学 A method of test prediction shale fracturing effect
CN110173251A (en) * 2019-06-26 2019-08-27 中国石油集团渤海钻探工程有限公司 Compact oil reservoir CO2Auxiliary energy-storage imbibition fracturing process
CN110441206A (en) * 2019-07-26 2019-11-12 中国石油大学(北京) Integrated shale imbibition device is cut in imbibition and imbibition efficiency parameters determine method
CN110472372A (en) * 2019-09-10 2019-11-19 中国石油大学(北京) Permeability Prediction method and system based on dual media

Non-Patent Citations (10)

* Cited by examiner, † Cited by third party
Title
BRABAZON, J.W. 等: "Spontaneous imbibition of a wetting fluid into a fracture with opposing fractal surfaces: theory and experimental validation.", 《FRACTALS》 *
FUYONG WANG 等: "A fractal permeability model for 2D complex tortuous fractured porous media", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 *
FUYONG WANG 等: "A mathematical model for co-current spontaneous water imbibition intooil-saturated tight sandstone Upscaling from pore-scale to core-scale withfractal approach", 《JOURNAL OF PETROLEUM SCIENCE AND ENGINEERING》 *
JIANCHAO CAI 等: "Generalized Modeling of Spontaneous Imbibition Basedon Hagen Poiseuille Flow in Tortuous Capillaries with Variably Shaped Apertures", 《LANGMUIR》 *
WANG JING 等: "Investigations on spontaneous imbibition and the influencing factors in tight oil reservoirs", 《FUEL》 *
WANG, F 等: "Fractal characterization of tight oil reservoir pore structure using nuclear magnetic resonance and mercury intrusion porosimetry", 《FRACTALS》 *
孟庆帮: "裂缝性油藏自发渗吸特征及影响因素研究", 《中国博士学位论文全文数据库 工程科技Ⅰ辑》 *
李帅 等: "考虑渗吸效应的致密油藏体积改造可行性分析", 《地质科技情报》 *
王世芳 等: "多孔介质渗透率的一种新分形模型", 《力学季刊》 *
蔡建超 等: "裂缝-孔隙型双重介质油藏渗吸机理的分形分析", 《物理学报》 *

Cited By (12)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111305805A (en) * 2020-03-02 2020-06-19 中国石油大学(北京) Reservoir fracture imbibition quality prediction method and system
CN111305805B (en) * 2020-03-02 2021-02-12 中国石油大学(北京) Reservoir fracture imbibition quality prediction method and system
CN114112771A (en) * 2020-08-25 2022-03-01 中国石油天然气股份有限公司 Rock core imbibition quality analysis method and device for oil and gas reservoir
CN114112771B (en) * 2020-08-25 2023-09-26 中国石油天然气股份有限公司 Core imbibition quality analysis method and device for oil and gas reservoir
CN112329975A (en) * 2020-09-04 2021-02-05 中国石油天然气股份有限公司 Oil-gas resource prediction method and device, computer equipment and storage medium
CN112329975B (en) * 2020-09-04 2024-04-30 中国石油天然气股份有限公司 Oil-gas resource prediction method, device, computer equipment and storage medium
CN112163360A (en) * 2020-10-13 2021-01-01 中国石油大学(北京) Crude oil recovery rate prediction method and system based on surfactant spontaneous imbibition
CN112163360B (en) * 2020-10-13 2022-08-23 中国石油大学(北京) Crude oil recovery rate prediction method and system based on surfactant spontaneous imbibition
CN113554613A (en) * 2021-07-21 2021-10-26 中国电子科技集团公司信息科学研究院 Image processing method and device based on fractal theory
CN113554613B (en) * 2021-07-21 2024-03-01 中国电子科技集团公司信息科学研究院 Image processing method and device based on fractal theory
CN113836767A (en) * 2021-09-18 2021-12-24 中国石油大学(华东) Optimization method for shut-in time after lamination of shale oil reservoir
CN114165205A (en) * 2021-12-06 2022-03-11 西安石油大学 Method for calculating cross-well flux of fracturing fluid considering imbibition

Also Published As

Publication number Publication date
CN111274528B (en) 2021-09-17

Similar Documents

Publication Publication Date Title
CN111274528B (en) Reservoir fracture imbibition quality prediction method and system
Daccord et al. Carbonate acidizing: toward a quantitative model of the wormholing phenomenon
CN104533370B (en) Pressure break horizontal well oil reservoir, crack, pit shaft coupled model method
US10731451B2 (en) Exploitation method and exploitation well pattern for coalbed methane in low-permeability coalbed
CN111305805B (en) Reservoir fracture imbibition quality prediction method and system
CN110306960B (en) Fracturing fluid imbibition production-increasing simulation method and device for double-hole double-permeability medium reservoir
Qian et al. A fluid-solid-chemistry coupling model for shale wellbore stability
Pagels et al. Measuring capillary pressure tells more than pretty pictures
CN105260543A (en) Double-pore model-based multi-medium oil gas flow simulation method and device
CN105738252A (en) Measurement method of flowable opening degree limit of thickened oil in cracks
Chuanliang et al. A study on mechanism of oil displacement by imbibition
Li et al. Tracer flowback based fracture network characterization in hydraulic fracturing
Shagiakhmetov et al. Investigation of selective properties of the gel-forming composition for the limitation of water inflow to carbonate reservoirs conditions
CN107451671B (en) Method and system for predicting initial capacity of shale formation after fracturing
CN111460651B (en) Method and system for predicting yield of crude oil absorbed by cracks
CN111335853B (en) Method and device for determining effective water injection amount of strong edge water reservoir during edge external water injection development
CN112163360B (en) Crude oil recovery rate prediction method and system based on surfactant spontaneous imbibition
Lakhtychkin et al. Modelling of transport of two proppant‐laden immiscible power‐law fluids through an expanding fracture
CN113586042B (en) Method and device for determining dynamic permeability of pore type carbonate rock
CN110080745A (en) Separate stratum fracfturing straight well PRODUCTION FORECASTING METHODS and device
CN114575812A (en) Method and device for determining fracturing scheme of shale reservoir gas well
CN111751879A (en) Reservoir porosity prediction method and device
CN114086933B (en) Well testing interpretation method and device, storage medium and computer equipment
CN112377186B (en) Dumbbell-shaped fracture-cavity type carbonate reservoir yield analysis model construction method and device
CN114060017B (en) Low-permeability oil reservoir vertical well volume transformation method and transformation system

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