CN115184990A - Microseism monitoring and observing method - Google Patents

Microseism monitoring and observing method Download PDF

Info

Publication number
CN115184990A
CN115184990A CN202210840925.1A CN202210840925A CN115184990A CN 115184990 A CN115184990 A CN 115184990A CN 202210840925 A CN202210840925 A CN 202210840925A CN 115184990 A CN115184990 A CN 115184990A
Authority
CN
China
Prior art keywords
monitoring
observation
well
effective receiving
ground
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
CN202210840925.1A
Other languages
Chinese (zh)
Other versions
CN115184990B (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.)
Jianghan Branch Of Sinopec Petroleum Engineering Geophysics Co ltd
Oil & Gas Survey Cgs
Original Assignee
Jianghan Branch Of Sinopec Petroleum Engineering Geophysics Co ltd
Oil & Gas Survey Cgs
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 Jianghan Branch Of Sinopec Petroleum Engineering Geophysics Co ltd, Oil & Gas Survey Cgs filed Critical Jianghan Branch Of Sinopec Petroleum Engineering Geophysics Co ltd
Priority to CN202210840925.1A priority Critical patent/CN115184990B/en
Publication of CN115184990A publication Critical patent/CN115184990A/en
Application granted granted Critical
Publication of CN115184990B publication Critical patent/CN115184990B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/288Event detection in seismic signals, e.g. microseismics
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

The invention discloses a microseism monitoring and observing method, which relates to the technical field of microseism monitoring and observing, and is characterized in that according to geological parameters such as stratum inclination angle, interface depth, longitudinal and transverse wave velocity, density, absorption factor Q and the like, engineering parameters such as the relative position of a fracturing well monitoring well, well track and the like, through forward modeling, wave front diffusion compensation, stratum absorption attenuation, transmission coefficients and the like are combined, and the direct wave amplitude of a ground observation point and the first arrival wave of an observation point in the monitoring well are obtained; determining the effective receiving range of the ground detector by taking a preset threshold as a critical point according to the amplitude curve of the direct wave of the ground observation point; determining an effective receiving area of a geophone in a well according to the distribution condition of the first-motion waves; and arranging the detectors according to the effective receiving range of the ground detectors and/or the effective receiving area of the borehole detectors to carry out microseismic monitoring observation. The invention can simultaneously or respectively carry out micro-seismic ground observation and micro-seismic well observation, and improves the utilization rate of the geophone by reasonably designing the position of the geophone.

Description

