CN117970482B - Method, system, medium and terminal for judging microseism event dry-wet type - Google Patents
Method, system, medium and terminal for judging microseism event dry-wet type Download PDFInfo
- Publication number
- CN117970482B CN117970482B CN202410364680.9A CN202410364680A CN117970482B CN 117970482 B CN117970482 B CN 117970482B CN 202410364680 A CN202410364680 A CN 202410364680A CN 117970482 B CN117970482 B CN 117970482B
- Authority
- CN
- China
- Prior art keywords
- fluid
- dry
- wet
- type
- model
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 47
- 239000012530 fluid Substances 0.000 claims abstract description 144
- 239000011435 rock Substances 0.000 claims abstract description 71
- 238000002347 injection Methods 0.000 claims abstract description 52
- 239000007924 injection Substances 0.000 claims abstract description 52
- 238000004088 simulation Methods 0.000 claims abstract description 30
- 230000007246 mechanism Effects 0.000 claims abstract description 24
- 238000003384 imaging method Methods 0.000 claims abstract description 15
- 238000006073 displacement reaction Methods 0.000 claims abstract description 11
- 230000008878 coupling Effects 0.000 claims abstract description 8
- 238000010168 coupling process Methods 0.000 claims abstract description 8
- 238000005859 coupling reaction Methods 0.000 claims abstract description 8
- 238000002474 experimental method Methods 0.000 claims abstract description 8
- 238000004364 calculation method Methods 0.000 claims abstract description 7
- 239000007787 solid Substances 0.000 claims abstract description 7
- 239000002245 particle Substances 0.000 claims description 33
- 206010017076 Fracture Diseases 0.000 claims description 17
- 238000010586 diagram Methods 0.000 claims description 17
- 208000010392 Bone Fractures Diseases 0.000 claims description 16
- 230000008569 process Effects 0.000 claims description 15
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 9
- 238000011160 research Methods 0.000 claims description 7
- 238000010276 construction Methods 0.000 claims description 6
- 239000011159 matrix material Substances 0.000 claims description 6
- 238000012544 monitoring process Methods 0.000 claims description 5
- 230000009471 action Effects 0.000 claims description 4
- 238000005336 cracking Methods 0.000 claims description 4
- 208000002565 Open Fractures Diseases 0.000 claims description 3
- 239000011148 porous material Substances 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 239000000243 solution Substances 0.000 description 6
- 230000006399 behavior Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000006185 dispersion Substances 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 238000011065 in-situ storage Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000750 progressive effect Effects 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/282—Application of seismic models, synthetic seismograms
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/288—Event detection in seismic signals, e.g. microseismics
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/30—Assessment of water resources
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Acoustics & Sound (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Business, Economics & Management (AREA)
- Emergency Management (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a method, a system, a medium and a terminal for judging the dry and wet types of microseism events, which comprise the following steps: determining reservoir type and reservoir parameters, and fluid properties and parameters of the fluid injection; collecting a rock core of reservoir rock, and determining rock mechanical parameters by utilizing a triaxial rock mechanical experiment; establishing a stratum model and setting rock mass and various parameters; carrying out fluid injection on the stratum model, carrying out numerical simulation circulation under the drive of time step, and obtaining different types of cracks and changes of model stress fields, speed fields, displacement fields and fluid pressure fields through fluid-solid coupling calculation and imaging; analyzing and distinguishing the dry and wet crack type according to different performances of the stress field and the fluid pressure of the crack type and the model; performing dry and wet judgment on microseism events in engineering based on a seismic source mechanism; according to the method, rock physics and numerical simulation are combined, and the generation difference of the dry and wet cracks is obtained from the essence of a rock breaking mechanism, so that the type of the microseism dry and wet event in the injection engineering is judged.
Description
Technical Field
The invention relates to the technical field of microseism monitoring, in particular to a method, a system, a medium and a terminal for judging the dry and wet types of microseism events.
Background
The microseism event is generated in the process of injecting fluid into an underground reservoir, and the rock body vibrates and propagates to the periphery in the form of earthquake waves in the process of disturbance or even fracture of the underground rock body due to local stress instability, and finally is collected and identified by a ground or well intermediate detector.
According to the distribution condition of microseism events and the time-space sequence of occurrence events, the microseism events are divided into dry events and wet events, and the fundamental problem is that whether fluid directly contacts to participate in rock mass damage or not so that the events are generated; the dry-wet event type judgment is performed to better describe and describe the sweep range of injected reservoir fluids, and has direct help to engineering management and seismic risk assessment.
At present, the dry and wet events are judged mainly according to the distance between the event occurrence position (event positioning) and the fluid injection area; however, this method of judgment is too subjective and is similar to manually judging event types;
Or judging whether the fluid exists according to the difference of microseism forward wave forms in the injection engineering area; however, the judging method starts from the vibration signal, adopts a waveform forward modeling method, does not start from the reservoir rock essence, and omits the direct factor that the mutual coupling relation between the fluid and the rock is broken and cracked in the injection or fracturing process, so that the microseism signal is received.
Therefore, how to determine the type of microseismic event dry and wet is a matter of great concern to those skilled in the art from the nature of the rock fracture mechanism.
Disclosure of Invention
In view of the above, the present invention provides a method, a system, a medium and a terminal for determining the type of dry and wet microseism event, so as to solve some of the technical problems mentioned in the background art.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a method for judging the dry and wet type of a microseism event comprises the following steps:
S1, determining the type and the parameters of a reservoir injected with fluid, and determining the properties and the parameters of the fluid injected;
S2, collecting a rock core of the reservoir rock, and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
S3, building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix by particles, and represents the existing cracks in a discrete crack network mode;
S4, setting a fluid injection position for the established stratum model, and injecting the fluid according to a fluid injection mode of engineering implementation standards;
s5, under the drive of a time step, performing numerical simulation circulation operation on the stratum model in the fluid injection process, obtaining the changes of a stress field, a speed field, a displacement field and a fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different performances of a crack pattern, the stress field of the model and the fluid pressure;
s6, analyzing and distinguishing the type of the wet and dry cracks according to the crack patterns in the imaging diagram and different performances of the stress field and the fluid pressure of the model.
Preferably, the reservoir parameters in step S1 include formation thickness, formation density, and porosity, and the fluid parameters include fluid density, fluid viscosity.
Preferably, the rock mechanical parameters in step S2 include modulus of elasticity, cohesion and slip angle.
Preferably, the fluid injection mode is constant pressure injection or constant rate injection.
Preferably, the flow of the numerical simulation cyclic operation is specifically:
s51, detecting the stratum model before the action of force and fluid in the current time step, calculating whether each particle and the whole model exceed a preset stress value F max, and detecting whether the model is damaged;
S52, calculating the interrelation of the inter-particle forces of the stratum model, calculating to obtain a stress result of an initial point based on a constitutive equation according to the fluid injection quantity and the pore channels among the particles, transmitting the stress result to the periphery through the grid from the initial point, and further calculating the stress magnitude born by other particles to obtain a normal force and a tangential force updated along with the time step;
s53, judging the force and fracture criteria, wherein when the stress is larger than the inherent bonding strength, namely the cohesive force, a channel is generated, and two crack types are obtained according to the size relation of the stress direction among particles, wherein the normal force is larger than the tangential force and is the tensile crack, and the tangential force is larger than the normal force and is the shear crack;
s54, updating a stress field, a speed field, a displacement field and a fluid pressure field to the global model to obtain a stress field of the whole model;
s55, running to the next time step, repeating the S51-S55 cycle, or reaching the preset time step, or completely destroying the model overall.
Preferably, the specific content of step S52 is:
the magnitude of the stress to which the particles are subjected is determined by the flow rate Q of the injected fluid:
Wherein, To expand the channel width for a fluid, when two particles are in contact,/>Equal to 0,/>Open channel length for fluid,/>Expressed as the average radius of two particles,/>To inject fluid viscosity,/>Is the pressure difference;
Wherein, Total flow of local channel fluid per time step,/>Open fracture surface area per time step,/>Is the total volume of cracks in the particle domain;
Viscous fluids induce shear forces at grain boundaries The method comprises the following steps:
the expression modes of the normal force and the tangential force updated along with the time step are as follows:
。
preferably, in step S6, the specific content of the analysis and identification of the type of the dry and wet cracks is as follows:
the fracture type is a tensile fracture and has fluid pressure display, and is a wet fracture;
the type of the crack is a shear crack, no fluid pressure is displayed, and the crack is a dry crack;
in step S7, the specific content of the type of the dry-wet event is determined as follows:
Waveform data in microseism monitoring is collected and processed to obtain focus positioning information, an inversion focus mechanism characterizes the type of cracks generated by cracking of a focus position through a focus mechanism ball, wet cracks are wet events, and dry cracks are dry events.
A system for judging the dry and wet type of a microseism event is based on a method for judging the dry and wet type of the microseism event, and comprises the following steps: the system comprises a reservoir and fluid parameter acquisition module, a rock mechanical parameter acquisition device, a stratum model construction module, a fluid injection module, a numerical simulation module, an analysis and discrimination module and a seismic source mechanism inversion module;
The reservoir and fluid parameter acquisition module is used for determining the reservoir type and reservoir parameters of fluid injection and determining the injected fluid properties and fluid parameters;
The rock mechanical parameter acquisition device is used for acquiring a rock core of the reservoir rock and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
The stratum model construction module is used for building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix in the form of particles, and represents the existing cracks in the form of a discrete crack network;
The fluid injection module is used for setting a fluid injection position for the established stratum model and injecting the fluid according to a fluid injection mode of engineering implementation standards;
The numerical simulation module is used for carrying out numerical simulation circulation operation on the stratum model in the fluid injection process under the drive of the time step, obtaining different types of cracks and the changes of the stress field, the speed field, the displacement field and the fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different types of cracks and different manifestations of the stress field and the fluid pressure of the model;
The analysis and discrimination module is used for analyzing and discriminating the type of the dry and wet cracks according to the type of the cracks in the imaging diagram and the different performances of the stress field and the fluid pressure of the model;
the seismic source mechanism inversion module is used for acquiring the seismic source information of the microseism event based on the seismic source mechanism and judging the type of the dry-wet event based on the type of the dry-wet crack.
A computer readable storage medium having stored thereon a computer program which when executed by a processor implements the method of determining a microseismic event dry and wet type.
The processing terminal comprises a memory and a processor, wherein a computer program capable of running on the processor is stored in the memory, and the processor realizes the method for judging the dry and wet type of the microseism event when executing the computer program.
Compared with the prior art, the method, the system, the medium and the terminal for judging the dry and wet type of the microseism event are disclosed, the dry and wet type of the event is not judged by means of the event dispersion and the distance alone, the method combining rock physics and numerical simulation is adopted, from the essence of a rock breaking mechanism, the difference between the participation of fluid and the participation of no fluid is obtained by injecting the fluid into the rock to simulate the engineering injection result, and further the generation difference of the dry and wet cracks is obtained, and the continuous behavior after the fluid is injected can be reflected better, so that the problem of judging the dry and wet event type in microseism monitoring data of the injection engineering is solved; the invention can be adjusted according to different projects to meet the implementation specificity of the corresponding fluid injection project.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings that are required to be used in the embodiments or the description of the prior art will be briefly described below, and it is obvious that the drawings in the following description are only embodiments of the present invention, and that other drawings can be obtained according to the provided drawings without inventive effort for a person skilled in the art.
FIG. 1 is a schematic diagram of a method for determining the type of the dry and wet of a microseism event according to the present invention;
FIG. 2 is a schematic diagram of formation model creation provided by the present invention;
FIG. 3 is a schematic diagram of model detection after various model parameters are input;
FIG. 4 is a schematic diagram of the primary dynamic stress field and primary dynamic velocity field prior to fluid injection in accordance with the present invention;
FIG. 5 is a schematic diagram of a numerical simulation flow provided by the present invention;
FIG. 6 is a diagram of a numerical simulation process imaging provided by the present invention;
FIG. 7 is a diagram showing the numerical simulation results provided by the present invention;
FIG. 8 is a diagram showing the results of classifying wet and dry events according to the present invention.
Detailed Description
The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The embodiment of the invention discloses a method for judging the dry and wet type of a microseism event, as shown in fig. 1, comprising the following steps:
s1, collecting basic data of a research target area: determining the type and parameters of the reservoir into which the fluid is injected, and determining the nature and parameters of the fluid to be injected;
S2, collecting a rock core of the reservoir rock, and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
In the embodiment, the core of reservoir rock is collected by drilling coring for rock experiments, different experimental instruments are adapted, and samples with different sizes are prepared;
S3, building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix by particles, and represents pre-existing cracks in the form of a discrete crack network, as shown in figures 2 and 3;
S4, setting a fluid injection position for the established stratum model, and injecting the fluid according to a fluid injection mode of engineering implementation standards;
s5, under the drive of a time step, performing numerical simulation circulation operation on the stratum model in the fluid injection process, obtaining the changes of a stress field, a speed field, a displacement field and a fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different performances of a crack pattern, the stress field of the model and the fluid pressure;
s6, analyzing and distinguishing the type of the wet and dry cracks according to the crack patterns in the imaging diagram and different performances of the stress field and the fluid pressure of the model.
To further implement the above technical solution, the reservoir parameters in step S1 include formation thickness, formation density and porosity, and the fluid parameters include fluid density and fluid viscosity.
In order to further implement the above technical solution, the rock mechanical parameters in step S2 include elastic modulus, cohesion and sliding angle.
Preferably, the fluid injection mode is constant pressure injection or constant rate injection.
In practical application, the fluid injection position is usually a round hole or a cylinder, so that the integral change of the model is convenient to observe; FIG. 2 (a) is a diagram of an initially created formation model, and FIG. 2 (b) is a diagram of a formation model center borehole;
Fig. 3 (a) and fig. 3 (b) are diagrams showing the stress field changes before and after the input of the parameters in steps S1 and S2, respectively.
In the actual engineering and numerical simulation process, with continuous fluid injection, stress field changes of in-situ reservoir rock bodies are affected, local speed and displacement are changed, corresponding rock body changes also act on injected fluid to interact until a dynamic balance state is achieved, and in the process, the behavior of the fluid is observed in real time by numerical simulation, which cannot be achieved by the actual engineering.
In order to further implement the above technical solution, as shown in fig. 4, the flow of the numerical simulation cycle operation specifically includes:
S51, detecting the stratum model before the action of force and fluid in the current time step, calculating whether each particle and the whole model exceed a preset stress value F max, and detecting whether the model is damaged, as shown in FIG. 5;
S52, calculating the interrelation of the inter-particle forces of the stratum model, calculating to obtain a stress result of an initial point based on a constitutive equation according to the fluid injection quantity and the pore channels among the particles, transmitting the stress result to the periphery through the grid from the initial point, and further calculating the stress magnitude born by other particles to obtain a normal force and a tangential force updated along with the time step;
s53, judging the force and fracture criteria, wherein when the stress is larger than the inherent bonding strength, namely the cohesive force, a channel is generated, and two crack types are obtained according to the size relation of the stress direction among particles, wherein the normal force is larger than the tangential force and is the tensile crack, and the tangential force is larger than the normal force and is the shear crack;
s54, updating a stress field, a speed field, a displacement field and a fluid pressure field to the global model to obtain a stress field of the whole model;
s55, running to the next time step, repeating the S51-S55 cycle, or reaching the preset time step, or completely destroying the model overall.
In this example, imaging the log simulation process, as in fig. 6, can yield a fracture pattern (red line in fig. 6) as well as the stress field (background change in fig. 6) and fluid pressure (colored dots in fig. 6) of the model, these different manifestations being used for the type analysis of wet and dry fractures.
In order to further implement the above technical solution, the specific content of step S52 is:
the magnitude of the stress to which the particles are subjected is determined by the flow rate Q of the injected fluid:
Wherein, To expand the channel width for a fluid, when two particles are in contact,/>Equal to 0,/>Open channel length for fluid,/>Expressed as the average radius of two particles,/>To inject fluid viscosity,/>Is the pressure difference;
Wherein, Total flow of local channel fluid per time step,/>Open fracture surface area per time step,/>Is the total volume of cracks in the particle domain;
Viscous fluids induce shear forces at grain boundaries The method comprises the following steps:
the expression modes of the normal force and the tangential force updated along with the time step are as follows:
。
In order to further implement the above technical solution, in step S6, the specific contents of analyzing and distinguishing the type of the wet and dry crack are as follows:
the fracture type is a tensile fracture and has fluid pressure display, and is a wet fracture;
the type of the crack is a shear crack, no fluid pressure is displayed, and the crack is a dry crack;
in step S7, the specific content of the type of the dry-wet event is determined as follows:
Waveform data in microseism monitoring is collected and processed to obtain focus positioning information, an inversion focus mechanism characterizes the type of cracks generated by cracking of a focus position through a focus mechanism ball, wet cracks are wet events, and dry cracks are dry events.
The specific principle is as follows: in the process of cracking, the interrelation between the fluid and the crack can be observed, and then the dry and wet cracks can be distinguished, the fundamental difference is the influence range of the fluid, whether the fluid pressure exists in the crack expansion or not is observed from numerical simulation, in addition, the crack expansion types of the fluid are different from the seismic source mechanism of the microseism event in the numerical simulation, as shown in fig. 7, the numerical simulation result shows the relation between different crack types and the fluid, the black is the wet crack, and the red is the dry crack.
In practical application, in practical engineering, according to the numerical simulation result, dry and wet events are correspondingly classified; firstly, according to the distance between the event occurrence position and the fluid, the other is the action result of different cracks in the numerical simulation process according to the event source mechanism, the dry cracks are mainly shear sliding or dislocation because no fluid participates, and the wet cracks have stretching components because of the fluid parameters and the source mechanism, as shown in fig. 8, the dry and wet event classification result (blue is wet event and black is dry event) of the combined numerical simulation and event mechanism.
A microseism event dry and wet type judging system is based on a microseism event dry and wet type judging method, comprising the following steps: the system comprises a reservoir and fluid parameter acquisition module, a rock mechanical parameter acquisition device, a stratum model construction module, a fluid injection module, a numerical simulation module, an analysis and discrimination module and a seismic source mechanism inversion module;
The reservoir and fluid parameter acquisition module is used for determining the reservoir type and reservoir parameters of fluid injection and determining the injected fluid properties and fluid parameters;
The rock mechanical parameter acquisition device is used for acquiring a rock core of the reservoir rock and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
The stratum model construction module is used for building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix in the form of particles, and represents the existing cracks in the form of a discrete crack network;
The fluid injection module is used for setting a fluid injection position for the established stratum model and injecting the fluid according to a fluid injection mode of engineering implementation standards;
The numerical simulation module is used for carrying out numerical simulation circulation operation on the stratum model in the fluid injection process under the drive of the time step, obtaining different types of cracks and the changes of the stress field, the speed field, the displacement field and the fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different types of cracks and different manifestations of the stress field and the fluid pressure of the model;
The analysis and discrimination module is used for analyzing and discriminating the type of the dry and wet cracks according to the type of the cracks in the imaging diagram and the different performances of the stress field and the fluid pressure of the model;
the seismic source mechanism inversion module is used for acquiring the seismic source information of the microseism event based on the seismic source mechanism and judging the type of the dry-wet event based on the type of the dry-wet crack.
A computer readable storage medium having stored thereon a computer program which when executed by a processor implements a method of determining a type of microseismic event dry and wet.
A processing terminal comprises a memory and a processor, wherein a computer program capable of running on the processor is stored in the memory, and the processor realizes a method for judging the dry and wet type of a microseism event when executing the computer program.
In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.
The previous description of the disclosed embodiments is provided to enable any person skilled in the art to make or use the present invention. Various modifications to these embodiments will be readily apparent to those skilled in the art, and the generic principles defined herein may be applied to other embodiments without departing from the spirit or scope of the invention. Thus, the present invention is not intended to be limited to the embodiments shown herein but is to be accorded the widest scope consistent with the principles and novel features disclosed herein.
Claims (8)
1. The method for judging the dry and wet type of the microseism event is characterized by comprising the following steps of:
S1, determining the type and the parameters of a reservoir injected with fluid, and determining the properties and the parameters of the fluid injected;
S2, collecting a rock core of the reservoir rock, and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
S3, building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix by particles, and represents the existing cracks in a discrete crack network mode;
S4, setting a fluid injection position for the established stratum model, and injecting the fluid according to a fluid injection mode of engineering implementation standards;
S5, under the drive of a time step, performing numerical simulation circulation operation on the stratum model in the fluid injection process, obtaining different types of cracks, and the changes of a stress field, a speed field, a displacement field and a fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different types of crack patterns and different performances of the stress field and the fluid pressure of the model;
s6, analyzing and distinguishing the type of the wet and dry cracks according to the type of the cracks in the imaging diagram and different performances of stress fields and fluid pressure of the model;
S7, acquiring the source information of the microseism event based on a source mechanism and judging the type of the dry and wet event based on the type of the dry and wet crack;
The numerical simulation cycle operation flow is specifically as follows:
s51, detecting the stratum model before the action of force and fluid in the current time step, calculating whether each particle and the whole model exceed a preset stress value F max, and detecting whether the model is damaged;
S52, calculating the interrelation of the inter-particle forces of the stratum model, calculating to obtain a stress result of an initial point based on a constitutive equation according to the fluid injection quantity and the pore channels among the particles, transmitting the stress result to the periphery through the grid from the initial point, and further calculating the stress magnitude born by other particles to obtain a normal force and a tangential force updated along with the time step;
s53, judging the force and fracture criteria, wherein when the stress is larger than the inherent bonding strength, namely the cohesive force, a channel is generated, and two crack types are obtained according to the size relation of the stress direction among particles, wherein the normal force is larger than the tangential force and is the tensile crack, and the tangential force is larger than the normal force and is the shear crack;
s54, updating a stress field, a speed field, a displacement field and a fluid pressure field to the global model to obtain a stress field of the whole model;
S55, running to the next time step, repeating the S51-S55 cycle, or reaching the preset time step, or completely destroying the overall model;
the specific content of step S52 is:
the magnitude of the stress to which the particles are subjected is determined by the flow rate Q of the injected fluid:
Wherein w is the fluid expansion channel width, when two particles are in contact, w is equal to 0, L p is the fluid opening channel length, L p is represented by the average radius of the two particles, μ is the injection fluid viscosity, ΔP is the pressure difference;
ΔP＝(∑Q-dVf)Vf
Wherein Σq is the total flow of local channel fluid in unit time step, dV f is the open fracture surface area in unit time step, V f is the total volume of the fracture in the particle domain;
The viscous fluid causes shear forces F s at the grain boundaries:
the expression modes of the normal force and the tangential force updated along with the time step are as follows:
Fn(t+Δt)＝Fn(t)n+ΔFn
Fs(t+Δt)＝Fs(t)n+ΔFs。
2. The method of claim 1, wherein the reservoir parameters in step S1 include formation thickness, formation density, and porosity, and the fluid parameters include fluid density and fluid viscosity.
3. The method according to claim 1, wherein the rock mechanical parameters in step S2 include elastic modulus, cohesion and slip angle.
4. The method of claim 1, wherein the fluid injection is constant pressure injection or constant rate injection.
5. The method for judging the type of the dry and wet microseism event according to claim 1, wherein the specific content of the analysis and identification of the type of the dry and wet cracks in the step S6 is as follows:
the fracture type is a tensile fracture and has fluid pressure display, and is a wet fracture;
the type of the crack is a shear crack, no fluid pressure is displayed, and the crack is a dry crack;
in step S7, the specific content of the type of the dry-wet event is determined as follows:
Waveform data in microseism monitoring is collected and processed to obtain focus positioning information, an inversion focus mechanism characterizes the type of cracks generated by cracking of a focus position through a focus mechanism ball, wet cracks are wet events, and dry cracks are dry events.
6. A system for determining a microseismic event dry-wet type, wherein the method for determining a microseismic event dry-wet type according to any one of claims 1-5 comprises: the system comprises a reservoir and fluid parameter acquisition module, a rock mechanical parameter acquisition device, a stratum model construction module, a fluid injection module, a numerical simulation module, an analysis and discrimination module and a seismic source mechanism inversion module;
The reservoir and fluid parameter acquisition module is used for determining the reservoir type and reservoir parameters of fluid injection and determining the injected fluid properties and fluid parameters;
The rock mechanical parameter acquisition device is used for acquiring a rock core of the reservoir rock and determining rock mechanical parameters of the reservoir rock by utilizing a triaxial rock mechanical experiment;
The stratum model construction module is used for building a stratum model according to actual research stratum, adding and setting rock mass and various parameters in the steps S1 and S2, wherein the stratum model represents a reservoir rock matrix in the form of particles, and represents the existing cracks in the form of a discrete crack network;
The fluid injection module is used for setting a fluid injection position for the established stratum model and injecting the fluid according to a fluid injection mode of engineering implementation standards;
The numerical simulation module is used for carrying out numerical simulation circulation operation on the stratum model in the fluid injection process under the drive of the time step, obtaining different types of cracks and the changes of the stress field, the speed field, the displacement field and the fluid pressure field of the stratum model through fluid-solid coupling calculation, and imaging different types of cracks and different manifestations of the stress field and the fluid pressure of the model;
The analysis and discrimination module is used for analyzing and discriminating the type of the dry and wet cracks according to the type of the cracks in the imaging diagram and the different performances of the stress field and the fluid pressure of the model;
the seismic source mechanism inversion module is used for acquiring the seismic source information of the microseism event based on the seismic source mechanism and judging the type of the dry-wet event based on the type of the dry-wet crack.
7. A computer readable storage medium having stored thereon a computer program, which when executed by a processor implements a method of determining the type of microseismic event dry or wet as claimed in any one of claims 1 to 5.
8. A processing terminal comprising a memory and a processor, wherein the memory stores a computer program executable on the processor, wherein the processor implements a method for determining the type of microseismic event dry and wet as claimed in any one of claims 1 to 5 when the computer program is executed by the processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410364680.9A CN117970482B (en) | 2024-03-28 | 2024-03-28 | Method, system, medium and terminal for judging microseism event dry-wet type |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410364680.9A CN117970482B (en) | 2024-03-28 | 2024-03-28 | Method, system, medium and terminal for judging microseism event dry-wet type |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117970482A CN117970482A (en) | 2024-05-03 |
| CN117970482B true CN117970482B (en) | 2024-06-21 |
Family
ID=90862847
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410364680.9A Active CN117970482B (en) | 2024-03-28 | 2024-03-28 | Method, system, medium and terminal for judging microseism event dry-wet type |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117970482B (en) |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115616659A (en) * | 2022-10-10 | 2023-01-17 | 中国矿业大学(北京) | Method and device for determining type of microseism event and electronic equipment |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20020013687A1 (en) * | 2000-03-27 | 2002-01-31 | Ortoleva Peter J. | Methods and systems for simulation-enhanced fracture detections in sedimentary basins |
| CN105089657B (en) * | 2015-06-15 | 2018-05-04 | 中国石油天然气股份有限公司 | Physical simulation method and experimental device for oil-gas filling of fracture-cavity carbonate reservoir |
| CN205787165U (en) * | 2016-07-07 | 2016-12-07 | 西南石油大学 | A kind of seismic prospecting analogue experiment installation |
| CN109061734A (en) * | 2018-05-24 | 2018-12-21 | 中国石油天然气集团有限公司 | A kind of processing method and processing device of seismic data |
| US11351682B2 (en) * | 2019-06-19 | 2022-06-07 | International Business Machines Corporation | Environment monitoring and associated monitoring device |
| CN112067482B (en) * | 2020-09-02 | 2023-11-17 | 内蒙古大学 | Intelligent simulated earthquake loading test box |
| CN115853488A (en) * | 2022-12-06 | 2023-03-28 | 吉林大学 | Multistage fracturing method for reducing cracking pressure of dry hot rock reservoir by using supercritical water |
-
2024
- 2024-03-28 CN CN202410364680.9A patent/CN117970482B/en active Active
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115616659A (en) * | 2022-10-10 | 2023-01-17 | 中国矿业大学(北京) | Method and device for determining type of microseism event and electronic equipment |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117970482A (en) | 2024-05-03 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN105319603B (en) | The Forecasting Methodology of tight sandstone reservoir complexity chicken-wire cracking | |
| Meng et al. | Experimental and numerical investigation on the effects of bedding plane properties on the mechanical and acoustic emission characteristics of sandy mudstone | |
| Wu et al. | Development and verification of three-dimensional equivalent discrete fracture network modelling based on the finite element method | |
| Verdon et al. | Examining the capability of statistical models to mitigate induced seismicity during hydraulic fracturing of shale gas reservoirs | |
| Ispas et al. | Laboratory testing and numerical modelling of fracture propagation from deviated wells in poorly consolidated formations | |
| Das et al. | A mesh-free approach for fracture modelling of gravity dams under earthquake | |
| Li et al. | Survey of measurement-while-drilling technology for small-diameter drilling machines | |
| CN105467429B (en) | A kind of method of hard-rock mine mountainous region pressure coupling prediction | |
| Yang et al. | Numerical simulation on slope stability analysis considering anisotropic properties of layered fractured rocks: a case study | |
| Haustveit et al. | Empirical Meets Analytical-Novel Case Study Quantifies Fracture Stress Shadowing and Net Pressure Using Bottom Hole Pressure and Optical Fiber | |
| Jiang et al. | Rock crack type identification by Gaussian process learning on acoustic emission | |
| AbuAisha et al. | Hydro-mechanically coupled FDEM framework to investigate near–wellbore hydraulic fracturing in homogeneous and fractured rock formations | |
| CN117970482B (en) | Method, system, medium and terminal for judging microseism event dry-wet type | |
| Patil et al. | Numerical prediction of ice rubble field loads on the Norströmsgrund lighthouse using cohesive element formulation | |
| Tang et al. | Discrete element simulation for investigating fragmentation mechanism of hard rock under ultrasonic vibration loading | |
| Lou et al. | Study on fracture mechanism of hydraulic fracturing in sandstone by acoustic emission parameters | |
| Brown | Seismic hazard evaluation using apparent stress ratio for mining-induced seismic events | |
| CN117309639A (en) | Method and equipment for confirming stress critical value of coal body at rock burst deep base point | |
| CN104792965B (en) | Wall rock loosening ring method of testing based on boring energy | |
| Halleck | Recent advances in understanding perforator penetration and flow performance | |
| CN113885073B (en) | Rock fracture and fracturing acoustic emission fine analysis method based on original waveform | |
| Luo et al. | Application of acoustic emission on crushing monitoring of individual soil particles in uniaxial compression test | |
| Wang et al. | Dynamic crack growth characteristics of half-circular disk PMMA specimens with prefabricated penetration defects under impact loading | |
| Liu et al. | A crushing index for granular soils based on comminution energy consumption theory | |
| CN113982571A (en) | Evaluation method for tight reservoir fracturing affected zone |
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 |