CN117637064A - Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal - Google Patents
Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal Download PDFInfo
- Publication number
- CN117637064A CN117637064A CN202311616181.6A CN202311616181A CN117637064A CN 117637064 A CN117637064 A CN 117637064A CN 202311616181 A CN202311616181 A CN 202311616181A CN 117637064 A CN117637064 A CN 117637064A
- Authority
- CN
- China
- Prior art keywords
- point
- seismic
- reinjection
- points
- micro
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000009826 distribution Methods 0.000 title claims abstract description 104
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 71
- 238000004458 analytical method Methods 0.000 title claims abstract description 38
- 230000000007 visual effect Effects 0.000 title claims abstract description 38
- 230000035939 shock Effects 0.000 claims abstract description 52
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000005457 optimization Methods 0.000 claims abstract description 19
- 238000012544 monitoring process Methods 0.000 claims description 22
- 230000006870 function Effects 0.000 claims description 16
- 238000004590 computer program Methods 0.000 claims description 9
- 230000008569 process Effects 0.000 claims description 8
- 238000005516 engineering process Methods 0.000 claims description 4
- 238000012800 visualization Methods 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- 238000012545 processing Methods 0.000 abstract description 6
- 238000010586 diagram Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 3
- 230000007547 defect Effects 0.000 description 2
- 238000002372 labelling Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- 238000003491 array Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000003672 gas field water Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/80—Data visualisation
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
- G06Q10/047—Optimisation of routes or paths, e.g. travelling salesman problem
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Chemical & Material Sciences (AREA)
- Theoretical Computer Science (AREA)
- Computing Systems (AREA)
- Economics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Strategic Management (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Crystallography & Structural Chemistry (AREA)
- Physics & Mathematics (AREA)
- Operations Research (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Data Mining & Analysis (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a reinjection well reinjection water distribution path visual analysis method, a reinjection well reinjection water distribution path visual analysis system and a reinjection well reinjection water distribution path visual analysis terminal, which relate to the technical field of data processing and have the technical scheme that: acquiring space positioning information of each micro-seismic point in a target area, and establishing an initial seismic point space model; taking a plurality of microseismic points which are in continuous distribution trend as effective shock points, and taking the rest microseismic points as ineffective shock points; carrying out distributed path optimization on all effective shock points in the initial shock point space model through a path optimization function; and visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model. According to the method, the distribution distance condition and the distribution direction condition among the micro-seismic points are considered, and a path optimization function is adopted to conduct distribution path optimization on all the effective seismic points in the initial seismic point space model, so that a reinjection well reinjection water distribution path capable of simulating real conditions in detail is obtained, and accurate analysis on reinjection effect of reinjection well reinjection water is facilitated.
Description
Technical Field
The invention relates to the technical field of data processing, in particular to a reinjection well reinjection water distribution path visual analysis method, a reinjection well reinjection water distribution path visual analysis system and a reinjection well reinjection water distribution terminal.
Background
Reinjection well reinjection water is an environmental protection measure for reinjecting oil and gas field produced water into a stratum through the reinjection well, so that stratum energy is supplemented, stratum pressure is kept, and environmental pollution is reduced. The distribution of reinjection well reinjection water in the stratum is monitored, so that the reinjection effect is effectively known.
The method mainly adopts microseism monitoring technology to collect the positions of microseism points, a plurality of monitoring substations distributed in arrays are arranged around a shaft of a reinjection well, the monitoring substations are utilized to receive the time difference of microseism waves emitted by the microseism points, and an inversion positioning method is utilized to solve the three-dimensional space coordinates of the microseism points, so that the overall distribution condition of the microseism points in the three-dimensional space is obtained. At present, the prior art records that gas field water waves and ranges below the surface of a reinjection well are analyzed according to the overall distribution condition of a plurality of microseismic points in a three-dimensional space, but the method cannot obtain the specific flow condition of reinjection water of the reinjection well in a stratum, and cannot accurately analyze the reinjection effect of reinjection water of the reinjection well, such as blockage condition.
Therefore, how to study and design a visual analysis method, a visual analysis system and a visual analysis terminal for a reinjection water distribution path of a reinjection well, which can overcome the defects, are the problems which are urgently needed to be solved at present.
Disclosure of Invention
In order to solve the defects in the prior art, the invention aims to provide a reinjection well reinjection water distribution path visual analysis method, a reinjection well reinjection water distribution path visual analysis system and a reinjection terminal, partial non-reinjection well reinjection water caused microseismic points are filtered in consideration of distribution distance conditions and distribution direction conditions among microseismic points, and then a path optimization function is adopted to conduct distribution path optimization on all effective seismic points in an initial seismic point space model, so that a reinjection well reinjection water distribution path capable of simulating real conditions in detail is obtained, and the reinjection well reinjection water distribution path is subjected to visual labeling in the initial seismic point space model, so that accurate analysis on reinjection effect of reinjection well reinjection water is facilitated.
The technical aim of the invention is realized by the following technical scheme:
in a first aspect, a method for visual analysis of a reinjection water distribution path of a reinjection well is provided, comprising the steps of:
acquiring the space positioning information of each micro-seismic point in the target area, and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point;
traversing and searching an initial seismic point space model around by taking a shaft as a monitoring center, taking a plurality of microseismic points with continuous distribution trend as effective seismic points, and taking the rest microseismic points as ineffective seismic points;
carrying out distribution path optimization on all effective shock points in the initial shock point space model through a path optimization function to obtain a reinjection water distribution path of the reinjection well;
and visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
Further, the spatial positioning information of the microseism points is acquired by utilizing a microseism monitoring technology, and the specific process is as follows:
a well bore is taken as a monitoring center to determine a target area, and a plurality of monitoring substations distributed in an array are arranged in the target area;
collecting the time when a plurality of monitoring substations receive the micro-seismic waves sent by the micro-seismic points;
and establishing an equation set according to the time difference of the micro seismic waves sent by the micro seismic points received by different monitoring substations, and solving the three-dimensional space coordinates of the micro seismic points by an inversion positioning method.
Further, the difference of the micro-seismic wave velocities corresponding to stratum mediums at different points in the target area is not larger than a wave velocity threshold;
the method comprises the steps of solving the micro-seismic wave velocity applied in the three-dimensional space coordinates of the micro-seismic points by an inversion positioning method, wherein the micro-seismic wave velocity is the average value of the micro-seismic wave velocity corresponding to stratum media of each point in a target area.
Further, the process of traversing and searching the initial seismic point space model around by taking the shaft as a monitoring center specifically comprises the following steps:
determining a search radius and a search bias angle of the traversal search;
selecting a micro-seismic point with a distance smaller than the searching radius and a searching path between the micro-seismic point B and the micro-seismic point B to obtain a seismic point set;
if the distance between the micro-seismic point C and the micro-seismic point B is smaller than the searching radius, and one micro-seismic point A is arbitrarily selected from the seismic point set, the following conditions are satisfied: if the space included angle between the first vector of the micro-seismic point A pointing to the micro-seismic point B and the second vector of the micro-seismic point B pointing to the micro-seismic point C does not exceed the search offset angle, the search radius exists between the micro-seismic point C and the micro-seismic point B, and the micro-seismic point C is taken as an effective seismic point;
and the micro-seismic points with the distance from the monitoring center smaller than the searching radius are taken as effective seismic points.
Further, the search bias angle is not more than 90 °.
Further, the searching radius is not smaller than the plane distribution interval of the microseismic points in the target area;
the plane distribution interval is specifically as follows:
calculating to obtain a unit area according to the ratio of the area of the horizontal plane in the target area to the number of all acquired microseismic points;
and establishing an equivalent circle with the area equal to the unit area, and taking the diameter of the equivalent circle as the plane distribution interval.
Further, the expression of the path optimization function is specifically:
wherein,representing the jth effective shock point in the kth path; />Representing the jth effective shock point in the (k+1) th path; e represents a set of all effective shock points in the initial shock point space model; d (D) i Representing the ith effective seismic point randomly selected from the set E of effective seismic points, and the effective seismic point D i And effective shock point->A search path exists between the two;representing the effective vibration point->Relative to the effective shock point D i And (D) represents the effective shock point D i As effective shock point->The effective vibration point D is selected as the starting point of the flow direction when the offset is calculated i-1 Representing the effective shock point D i The effective vibration points which are selected as the flow direction starting points when the offset is calculated; n (N) k Representing the number of effective shock points in the kth path; m represents the number of paths in the distribution path; and n represents an intersection symbol.
Further, the calculation formula of the offset is specifically:
β＝(x i -x i-1 ，y i -y i-1 ，z i -z i-1 )
wherein,representing the effective vibration point->Is a three-dimensional space coordinate of (2); (x) i ，y i ，z i ) Representing the effective shock point D i Is a three-dimensional space coordinate of (2); (x) i-1 ，y i-1 ，z i-1 ) Representing the effective shock point D i-1 Is a three-dimensional space coordinate of (2); alpha represents the effective shock point->And the effective vibration point D i A space vector between; beta represents the effective shock point D i And the effective vibration point D i-1 Spatial vector between them.
In a second aspect, a reinjection well reinjection water distribution path visualization analysis system is provided, comprising:
the monitoring and positioning module is used for acquiring the space positioning information of each micro-seismic point in the target area and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point;
the earthquake point searching module is used for searching the initial earthquake point space model in a traversing way around by taking the shaft as a monitoring center, taking a plurality of micro earthquake points with continuous distribution trend as effective earthquake points and taking the rest micro earthquake points as ineffective earthquake points;
the path optimizing module is used for optimizing the distribution paths of all the effective vibration points in the initial vibration point space model through a path optimizing function to obtain a reinjection water distribution path of the reinjection well;
and the path marking module is used for visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
In a third aspect, a computer terminal is provided, comprising a memory, a processor and a computer program stored in the memory and executable on the processor, the processor implementing a reinjection well reinjection water distribution path visualization analysis method according to any one of the first aspects when executing the program.
Compared with the prior art, the invention has the following beneficial effects:
1. according to the reinjection well reinjection water distribution path visual analysis method provided by the invention, partial microseismic points caused by reinjection water of the non-reinjection well are filtered in consideration of the distribution distance condition and the distribution direction condition among the microseismic points, and then a path optimization function is adopted to perform distribution path optimization on all effective earthquake points in the initial earthquake point space model, so that a reinjection well reinjection water distribution path capable of simulating real conditions in detail is obtained, and the reinjection well reinjection water distribution path is subjected to visual marking in the initial earthquake point space model, so that accurate analysis on reinjection effect of reinjection water of the reinjection well is facilitated;
2. when the micro-seismic points are searched in a traversing manner, searching and analyzing are carried out on the single micro-seismic points from multiple directions, so that effective seismic points can be comprehensively searched, the situation that large differences exist in reinjection water distribution paths of reinjection wells due to the loss of the effective seismic points is avoided, and the simulation reliability and accuracy of reinjection water distribution paths of the reinjection wells are enhanced;
3. the invention considers the inertia of reinjection water of the reinjection well as fluid in flowing, takes the scanned area of the directional offset as the offset, and takes the minimum sum of the offsets of all effective vibration points as the target to determine the reinjection water distribution path of the reinjection well, which is composed of a plurality of paths, so that the efficiency of the whole optimizing process is higher, and the error between the reinjection water distribution path of the reinjection well and the real situation is smaller.
Drawings
The accompanying drawings, which are included to provide a further understanding of embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiments of the invention. In the drawings:
FIG. 1 is a flow chart in embodiment 1 of the present invention;
FIG. 2 is a partial schematic diagram of a traversal search according to embodiment 1 of the invention;
FIG. 3 is a graph showing the analysis of the offset in example 1 of the present invention;
fig. 4 is a system block diagram in embodiment 2 of the present invention.
Detailed Description
For the purpose of making apparent the objects, technical solutions and advantages of the present invention, the present invention will be further described in detail with reference to the following examples and the accompanying drawings, wherein the exemplary embodiments of the present invention and the descriptions thereof are for illustrating the present invention only and are not to be construed as limiting the present invention.
Example 1: a reinjection well reinjection water distribution path visual analysis method, as shown in fig. 1, comprises the following steps:
s1: acquiring the space positioning information of each micro-seismic point in the target area, and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point;
s2: traversing and searching an initial seismic point space model around by taking a shaft as a monitoring center, taking a plurality of microseismic points with continuous distribution trend as effective seismic points, and taking the rest microseismic points as ineffective seismic points;
s3: carrying out distribution path optimization on all effective shock points in the initial shock point space model through a path optimization function to obtain a reinjection water distribution path of the reinjection well;
s4: and visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
In this embodiment, the spatial positioning information of the microseismic points is collected by using a microseism monitoring technology, and the specific process is as follows: a well bore is taken as a monitoring center to determine a target area, and a plurality of monitoring substations distributed in an array are arranged in the target area; collecting the time when a plurality of monitoring substations receive the micro-seismic waves sent by the micro-seismic points; and establishing an equation set according to the time difference of the micro seismic waves sent by the micro seismic points received by different monitoring substations, and solving the three-dimensional space coordinates of the micro seismic points by an inversion positioning method. For example, 24 monitoring substations are used for monitoring, and the target area can be a circular area or a square area, and the monitoring substations are arranged on the ground surface. The micro-seismic waves may employ P-waves.
Because different geological structures can cause certain difference of the wave velocities of the microwave waves, particularly a section area, in order to ensure that the positioning of the micro-seismic points is more accurate, the difference of the wave velocities of the micro-seismic waves corresponding to stratum mediums at different points in a target area is not more than a wave velocity threshold; the method comprises the steps of solving the micro-seismic wave velocity applied in the three-dimensional space coordinates of the micro-seismic points by an inversion positioning method, wherein the micro-seismic wave velocity is the average value of the micro-seismic wave velocity corresponding to stratum media of each point in a target area.
In this embodiment, as shown in fig. 2, the process of traversing and searching the initial seismic point space model around by using the wellbore as a monitoring center specifically includes: determining a search radius and a search bias angle of the traversal search; selecting a micro-seismic point with a distance smaller than the searching radius and a searching path between the micro-seismic point B and the micro-seismic point B to obtain a seismic point set; if the distance between the micro-seismic point C and the micro-seismic point B is smaller than the searching radius, and one micro-seismic point A is arbitrarily selected from the seismic point set, the following conditions are satisfied: if the space included angle between the first vector of the micro-seismic point A pointing to the micro-seismic point B and the second vector of the micro-seismic point B pointing to the micro-seismic point C does not exceed the search offset angle, the search radius exists between the micro-seismic point C and the micro-seismic point B, and the micro-seismic point C is taken as an effective seismic point; and the micro-seismic points with the distance from the monitoring center smaller than the searching radius are taken as effective seismic points.
The traversal search is not limited to searching from a region with a smaller radius to a region with a larger radius, but may be searching from a region with a larger radius to a region with a smaller radius, or may be searching laterally, and is not limited thereto.
The search bias angle is not more than 90 °. And the searching radius is not smaller than the plane distribution interval of the microseismic points in the target area.
In this embodiment, the planar distribution pitch is specifically: calculating to obtain a unit area according to the ratio of the area of the horizontal plane in the target area to the number of all acquired microseismic points; and establishing an equivalent circle with the area equal to the unit area, and taking the diameter of the equivalent circle as the plane distribution interval.
When the micro-seismic points are searched in a traversing manner, the single micro-seismic points are searched and analyzed from multiple directions, so that the effective seismic points can be comprehensively searched, the condition that large differences exist in the reinjection water distribution paths of the reinjection wells due to the loss of the effective seismic points is avoided, and the simulation reliability and accuracy of the reinjection water distribution paths of the reinjection wells are enhanced.
In this embodiment, the expression of the path optimization function is specifically:
wherein,representing the jth effective shock point in the kth path; />Representing the jth effective shock point in the (k+1) th path; e represents a set of all effective shock points in the initial shock point space model; d (D) i Representing the ith effective seismic point randomly selected from the set E of effective seismic points, and the effective seismic point D i And effective shock point->A search path exists between the two;representing the effective vibration point->Relative to the effective shock point D i And (D) represents the effective shock point D i As effective shock point->Selected as the origin of the flow direction when calculating the offsetVibration point, D i-1 Representing the effective shock point D i The effective vibration points which are selected as the flow direction starting points when the offset is calculated; n (N) k Representing the number of effective shock points in the kth path; m represents the number of paths in the distribution path; and n represents an intersection symbol.
As shown in fig. 3, the calculation formula of the offset is:
β＝(x i -x i-1 ，y i -y i-1 ，z i -z i-1 )
wherein,representing the effective vibration point->Is a three-dimensional space coordinate of (2); (x) i ，y i ，z i ) Representing the effective shock point D i Is a three-dimensional space coordinate of (2); (x) i-1 ，y i-1 ，z i-1 ) Representing the effective shock point D i-1 Is a three-dimensional space coordinate of (2); alpha represents the effective shock point->And the effective vibration point D i A space vector between; beta represents the effective shock point D i And the effective vibration point D i-1 Spatial vector between them.
Example 2: the reinjection well reinjection water distribution path visual analysis system is used for realizing the reinjection well reinjection water distribution path visual analysis method described in the embodiment 1, and comprises a monitoring and positioning module, a vibration point searching module, a path optimizing module and a path labeling module as shown in fig. 4.
The monitoring and positioning module is used for acquiring the space positioning information of each micro-seismic point in the target area and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point; the earthquake point searching module is used for searching the initial earthquake point space model in a traversing way around by taking the shaft as a monitoring center, taking a plurality of micro earthquake points with continuous distribution trend as effective earthquake points and taking the rest micro earthquake points as ineffective earthquake points; the path optimizing module is used for optimizing the distribution paths of all the effective vibration points in the initial vibration point space model through a path optimizing function to obtain a reinjection water distribution path of the reinjection well; and the path marking module is used for visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
Working principle: according to the method, the distribution distance condition and the distribution direction condition among the micro-seismic points are considered, the micro-seismic points caused by the reinjection water of part of non-reinjection wells are filtered, and then a path optimization function is adopted to conduct distribution path optimization on all the effective seismic points in the initial seismic point space model, so that a reinjection well reinjection water distribution path capable of simulating real conditions in detail is obtained, and the reinjection well reinjection water distribution path is subjected to visual marking in the initial seismic point space model, so that accurate analysis on reinjection effect of reinjection water of the reinjection wells is facilitated.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The foregoing detailed description of the invention has been presented for purposes of illustration and description, and it should be understood that the invention is not limited to the particular embodiments disclosed, but is intended to cover all modifications, equivalents, alternatives, and improvements within the spirit and principles of the invention.
Claims (10)
1. The reinjection well reinjection water distribution path visual analysis method is characterized by comprising the following steps of:
acquiring the space positioning information of each micro-seismic point in the target area, and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point;
traversing and searching an initial seismic point space model around by taking a shaft as a monitoring center, taking a plurality of microseismic points with continuous distribution trend as effective seismic points, and taking the rest microseismic points as ineffective seismic points;
carrying out distribution path optimization on all effective shock points in the initial shock point space model through a path optimization function to obtain a reinjection water distribution path of the reinjection well;
and visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
2. The visual analysis method of reinjection water distribution paths of reinjection wells according to claim 1, wherein the spatial positioning information of the microseismic points is acquired by utilizing a microseismic monitoring technology, and the specific process is as follows:
a well bore is taken as a monitoring center to determine a target area, and a plurality of monitoring substations distributed in an array are arranged in the target area;
collecting the time when a plurality of monitoring substations receive the micro-seismic waves sent by the micro-seismic points;
and establishing an equation set according to the time difference of the micro seismic waves sent by the micro seismic points received by different monitoring substations, and solving the three-dimensional space coordinates of the micro seismic points by an inversion positioning method.
3. The visual analysis method of reinjection well reinjection water distribution paths according to claim 2, wherein the difference of microseismic wave velocities corresponding to stratum mediums at different points in the target area is not more than a wave velocity threshold;
the method comprises the steps of solving the micro-seismic wave velocity applied in the three-dimensional space coordinates of the micro-seismic points by an inversion positioning method, wherein the micro-seismic wave velocity is the average value of the micro-seismic wave velocity corresponding to stratum media of each point in a target area.
4. The visual analysis method of reinjection water distribution paths of reinjection wells according to claim 1, wherein the process of searching the initial seismic point space model by traversing the well bore around as a monitoring center is specifically as follows:
determining a search radius and a search bias angle of the traversal search;
selecting a micro-seismic point with a distance smaller than the searching radius and a searching path between the micro-seismic point B and the micro-seismic point B to obtain a seismic point set;
if the distance between the micro-seismic point C and the micro-seismic point B is smaller than the searching radius, and one micro-seismic point A is arbitrarily selected from the seismic point set, the following conditions are satisfied: if the space included angle between the first vector of the micro-seismic point A pointing to the micro-seismic point B and the second vector of the micro-seismic point B pointing to the micro-seismic point C does not exceed the search offset angle, the search radius exists between the micro-seismic point C and the micro-seismic point B, and the micro-seismic point C is taken as an effective seismic point;
and the micro-seismic points with the distance from the monitoring center smaller than the searching radius are taken as effective seismic points.
5. A method of visual analysis of a reinjection well reinjection water distribution path according to claim 4, characterized in that the search bias angle is not more than 90 °.
6. The visual analysis method of reinjection well reinjection water distribution paths according to claim 4, wherein the search radius is not smaller than the plane distribution interval of microseismic points in a target area;
the plane distribution interval is specifically as follows:
calculating to obtain a unit area according to the ratio of the area of the horizontal plane in the target area to the number of all acquired microseismic points;
and establishing an equivalent circle with the area equal to the unit area, and taking the diameter of the equivalent circle as the plane distribution interval.
7. The reinjection well reinjection water distribution path visualization analysis method according to claim 1, wherein the expression of the path optimization function is specifically:
wherein,representing the jth effective shock point in the kth path; />Representing the jth effective shock point in the (k+1) th path; e represents a set of all effective shock points in the initial shock point space model; d (D) i Representing the ith effective seismic point randomly selected from the set E of effective seismic points, and the effective seismic point D i And effective shock point->A search path exists between the two; />Representing the effective vibration point->Relative to the effective shock point D i And (D) represents the effective shock point D i As effective shock point->The effective vibration point D is selected as the starting point of the flow direction when the offset is calculated i-1 Representing the effective shock point D i The effective vibration points which are selected as the flow direction starting points when the offset is calculated; n (N) k Representing the number of effective shock points in the kth path; m represents the number of paths in the distribution path; and n represents an intersection symbol.
8. The visual analysis method of the reinjection water distribution path of the reinjection well according to claim 7, wherein the calculation formula of the offset is specifically as follows:
β＝(x i -x i-1 ，y i -y i-1 ，z i -z i-1 )
wherein,representing the effective vibration point->Is a three-dimensional space coordinate of (2); (x) i ，y i ，z i ) Representing the effective shock point D i Is a three-dimensional space coordinate of (2); (x) i-1 ，y i-1 ，z i-1 ) Representing the effective shock point D i-1 Is a three-dimensional space coordinate of (2); alpha represents the effective shock point->And the effective vibration point D i A space vector between; beta represents the effective shock point D i And the effective vibration point D i-1 Spatial vector between them.
9. A reinjection well reinjection water distribution path visualization analysis system, comprising:
the monitoring and positioning module is used for acquiring the space positioning information of each micro-seismic point in the target area and establishing an initial seismic point space model according to the space positioning information of each micro-seismic point;
the earthquake point searching module is used for searching the initial earthquake point space model in a traversing way around by taking the shaft as a monitoring center, taking a plurality of micro earthquake points with continuous distribution trend as effective earthquake points and taking the rest micro earthquake points as ineffective earthquake points;
the path optimizing module is used for optimizing the distribution paths of all the effective vibration points in the initial vibration point space model through a path optimizing function to obtain a reinjection water distribution path of the reinjection well;
and the path marking module is used for visually marking the reinjection water distribution path of the reinjection well in the initial seismic point space model to obtain a seismic point visual distribution model.
10. A computer terminal comprising a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor implements a method for visual analysis of a reinjection well reinjection water distribution path according to any one of claims 1 to 8 when the program is executed by the processor.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311616181.6A CN117637064A (en) | 2023-11-29 | 2023-11-29 | Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311616181.6A CN117637064A (en) | 2023-11-29 | 2023-11-29 | Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN117637064A true CN117637064A (en) | 2024-03-01 |
Family
ID=90019524
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311616181.6A Pending CN117637064A (en) | 2023-11-29 | 2023-11-29 | Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117637064A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118884520A (en) * | 2024-07-18 | 2024-11-01 | 中国地震局地球物理研究所 | Method and system for identifying structural surfaces and flow paths induced by fluid injection |
-
2023
- 2023-11-29 CN CN202311616181.6A patent/CN117637064A/en active Pending
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118884520A (en) * | 2024-07-18 | 2024-11-01 | 中国地震局地球物理研究所 | Method and system for identifying structural surfaces and flow paths induced by fluid injection |
| CN118884520B (en) * | 2024-07-18 | 2025-02-11 | 中国地震局地球物理研究所 | Method and system for identifying structural surfaces and flow paths induced by fluid injection |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108254784B (en) | Fault identification method, device and system based on two-dimensional seismic data | |
| US10345463B2 (en) | Methods and systems for using known source events in seismic data processing | |
| CN109828236A (en) | A kind of microseism/acoustic emission source locating method in labyrinth containing dead zone | |
| CN112305591B (en) | Tunnel advance geological prediction method, computer readable storage medium | |
| CN111221036B (en) | A method and system for source location of target area containing unknown voids | |
| CN103105624A (en) | Longitudinal and transversal wave time difference positioning method based on base data technology | |
| Chamarczuk et al. | Unsupervised learning used in automatic detection and classification of ambient‐noise recordings from a large‐N array | |
| CN117637064A (en) | Reinjection well reinjection water distribution path visual analysis method, reinjection well reinjection water distribution path visual analysis system and reinjection well reinjection water distribution path visual analysis terminal | |
| CN103105622B (en) | Based on the homotype ripple time difference positioning method of database technology | |
| CN108845358B (en) | Tomography and the recognition methods of structural anomaly body and device | |
| US20240230942A9 (en) | Method and System for Real-Time Calculating a Microseismic Focal Mechanism Based on Deep Learning | |
| CN119418066A (en) | A fracture contour extraction method for well logging images | |
| CN113568042B (en) | Streamer acquisition method and system for improving repeatability of time-lapse seismic data | |
| CN109709602B (en) | Remote detection sound wave offset imaging method, device and system | |
| CN107664771B (en) | A kind of microseism Full wave shape localization method based on likeness coefficient | |
| US7206444B2 (en) | System and method for locating multiple peak summits in three-dimensional data | |
| Alghalandis et al. | Similarity analysis of discrete fracture networks | |
| CN110687598A (en) | A method and device for accelerating microseismic numerical simulation | |
| CN105911138A (en) | Impact imaging method and system based on average response energy | |
| CN113568041B (en) | Repeatability analysis method and system for time-lapse seismic three-dimensional towing cable acquired data | |
| CN116660980A (en) | A Microseismic Location Method Based on Improved Equation | |
| CN104111476A (en) | Method and device for building formation velocity field | |
| CN104597504A (en) | High-efficiency microseism seismic focus locating method | |
| CN111189926B (en) | Method and system for identifying structure hole position based on global search | |
| CN116184500B (en) | Real-time inversion method and device for ground stress of tunnel based on microseismic information |
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 |