Microseism monitoring and observing method
Technical Field
The invention relates to the technical field of microseism monitoring observation, in particular to a microseism monitoring observation method.
Background
The micro-seismic monitoring technology is a new geophysical technology developed in the 20 th century, and is a geophysical technology for monitoring the influence, effect and underground state of production activities by observing and analyzing micro-seismic events generated in the production activities. The method can be applied to a plurality of fields such as oil gas development, coal mine monitoring, mine pressure monitoring, geological disaster monitoring and the like.
According to the arrangement position of the geophone, the microseism monitoring technology is divided into borehole microseism monitoring and ground microseism monitoring, when the geophone is arranged, only two parameters of the depth of a target layer and the length of a horizontal well section are usually referred, and factors influencing seismic wave propagation distance and energy such as stratum speed, dip angle and absorption attenuation are not considered, so that more ineffective geophones are generated, resource waste is caused, and meanwhile, the detection accuracy is influenced.
Therefore, how to improve the effectiveness of the detector setting, improve the utilization rate of the detector, and more scientifically perform microseism monitoring observation is a problem that needs to be solved urgently by the technical personnel in the field.
Disclosure of Invention
In view of the above, the present invention provides a microseism monitoring and observation method to overcome the above technical problems.
In order to achieve the above purpose, the invention provides the following technical scheme:
a microseism monitoring observation method comprises the following steps:
collecting geological parameters and engineering parameters;
according to geological parameters and engineering parameters, a work area geological model is established, the position of a micro-seismic excitation point is determined, and the effective receiving range of a ground detector and/or the effective receiving area of a borehole detector are/is determined through forward simulation;
and setting the detectors to carry out microseism monitoring observation according to the effective receiving range of the ground detectors and/or the effective receiving area of the borehole detectors.
Optionally, the method for determining the effective receiving range of the ground detector comprises: according to geological parameters and engineering parameters, a work area geological model is established, the position of a micro-seismic excitation point is determined, and through forward modeling, wave front diffusion compensation, stratum absorption attenuation and transmission coefficients are combined to obtain the direct wave amplitude of a ground observation point; and determining the effective receiving range of the ground detector by taking a preset threshold as a critical point according to the amplitude curve of the direct wave of the ground observation point.
Optionally, the method for determining the effective receiving area of the borehole detector comprises: according to geological parameters and engineering parameters, establishing a work area geological model, determining the position of a micro-seismic excitation point and the relative position of an observation point and the micro-seismic excitation point in a monitoring well, and acquiring a first-arrival wave of the observation point in the monitoring well through forward simulation; and determining the effective receiving area of the detector in the well according to the distribution of the first arrival waves.
Optionally, the geological parameters include a formation dip angle, an interface depth, a longitudinal and transverse wave velocity, a density, an absorption factor Q, and the like, and the engineering parameters include a relative position of a fracturing well monitoring well, a well track, and the like.
Through the geological parameters, factors such as the wave front diffusion effect of seismic waves, stratum absorption attenuation, the polarization characteristic of the seismic waves, stratum inclination angles and the like are considered during forward simulation, and the amplitude energy can be reflected more accurately.
Optionally, the preset threshold is 50% of the energy of the microseism excitation point, and a range in which the amplitude of the direct wave of the ground observation point is greater than 50% of the energy of the microseism excitation point is determined as an effective receiving range of the ground detector.
Optionally, the forward simulation is performed by ray tracing.
Optionally, the method for determining the effective receiving area of the borehole geophone according to the distribution of the first-motion waves comprises: and determining the region taking the initial wave as the direct wave as the effective receiving region of the borehole geophone.
Optionally, the detectors are arranged in the effective receiving range of the ground detector and/or the effective receiving area of the borehole detector at equal intervals, the detection signals of the detectors are remotely transmitted to the terminal, the positions of the micro-seismic event points are obtained through data processing, and the positions are three-dimensionally displayed, so that micro-seismic monitoring and observation are realized.
According to the technical scheme, the invention discloses and provides a microseism monitoring and observing method, and compared with the prior art, the microseism monitoring and observing method has the following beneficial effects:
(1) The invention can determine the effective receiving length of the ground micro-seismic monitoring survey line in different directions through forward simulation according to a geological model (comprising parameters such as a stratum inclination angle, an interface depth, a longitudinal wave velocity, a density and an absorption factor Q), forms a more targeted observation effect, breaks through the current situation of mechanical uniform survey line distribution, enables the observation effect to be more consistent with the actual geological model, and improves the utilization rate of the detector.
(2) When the micro-earthquake monitoring and observation device is used for monitoring and observing the micro-earthquake, the influence of various parameters on the earthquake wave propagation process is fully considered, the longitudinal and transverse wave propagation rules are simulated, the effective receiving range is determined, and the geophones are arranged in the effective range, so that on one hand, the geophone waste is avoided, and on the other hand, the interference of the geophones on positioning is reduced.
(3) The invention can simultaneously or respectively carry out micro-seismic ground observation and micro-seismic well observation, and provides a basis for the perfection of the micro-seismic monitoring construction design.
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 or the prior art will be briefly described below, it is obvious that the drawings in the following description are only embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the provided drawings without creative efforts.
FIG. 1 is a schematic diagram showing the attenuation of seismic wave energy in example 1;
FIG. 2 (a) is a schematic diagram of total reflection of an inclined layered medium in a short range of a detector arrangement in example 2;
FIG. 2 (b) is a schematic diagram of total reflection of a tilted layered medium in the case of a longer range of the detector arrangement in example 2;
FIG. 3 is a schematic diagram showing the attenuation of seismic wave energy in a dipping formation in example 3;
FIG. 4 is a graph showing the ray tracing result at the nearest monitoring distance;
FIG. 5 is a schematic diagram of a forward modeling result of a waveform at a closest monitoring distance;
FIG. 6 is a diagram illustrating ray tracing results at the farthest monitoring distance;
FIG. 7 is a schematic diagram of a waveform forward modeling result at the farthest monitoring distance;
FIG. 8 is a software function interface diagram;
FIG. 9 is a diagram of a preset interface for micro-seismic ground observation;
FIG. 10 is a ray tracing and ray energy diagram for a microseismic ground survey;
FIG. 11 is a schematic view of observation and ray analysis in a microseismic well;
FIG. 12 is a graph showing the results of single interface ray energy analysis;
FIG. 13 is a schematic representation of the method steps of the present invention;
FIG. 14 is a Fermat principle schematic;
FIG. 15 is a schematic diagram of the variation of the wavefront in a homogeneous medium.
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 obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
The embodiment of the invention discloses a microseism monitoring and observing method, which is shown in figure 13 and has the following design principle:
1. and collecting and determining geological parameters and engineering parameters.
The geological parameters comprise a stratum inclination angle, an interface depth, a longitudinal wave speed, a transverse wave speed, a density, an absorption factor Q and the like, and the engineering parameters comprise a relative position of a fracturing well monitoring well, a well track and the like.
2. And determining the effective receiving range of the ground detector.
A geological model of a work area is built, a stratum inclination angle, an interface depth, a longitudinal wave velocity, a transverse wave velocity, a density, an absorption factor Q and a relative position of a fracturing well monitoring well are input, forward modeling is carried out through ray tracing, wave front diffusion compensation, stratum absorption attenuation, a transmission coefficient and the like are considered, and the direct wave amplitude of a ground observation point is obtained; determining the effective receiving range of the ground detector by taking a preset threshold as a critical point according to the amplitude curve of the direct wave of the ground observation point;
the specific method for performing forward modeling through ray tracing comprises the following steps:
ray tracing is to obtain the seismic wave propagation path and propagation time from the excitation point to the demodulator probe according to the geologic model under the condition of setting the positions of the excitation point and the demodulator probe. The calculation method includes a pilot shot method and a bending method.
The trial injection method is based on Snell's law: assuming that the longitudinal and transverse wave velocities of different layers are respectively expressed as upsilon P1 ,υ S1 ,υ P2 ,υ S2 ,···,υ Pi ,υ Si The angle of incidence of the wave is represented by θ P1 ,θ S1 ,θ P2 ,θ S2 ,···,θ Pi ,θ Si Expressed, then snell's law can be expressed as:
Figure BDA0003751007130000051
wherein P is a ray parameter.
The trial injection method determines a ray path by modifying the emergent direction of rays at an excitation point and calculating a proper incident angle at which the emergent rays can reach a detection point.
The bending method is based on the fermat principle: the propagation path of a wave in various media satisfies the condition that the time taken is the shortest. According to the method, the initial path connecting the excitation point and the receiving point is disturbed, the intersection point coordinates of the ray path and each stratum interface are continuously and iteratively modified until the minimum travel time criterion is met, and the calculation efficiency of the method is improved.
The Fermat principle is as follows: in fig. 14, the solid line represents the actual propagation path of the seismic wave, the dotted line represents other imaginary paths, and the solid line represents the shortest propagation time among all paths.
1. Wave front diffusion compensation:
assuming that the medium is uniform, when the seismic wave propagates in the underground medium, as the propagation time increases, the wavefront surface can be regarded as a spherical surface with the seismic source as the center, which is continuously enlarged, so that under the condition that the total energy emitted by the seismic source is not changed, the energy per unit area on the wavefront surface is continuously reduced, and the amplitude of the received seismic wave is continuously reduced.
As shown in fig. 15, assuming the total energy of the source is E, the energy density on the wavefront surface at any time is:
Figure BDA0003751007130000061
wherein E is the total energy; r-wave propagation distance; v-the propagation velocity of the wave; t-the propagation time of the wave.
According to equation 2, the energy density of the wavefront at r =1 is:
Figure BDA0003751007130000062
from equation 2 and equation 3, equation 4 can be derived, which can represent the wavefront diffusion energy attenuation ratio:
Figure BDA0003751007130000063
2. absorption and attenuation of stratum:
since the formation medium is usually not completely elastic, in addition to the attenuation of the amplitude of the seismic wave due to wave front diffusion, the absorption of the formation also causes the amplitude of the seismic wave to be exponentially attenuated with increasing propagation distance. Assume that the initial amplitude produced when the seismic source excites a seismic wave is A 0 Where A is the amplitude of the seismic wave at a distance r from the seismic source, and α is the absorption coefficient of the medium, then:
A=A 0 e -αr (formula 5)
The amplitude attenuation factor due to formation absorption is:
Figure BDA0003751007130000064
in the formula, t is the travel time of the seismic wave propagation distance r; beta-the attenuation coefficient of the medium.
β = α ν (formula 7)
Where ν is the propagation velocity of seismic waves in the medium.
In the process of seismic data processing, the attenuation energy of the stratum to the seismic wave is described by using a quality factor Q, where the quality factor represents the ratio of the originally stored energy E to the consumed energy Δ E after the seismic wave propagates by a wavelength λ, that is:
Figure BDA0003751007130000071
and (3) removing high-order terms after the above formula is expanded to obtain the relation between the quality factor Q and the absorption coefficient:
Figure BDA0003751007130000072
in the formula, f represents the frequency of seismic waves;
the relationship between the quality factor and the attenuation factor can be obtained from equation 6 and equation 9:
Figure BDA0003751007130000073
thus in a non-fully elastic medium, the high frequency components of the seismic waves attenuate faster than the low frequency components.
3. Transmission coefficient:
according to the Zoeppritz equation, the reflection coefficient and the projection coefficient are related to the formation velocity, the density, the incident angle, and the like as follows:
Figure BDA0003751007130000074
wherein the content of the first and second substances,
Figure BDA0003751007130000075
in the formula, R PP ,R PS ,T PP ,T PS Reflection coefficients and transmission coefficients of reflected P-waves, reflected SV-waves, transmitted P-waves and transmitted SV-waves, respectively, expressed in displacement amplitudes; b is 1 ,B 2 ,B 3 ,B 4 ,B 5 Is the displacement amplitude of each wave. Rho 1 ,ρ 2 ,υ P1 ,υ P2 ,υ S1 ,υ S2 The density and the longitudinal and transverse wave propagation speeds of the medium are shown, alpha represents an incident angle, beta represents a reflection angle, and alpha And beta Is the transmission angle. According to snellThe reflection angle and the transmission angle can be expressed as a function of the incident angle α according to the law of 'er's law, and thus the reflection coefficient and the transmission coefficient are actually related to the incident angle, as well as the medium density and the velocity of the longitudinal and transverse waves.
In the following, taking an inclined stratum model as an example, the observation optimization is performed by considering factors such as a total reflection critical point and stratum attenuation, and the like as follows.
Example 1: observation optimization with formation attenuation factor as constraint
Due to different propagation paths of seismic waves, micro-seismic signals received by the ground detectors are attenuated to different degrees, and the amplitude energy received by each detector is also obviously different. When the attenuation degree of the amplitude energy is large, the detector cannot monitor an effective signal easily, and a plurality of factors causing the amplitude attenuation exist.
In this embodiment, as shown in fig. 1, a red explosion symbol indicates a fracture position of a destination layer, that is, a position of an excitation point, and an arrangement design in different directions is performed by using a projection of the position on the ground as a coordinate origin, where both sides of the design origin in fig. 1 have an arrangement length of 2000m, and energy of the excitation point is 1, and an amplitude of a direct wave received by each receiving point is calculated. When the wave front diffusion effect of the seismic waves is considered, the green curve in fig. 1 reflects the amplitude energy (the coordinate axis marks the amplitude) of the direct waves of different receiving points; after adding absorption attenuation to the formation, the amplitude energy is further reduced (see blue curve in fig. 1); further considering the polarization characteristics of the seismic wave, a red hyperbola is obtained through simulation, and the amplitude energy can be reflected more objectively. In a specific embodiment, the signal cannot be effectively received and identified when the amplitude is determined to be attenuated to half of the initial energy, and the effective receiving length of the array is determined to be about 1400m through numerical simulation of the geological model.
Example 2: observing and optimizing by taking total reflection critical point as constraint
When a surface micro-seismic monitoring design of a certain well is carried out, the embodiment determines the arrangement monitoring length to be 3200m according to the design criteria required by technical regulations. However, the shallow velocity of the well is lower than the deep velocity, the signal is propagated from the deep layer to the shallow interface and may be totally reflected at a certain critical point, and the arrangement beyond the critical point cannot receive the direct wave transmission signal and is regarded as invalid arrangement in microseism monitoring.
When a crack occurs, if the detector array is short, the ray emergence angle is small, and all the arranged detectors can receive microseismic signals (see fig. 2 (a)); if the alignment length is increased, total reflection occurs in the formation in the direction of the formation inclination (see fig. 2 (b)). Through ray tracing, it is determined that the formation updip direction can receive direct wave transmission signals within the range of 2800m of the ground surface projection point, no effective signals can be monitored within the range of 2900-3200m in an arrangement mode, and detector waste is avoided.
Example 3: observation optimization by taking stratum inclination angle as constraint
When the inclination angle of the stratum is large (see fig. 3), the difference between the energy of the signals received by the detectors at two sides perpendicular to the structure trend is large, and under the condition of the same offset distance, the energy received in the upward inclination direction is far smaller than that in the downward inclination direction. Under the geological model of FIG. 3, the amplitude energy received by a detector with a formation dip offset of 2000m is comparable to the amplitude energy received at 1600m dip. Therefore, when the design is carried out, the arrangement length can be distributed according to the proportion according to the amplitude energy change characteristics, and the ratio of the arrangement length in the ascending direction to the arrangement length in the descending direction from the central point is set to be 4:5.
3. An effective receiving area of the borehole geophone is determined.
In the process of microseismic monitoring in a well, the closer the detector is to a target layer, the better the monitoring effect is generally considered. However, in some specific geological environments, direct waves and refracted waves monitored by a detector close to a target layer interfere with each other, accurate picking up of the first arrival of the direct waves cannot be achieved, and positioning accuracy of a micro seismic event in a well is reduced, so that in-well observation optimization is performed by taking a development range of the refracted waves as a constraint, and detection efficiency can be effectively improved.
The shale gas exploration block target layer takes shale as a main part, the speed is low, high-speed layers are arranged above and below the target layer, when the monitoring distance is large, refracted waves can be generated on a high-speed layer interface, the first-arrival waves received by the detector are refracted waves, the direct-arrival first-arrival waves are difficult to pick up, and micro-seismic event positioning cannot be carried out.
Fig. 4 is a diagram illustrating the ray tracing result at the nearest monitoring distance. In fig. 4, the red solid line is the fractured well track, the black dots represent the position of each detector, and ray tracing simulation is carried out on the primary waves of the detectors at different depths. Simulation results show that when the target point A (the target point A is a fracture position closest to the detector) is excited, (1) the primary wave received by the detector in the area is a direct wave, and (2) the primary wave received by the detector in the area (near a target layer) is a refracted wave. Fig. 5 shows that the target point a is excited, and seismic waves received by the geophone are forward recorded, (1) direct waves in the area are first arrival waves, and the direct waves can be accurately picked up when arriving for positioning the micro-seismic event, but (2) refracted waves and the direct waves in the high-speed layer in the area are mixed up, so that the first arrival waves are difficult to pick up, and the positioning precision of the micro-seismic event is influenced, so that the positioning precision of the micro-seismic event is influenced.
This phenomenon is more obvious as the monitoring distance increases, and the development range of the refracted wave increases after the monitoring distance increases, as shown in fig. 6 and 7, which are the ray tracing result and the waveform forward simulation diagram at the farthest monitoring distance, respectively. Simulation results show that a target point B (a fracturing position farthest from a detector) performs seismic source excitation, and a primary wave received by the detector in the area (1) is a direct wave; (2) direct waves received by the area detector interfere with refracted wave waveforms, and after reverse superposition, signals are weakened and cannot be picked up in first arrivals; (3) the first arrival wave received by the area detector is a refracted wave, and direct wave pickup is difficult. Under such a monitoring distance, the closer the detector to the target layer is, the more difficult it is to pick up reliable direct wave first arrival time, and the positioning of the microseism event cannot be realized, so that in order to ensure that the microseism signal received by the target point B can be accurately positioned, the depth of the detector cannot exceed the area (1).
4. And carrying out microseism monitoring observation.
According to the position of the micro-seismic excitation point, the effective receiving range of the ground detector and the effective receiving area of the borehole detector, the detectors are arranged in the effective receiving range of the ground detector and the effective receiving area of the borehole detector at equal intervals or in other modes, then the detection signals of the detectors are remotely transmitted to a terminal and are displayed in a three-dimensional mode, and the micro-seismic monitoring observation is realized, referring to fig. 8.
The software of the microseism monitoring and observing system is designed according to the method, the basic functions of the software are shown in figures 8-13, and functions of displaying microseism events in segments and whole segments, self-media demonstration and the like are integrated in the software.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
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. A microseism monitoring and observing method is characterized by comprising the following steps:
collecting geological parameters and engineering parameters;
according to geological parameters and engineering parameters, a work area geological model is established, the position of a micro-seismic excitation point is determined, and the effective receiving range of a ground detector and/or the effective receiving area of a borehole detector are/is determined through forward simulation;
and setting the detectors to carry out microseism monitoring observation according to the effective receiving range of the ground detectors and/or the effective receiving area of the borehole detectors.
2. The microseism monitoring observation method of claim 1, wherein the method for determining the effective receiving range of the ground geophone comprises the following steps:
according to geological parameters and engineering parameters, a work area geological model is established, the position of a microseism excitation point is determined, and through forward modeling, wave front diffusion compensation, stratum absorption attenuation and transmission coefficients are combined to obtain the direct wave amplitude of a ground observation point; and determining the effective receiving range of the ground detector by taking a preset threshold as a critical point according to the amplitude curve of the direct wave of the ground observation point.
3. A microseismic monitoring observation method as set forth in claim 1 wherein the method of determining the effective receiving area of the geophone in a well is:
according to the geological parameters and the engineering parameters, a work area geological model is established, the position of a micro-seismic excitation point and the relative position of an observation point and the micro-seismic excitation point in a monitoring well are determined, and a first arrival wave of the observation point in the monitoring well is obtained through forward simulation; and determining the effective receiving area of the geophone in the well according to the distribution situation of the first-motion waves.
4. The microseism monitoring observation method of claim 1, wherein the geological parameters comprise formation dip, interface depth, longitudinal and transverse wave velocity, density and absorption factor Q, and the engineering parameters comprise relative position of a fracturing well monitoring well and well track.
5. The microseism monitoring and observation method according to claim 2, wherein the preset threshold is 50% of the energy of the microseism excitation point, and the range in which the direct wave amplitude of the ground observation point is greater than 50% of the energy of the microseism excitation point is determined as the effective receiving range of the ground geophone.
6. A microseismic monitoring observation method according to claim 1 wherein the forward modeling is performed by ray tracing.
7. A microseismic monitoring observation method according to claim 3 wherein the method of determining the effective receiving area of the geophone in the well based on the distribution of the first-arrival waves comprises: and determining the region of the first-motion wave as the direct wave as the effective receiving region of the borehole geophone.
8. The microseism monitoring and observation method according to claim 1, wherein the geophones are arranged at equal distances in the effective receiving range of the ground geophone and/or the effective receiving area of the geophone in a well, the detection signal of the geophone is remotely transmitted to a terminal, the position of a microseism event point is obtained, and the position is three-dimensionally displayed to realize microseism monitoring and observation.
CN202210840925.1A 2022-07-18 2022-07-18 Microseism monitoring and observing method Active CN115184990B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210840925.1A CN115184990B (en) 2022-07-18 2022-07-18 Microseism monitoring and observing method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210840925.1A CN115184990B (en) 2022-07-18 2022-07-18 Microseism monitoring and observing method

Publications (2)

Publication Number Publication Date
CN115184990A true CN115184990A (en) 2022-10-14
CN115184990B CN115184990B (en) 2023-03-10

Family

ID=83520238

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210840925.1A Active CN115184990B (en) 2022-07-18 2022-07-18 Microseism monitoring and observing method

Country Status (1)

Country Link
CN (1) CN115184990B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9225984D0 (en) * 1992-12-12 1993-02-10 Schlumberger Ltd Method for improving signal to noise ratio from a geophone array
US6182012B1 (en) * 1998-02-20 2001-01-30 Institut Francais Du Petrole Method and device intended for permanent monitoring of an underground formation
US20040037165A1 (en) * 2002-08-21 2004-02-26 Erwann Lemenager Method of accurately determining positions of deployed seismic geophones
CN104950327A (en) * 2015-06-30 2015-09-30 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for determining positions of geophones of ground microseismic observation system
CN105572725A (en) * 2016-02-17 2016-05-11 西南石油大学 Distribution design method of ground microseism monitoring station
CN106054239A (en) * 2016-05-23 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Microseismic pressing crack monitoring observation method
CN109655918A (en) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 Ground shallow well micro-seismic monitoring measuring platform station location determines method and system
CN112415601A (en) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 Method and device for determining surface quality factor Q value

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB9225984D0 (en) * 1992-12-12 1993-02-10 Schlumberger Ltd Method for improving signal to noise ratio from a geophone array
US6182012B1 (en) * 1998-02-20 2001-01-30 Institut Francais Du Petrole Method and device intended for permanent monitoring of an underground formation
US20040037165A1 (en) * 2002-08-21 2004-02-26 Erwann Lemenager Method of accurately determining positions of deployed seismic geophones
CN104950327A (en) * 2015-06-30 2015-09-30 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Method for determining positions of geophones of ground microseismic observation system
CN105572725A (en) * 2016-02-17 2016-05-11 西南石油大学 Distribution design method of ground microseism monitoring station
CN106054239A (en) * 2016-05-23 2016-10-26 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Microseismic pressing crack monitoring observation method
CN109655918A (en) * 2017-10-11 2019-04-19 中国石油化工股份有限公司 Ground shallow well micro-seismic monitoring measuring platform station location determines method and system
CN112415601A (en) * 2020-11-03 2021-02-26 中国石油天然气集团有限公司 Method and device for determining surface quality factor Q value

Also Published As

Publication number Publication date
CN115184990B (en) 2023-03-10

Similar Documents

Publication Publication Date Title
RU2319982C2 (en) Definition of characteristics of inclined cracks with the aid of exchanged reflected waves in combination with the process of consistent exclusion of influence of the superposed layers
Haldorsen et al. Borehole acoustic waves
CN102282481B (en) Data acquisition and prestack migration based on seismic visibility analysis
US7872944B2 (en) Method of monitoring microseismic events
US8072840B2 (en) Fracture clusters identification
US6832161B1 (en) Method for attenuating water layer multiples
EP3433643B1 (en) Method and device for estimating sonic slowness in a subterranean formation
US10670754B2 (en) System and method for processing microseismic data
US8045417B2 (en) Analyzing 2-D surface and/or borehole seismic data to locate subsurface diffractors
CN102985850B (en) Use passive seismic data acquisition and the processing of multiple level sensor array
US9329293B2 (en) Systems and methods for removing acquisition related effects from seismic data
US20170371057A1 (en) Method of and system for creating a seismic profile
CN105277971A (en) Micro-seismic monitoring system and method
US20100312480A1 (en) Method for monitoring fluid flow in a multi-layered system
Luo et al. Seismic inversion of shale reservoir properties using microseismic-induced guided waves recorded by distributed acoustic sensing
Langet et al. Joint focal mechanism inversion using downhole and surface monitoring at the Decatur, Illinois, CO2 injection site
US20150268365A1 (en) Method to characterize geological formations using secondary source seismic data
CN110687607B (en) Stoneley wave detection method and system
CN115184990B (en) Microseism monitoring and observing method
CN110579798B (en) Seismic acquisition observation method and system with equal reflection angle intervals
CN110579799B (en) Seismic acquisition observation method and system with equal travel time intervals
US11435490B2 (en) Seismic surveys using two-way virtual source redatuming
US20140286125A1 (en) Seismic acquisition method and apparatus
Djikpesse C 13 and Thomsen anisotropic parameter distributions for hydraulic fracture monitoring
US20230273333A1 (en) Simultaneous common-offset migration

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