CN117451857A - Shale gas reservoir space detection method and related equipment thereof - Google Patents
Shale gas reservoir space detection method and related equipment thereof Download PDFInfo
- Publication number
- CN117451857A CN117451857A CN202311770787.5A CN202311770787A CN117451857A CN 117451857 A CN117451857 A CN 117451857A CN 202311770787 A CN202311770787 A CN 202311770787A CN 117451857 A CN117451857 A CN 117451857A
- Authority
- CN
- China
- Prior art keywords
- dimensional model
- target
- shale
- data
- acquisition matrix
- 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
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 151
- 239000011159 matrix material Substances 0.000 claims abstract description 151
- 239000000523 sample Substances 0.000 claims abstract description 58
- 238000012216 screening Methods 0.000 claims abstract description 35
- 238000012545 processing Methods 0.000 claims abstract description 25
- 238000000034 method Methods 0.000 claims description 28
- 230000008569 process Effects 0.000 claims description 20
- 238000004891 communication Methods 0.000 claims description 15
- 238000004590 computer program Methods 0.000 claims description 8
- 230000001174 ascending effect Effects 0.000 claims description 3
- 238000012805 post-processing Methods 0.000 abstract description 5
- 230000000875 corresponding effect Effects 0.000 description 99
- 239000007789 gas Substances 0.000 description 74
- 230000005540 biological transmission Effects 0.000 description 8
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 8
- 230000001276 controlling effect Effects 0.000 description 7
- 230000006870 function Effects 0.000 description 5
- 239000003345 natural gas Substances 0.000 description 4
- 230000009471 action Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 3
- 239000003086 colorant Substances 0.000 description 3
- 230000008878 coupling Effects 0.000 description 3
- 238000010168 coupling process Methods 0.000 description 3
- 238000005859 coupling reaction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 238000012217 deletion Methods 0.000 description 2
- 230000037430 deletion Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 239000011435 rock Substances 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000012800 visualization Methods 0.000 description 2
- 235000012431 wafers Nutrition 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 230000002596 correlated effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4418—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a model, e.g. best-fit, regression analysis
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/02—Analysing fluids
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/02—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems using reflection of acoustic waves
- G01S15/06—Systems determining the position data of a target
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S15/00—Systems using the reflection or reradiation of acoustic waves, e.g. sonar systems
- G01S15/88—Sonar systems specially adapted for specific applications
- G01S15/89—Sonar systems specially adapted for specific applications for mapping or imaging
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/52—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00
- G01S7/539—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S15/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T11/00—2D [Two Dimensional] image generation
- G06T11/40—Filling a planar surface by adding surface attributes, e.g. colour or texture
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T17/00—Three dimensional [3D] modelling, e.g. data description of 3D objects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T3/00—Geometric image transformations in the plane of the image
- G06T3/40—Scaling of whole images or parts thereof, e.g. expanding or contracting
-
- 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)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Life Sciences & Earth Sciences (AREA)
- Theoretical Computer Science (AREA)
- Biochemistry (AREA)
- General Health & Medical Sciences (AREA)
- Analytical Chemistry (AREA)
- Immunology (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Computer Networks & Wireless Communication (AREA)
- Chemical & Material Sciences (AREA)
- Health & Medical Sciences (AREA)
- Geology (AREA)
- Computer Graphics (AREA)
- Signal Processing (AREA)
- Mining & Mineral Resources (AREA)
- Software Systems (AREA)
- Geometry (AREA)
- Geochemistry & Mineralogy (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Environmental & Geological Engineering (AREA)
- Architecture (AREA)
- Computer Hardware Design (AREA)
- General Engineering & Computer Science (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
The application belongs to the technical field of shale gas reservoir space detection, and discloses a shale gas reservoir space detection method and related equipment, wherein a phased array ultrasonic probe is used for detection, only one detection hole is needed to be arranged, and the phased array ultrasonic probe is enabled to transversely scan while moving along the detection hole, so that the distribution information of the shale gas reservoir space in the target shale in a three-dimensional space can be obtained, and the detection efficiency is high; calling a pre-established reference two-dimensional model and a pre-established reference three-dimensional model, and performing color filling on the called models according to the distribution condition of the shale gas reservoir space so as to obtain a shale gas reservoir space distribution model, so that on one hand, the detection result can be displayed more intuitively, on the other hand, the generation efficiency of the shale gas reservoir space distribution model can be improved, and the detection efficiency is further improved; the data processing amount of the post-processing of the data by the upper computer can be reduced by screening and zeroing the data of the acquisition matrix, so that the detection efficiency is further improved.
Description
Technical Field
The application relates to the technical field of shale gas reservoir space detection, in particular to a shale gas reservoir space detection method and related equipment thereof.
Background
Shale gas is an energy source stored inside shale, and belongs to an unconventional energy source. Because shale gas is natural gas existing in the middle of an underground rock layer, the three-dimensional shape of a reservoir space has large variability, and therefore, the shale gas has certain difficulty and complexity compared with the traditional natural gas energy exploitation. In the existing shale gas reservoir space detection method, CT fault scanning is firstly carried out on the ground, then a corresponding drilling tower is built after the approximate direction is determined, then vertical drilling is carried out through a drill bit, when the drill bit meets rock, the drill bit is rotated by 90 degrees and then horizontal drilling is carried out, after the drill bit is taken out, a detection device is penetrated into the ground to detect the three-dimensional storage space of shale, and then the volume of the shale gas reservoir space is approximately determined, but the accuracy of the determination method has larger error due to the self-variability of shale, and the real content of shale gas also has larger error due to the shape variability of shale.
On the other hand, the existing shale gas reservoir space detection method has two defects, namely, the result is not visual enough, and the representation of the local position of the natural gas storage in the shale is not accurate and obvious enough; and secondly, the overall speed of detection is slower due to the larger data volume. Therefore, how to realize the rapidity and visualization of the detection of the inside of the shale gas reservoir space is particularly important.
Disclosure of Invention
The purpose of the application is to provide a shale gas reservoir space detection method and related equipment thereof, which can improve the detection efficiency of shale gas reservoir space detection and the visualization degree of detection results.
In a first aspect, the application provides a shale gas reservoir space detection system, which comprises a collection device and an upper computer, wherein the collection device comprises a phased array ultrasonic probe; the acquisition device is in communication connection with the upper computer, and the upper computer stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the acquisition device is used for controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections, wherein the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data, an acquisition matrix of each scanning cross section is generated according to the reflected signal amplitude data of each scanning cross section, screening and zeroing processing is carried out on the data of each acquisition matrix, and the acquisition matrix and the corresponding reflected position data after the screening and zeroing processing are sent to the upper computer;
The upper computer is used for determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to the received acquisition matrix and the corresponding reflection position data, respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model, carrying out color filling on the target two-dimensional model corresponding to each scanning cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data to obtain corresponding two-dimensional color images, stacking and filling the two-dimensional color images into the target three-dimensional model, and generating a shale gas storage space distribution model of the target shale.
The phased array ultrasonic probe can scan and detect the whole scanning cross section of the target shale at one position, and only one detection hole is needed to be arranged when the shale gas reservoir space of the target shale is detected, and the phased array ultrasonic probe can transversely scan while moving along the detection hole, so that the distribution information of the shale gas reservoir space in the target shale in a three-dimensional space can be obtained, and the detection efficiency is high; calling a pre-established reference two-dimensional model and a pre-established reference three-dimensional model, and performing color filling on the called models according to the distribution condition of the shale gas reservoir space so as to obtain a shale gas reservoir space distribution model, so that on one hand, the detection result can be displayed more intuitively, on the other hand, the generation efficiency of the shale gas reservoir space distribution model can be improved, and the detection efficiency is further improved; the data processing amount of the post-processing of the data by the upper computer can be reduced by screening and zeroing the data of the acquisition matrix, so that the detection efficiency is further improved.
Preferably, the ultrasonic detection data of each scanning cross section comprises one-dimensional reflected signal amplitude data obtained by multi-frame detection along different detection angles and reflected position data corresponding to each amplitude data in the reflected signal amplitude data;
the acquisition device performs the following steps when generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section:
taking the frame number of the reflected signal amplitude data of the same scanning cross section as the column number of the corresponding acquisition matrix, and determining the row number of the corresponding acquisition matrix according to the reflected signal amplitude data of the scanning cross section;
generating the acquisition matrix according to the number of rows and the number of columns and initializing the acquisition matrix to be a zero matrix;
and filling the amplitude data of the reflected signals of each frame into the data of each column of the acquisition matrix in sequence.
Preferably, the acquisition device performs screening and zeroing processing on the data of each acquisition matrix:
ascending order sorting is carried out on all the numerical values in the same acquisition matrix;
selecting a numerical value of a preset ranking in the acquisition matrix as a screening threshold of the acquisition matrix according to the ranking result;
And setting the numerical value smaller than the screening threshold value in the acquisition matrix to zero.
By the method, the data volume can be reduced while effective characteristic data is reserved, so that the data transmission bandwidth requirement is reduced, the shale gas reservoir space distribution model is built quickly, and the detection result is ensured to have enough accuracy.
Preferably, the acquisition device performs when transmitting the acquisition matrix after screening and zeroing processing and the corresponding reflection position data to the upper computer:
and transmitting the acquisition matrix after screening and zero setting processing and reflection position data corresponding to each non-zero value in the acquisition matrix to the upper computer.
Therefore, the data volume can be further reduced, the burden of data transmission is reduced, and the data transmission bandwidth requirement is reduced.
Preferably, the upper computer determines a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to the received acquisition matrix and the corresponding reflection position data, and respectively marks the reference two-dimensional model and the target three-dimensional model as a target two-dimensional model and a target three-dimensional model, and then performs the following steps:
acquiring edge contour lines of the scanning cross sections according to the received acquisition matrix and the corresponding reflection position data;
Determining a reference two-dimensional model matched with each scanning cross section from the reference two-dimensional models stored by the upper computer according to the edge contour line of each scanning cross section, and marking the reference two-dimensional model as a target two-dimensional model;
scaling the target two-dimensional model corresponding to each scanning cross section according to the edge contour line of each scanning cross section;
determining a reference three-dimensional model matched with the target shale from the reference three-dimensional models stored by the upper computer according to the target two-dimensional model subjected to scaling processing corresponding to each scanning cross section, and marking the reference three-dimensional model as a target three-dimensional model;
and scaling the target three-dimensional model according to the scaled target two-dimensional model corresponding to each scanning cross section.
Preferably, the upper computer performs color filling on the target two-dimensional model corresponding to each scanning cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data to obtain a corresponding two-dimensional color map, and performs:
determining the color corresponding to each non-zero value according to the non-zero value of the acquisition matrix, and marking the color as a filling color;
and setting the color of the corresponding position point in the target two-dimensional model as the corresponding filling color according to the reflection position data corresponding to the non-zero numerical value of the acquisition matrix.
In a second aspect, the application provides a shale gas reservoir space detection method, which is applied to an acquisition device, wherein the acquisition device comprises a phased array ultrasonic probe and is in communication connection with an upper computer;
the shale gas reservoir space detection method comprises the following steps:
A1. controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections; the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data;
A2. generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section;
A3. screening and zeroing the data of each acquisition matrix;
A4. and sending the acquisition matrix after screening and zeroing treatment and corresponding reflection position data to the upper computer, so that the upper computer generates a shale gas reservoir space distribution model of the target shale based on a reference two-dimensional model and a reference three-dimensional model of various typical structures of the shale, which are established in advance, according to the acquisition matrix and the reflection position data.
In a third aspect, the application provides a shale gas reservoir space detection method, which is applied to an upper computer, wherein the upper computer is in communication connection with an acquisition device, and the acquisition device comprises a phased array ultrasonic probe; the upper computer stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the shale gas reservoir space detection method comprises the following steps:
B1. determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to an acquisition matrix and corresponding reflection position data from the acquisition device, and respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model; the acquisition matrix and the corresponding reflection position data are generated by the acquisition device according to ultrasonic detection data of a scanning cross section obtained by transversely scanning the target shale in the process that the phased array ultrasonic probe moves in a detection hole extending along the surface of the target shale;
B2. according to the numerical value of the acquisition matrix and the corresponding reflection position data, performing color filling on the target two-dimensional model corresponding to each scanning cross section to obtain a corresponding two-dimensional color map;
B3. and filling the two-dimensional color map stack into the target three-dimensional model to generate a shale gas reservoir space distribution model of the target shale.
In a fourth aspect, the present application provides an electronic device comprising a processor and a memory, the memory storing a computer program executable by the processor, which when executed, performs steps in a shale gas reservoir space detection method as hereinbefore described.
In a fifth aspect, the present application provides a computer readable storage medium having stored thereon a computer program which, when executed by a processor, performs steps in a shale gas reservoir space detection method as hereinbefore described.
The beneficial effects are that: according to the shale gas reservoir space detection method and the related equipment, as the phased array ultrasonic probe can scan and detect the whole scanning cross section of the target shale at one position, when the shale gas reservoir space of the target shale is detected, only one detection hole is needed to be arranged, and the phased array ultrasonic probe can transversely scan while moving along the detection hole, so that the distribution information of the shale gas reservoir space in the target shale in a three-dimensional space can be obtained, and the detection efficiency is high; calling a pre-established reference two-dimensional model and a pre-established reference three-dimensional model, and performing color filling on the called models according to the distribution condition of the shale gas reservoir space so as to obtain a shale gas reservoir space distribution model, so that on one hand, the detection result can be displayed more intuitively, on the other hand, the generation efficiency of the shale gas reservoir space distribution model can be improved, and the detection efficiency is further improved; the data processing amount of the post-processing of the data by the upper computer can be reduced by screening and zeroing the data of the acquisition matrix, so that the detection efficiency is further improved.
Drawings
Fig. 1 is a schematic structural diagram of a shale gas reservoir space detection system according to an embodiment of the present application.
Fig. 2 is a flowchart of a shale gas reservoir space detection method according to an embodiment of the present application.
Fig. 3 is a flowchart of another shale gas reservoir space detection method provided in an embodiment of the present application.
Fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Fig. 5 is a front view of an exemplary target shale.
Fig. 6 is a schematic scanning cross-section.
Description of the reference numerals: 1. a collection device; 101. a phased array ultrasonic probe; 2. an upper computer; 301. a processor; 302. a memory; 303. a communication bus; 90. target shale; 91. a detection hole; 92. and (5) erecting holes.
Detailed Description
The following description of the embodiments of the present application will be made clearly and completely with reference to the drawings in the embodiments of the present application, and it is apparent that the described embodiments are only some embodiments of the present application, not all embodiments. The components of the embodiments of the present application, which are generally described and illustrated in the figures herein, may be arranged and designed in a wide variety of different configurations. Thus, the following detailed description of the embodiments of the present application, as provided in the accompanying drawings, is not intended to limit the scope of the application, as claimed, but is merely representative of selected embodiments of the application. All other embodiments, which can be made by those skilled in the art based on the embodiments of the present application without making any inventive effort, are intended to be within the scope of the present application.
It should be noted that: like reference numerals and letters denote like items in the following figures, and thus once an item is defined in one figure, no further definition or explanation thereof is necessary in the following figures. Meanwhile, in the description of the present application, the terms "first", "second", and the like are used only to distinguish the description, and are not to be construed as indicating or implying relative importance.
Referring to fig. 1, fig. 1 is a shale gas reservoir space detection system in some embodiments of the present application, including a collecting device 1 and an upper computer 2, where the collecting device 1 includes a phased array ultrasonic probe 101; the acquisition device 1 is in communication connection with the upper computer 2, and the upper computer 2 stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the acquisition device 1 is used for controlling the phased array ultrasonic probe 101 to move in a detection hole 91 extending along the surface of the target shale 90, and transversely scanning the target shale 90 in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections, wherein the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data, generating acquisition matrixes of each scanning cross section according to the reflected signal amplitude data of each scanning cross section, carrying out screening and zeroing treatment on the data of each acquisition matrix, and transmitting the acquisition matrixes after the screening and zeroing treatment and the corresponding reflected position data to the upper computer 2;
The upper computer 2 is configured to determine a reference two-dimensional model matched with each scan cross section (corresponding reference two-dimensional model is determined for each scan cross section) and a reference three-dimensional model matched with the target shale according to the received acquisition matrix and corresponding reflection position data, record the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model respectively, perform color filling on the target two-dimensional model corresponding to each scan cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data, obtain a corresponding two-dimensional color map, stack and fill the two-dimensional color map into the target three-dimensional model, and generate a shale gas reservoir space distribution model of the target shale.
Since the phased array ultrasonic probe 101 can scan and detect the whole scanning cross section of the target shale 90 at one position, when the shale gas reservoir space of the target shale 90 is detected, only one detection hole 91 is required to be arranged, and the phased array ultrasonic probe 101 can transversely scan while moving along the detection hole 91, so that the distribution information of the shale gas reservoir space in the target shale 90 in a three-dimensional space (for example, a in fig. 5 and 6 is the shale gas reservoir space) can be obtained, and the detection efficiency is high; calling a pre-established reference two-dimensional model and a pre-established reference three-dimensional model, and performing color filling on the called models according to the distribution condition of the shale gas reservoir space so as to obtain a shale gas reservoir space distribution model, so that on one hand, the detection result can be displayed more intuitively, on the other hand, the generation efficiency of the shale gas reservoir space distribution model can be improved, and the detection efficiency is further improved; the data processing amount of the post-processing of the data by the upper computer 2 can be reduced by screening and zeroing the data of the acquisition matrix, so that the detection efficiency is further improved.
Preferably, the ultrasonic detection data of each scanning cross section comprises one-dimensional reflected signal amplitude data obtained by multi-frame detection along different detection angles and reflection position data corresponding to each amplitude data in each frame of reflected signal amplitude data.
The phased array ultrasonic probe 101 includes a plurality of piezoelectric vibration wafers arranged in an array, and by controlling a vibration phase difference of each piezoelectric vibration wafer, adjustment of an ultrasonic beam direction emitted by the phased array ultrasonic probe 101 (this is the prior art) can be achieved, so that adjustment of a detection angle is achieved, and the phased array ultrasonic probe 101 can achieve scanning detection of the whole scanning cross section of the target shale 90 by adjusting the ultrasonic beam direction at the same position. After the phased array ultrasonic probe 101 emits an ultrasonic wave beam in one direction, when the ultrasonic wave beam encounters an interface of two different mediums (for example, the boundary of shale gas reservoir space inside the target shale 90 and the boundary of the target shale 90), part of ultrasonic waves are reflected to form a reflected signal, the amplitude of the reflected signal is positively correlated with the size of the interface of the different mediums encountered by the ultrasonic wave beam, the reflected signal is received by the phased array ultrasonic probe 101, and the amplitude of the reflected signal received by the phased array ultrasonic probe 101 at the same detection angle is arranged according to the receiving time to obtain a frame of reflected signal amplitude data with time information. When the same scanning cross section is scanned and detected, multi-frame reflected signal amplitude data are obtained.
For example, in fig. 6, the center lines b represent the center lines of the ultrasonic beam, each center line b represents the detection at a detection angle, the angle between the adjacent center lines b is a swing angle, and the smaller the swing angle is, the more the number of frames of the reflected signal amplitude data obtained at the same scanning cross section is, the finer the scanning is, the more accurate the detection result is, but the larger the data processing amount is, and the swing angle can be specifically adjusted according to the actual needs.
At the same detection angle, the emission time of the ultrasonic wave beam and the receiving time of the reflected signals are recorded, and the detection angle is also recorded, so that the positions of the reflection points of the reflected signals can be calculated according to the following formula to obtain reflection position data:
q=（r，θ）；
r=（t1-t0）*v/2；
where q is the position of the reflection point expressed in polar coordinates, r is the distance between the reflection point and the phased array ultrasonic probe 101, θ is the detection angle, t1 is the reception time of the reflected signal, t0 is the transmission time of the ultrasonic beam, and v is the propagation speed of the ultrasonic wave (which is a known value).
In practice, the position of the reflection point expressed in polar coordinates can be converted into coordinates in the ultrasonic probe coordinate system, and the conversion formula is as follows:
x=r*sinθ；
y=r*cosθ；
wherein x is the abscissa of the position of the reflection point under the ultrasonic probe coordinate system, y is the ordinate of the position of the reflection point under the ultrasonic probe coordinate system, and the detection angle θ is the included angle between the central line of the ultrasonic beam and the central axis of the phased array ultrasonic probe 101 (the central axis is perpendicular to the array plane of the phased array ultrasonic probe 101 and passes through the center of the array plane).
Referring to fig. 5, during the detection, a hole needs to be drilled downwards from the ground to obtain a vertical hole 92 extending to the upper surface of the target shale 90, then a drill bit is turned to drill a detection hole 91 along the upper surface of the target shale 90 (this is the prior art), after the drill bit is taken out, the phased array ultrasonic probe 101 is sent into the detection hole 91 and moves along the detection hole 91, and during the movement, the scanning detection of a plurality of scanning cross sections is realized by adjusting the detection angle; the phased array ultrasonic probe 101 can move in a stepping manner, and scanning detection is performed once every moving step; because the speed of scanning detection is high, the phased array ultrasonic probe 101 can also move at a constant speed and perform scanning detection in the moving process, the moving distance of the phased array ultrasonic probe 101 in the scanning detection process is small, the scanning detection process can be considered to be scanning detection on the same scanning cross section, and the next scanning detection can be performed at intervals of preset time (which can be set according to actual needs) after each scanning detection. The denser the scanned cross sections detected by scanning are, the more accurate the finally obtained shale gas reservoir space distribution model is, but the larger the data processing capacity is, the number of the scanned cross sections can be set according to actual needs. For example, in fig. 5, each line c represents the position of one scan cross section.
Specifically, the acquisition device 1 performs, when generating an acquisition matrix for each scan cross section from the reflected signal amplitude data for each scan cross section:
taking the frame number of the reflected signal amplitude data of the same scanning cross section as the column number of the corresponding acquisition matrix, and determining the row number of the corresponding acquisition matrix according to the reflected signal amplitude data of the scanning cross section;
generating an acquisition matrix according to the number of rows and the number of columns and initializing the acquisition matrix into a zero matrix;
and sequentially filling the amplitude data of the reflected signals of each frame into the data of each column of the acquisition matrix.
For example, assuming that a scan cross section performs detection at n detection angles during scan detection, n frames of reflected signal amplitude data are obtained, so that the number of columns of the acquisition matrix corresponding to the scan cross section is n.
When the ultrasonic beam reaches the lower surface of the target shale 90, a reflected signal with a larger amplitude is formed, the reflected signal is hereinafter referred to as a termination signal, all data located after the termination signal in the amplitude data of the reflected signal can be deleted, and then the maximum value of the amplitude data quantity of the amplitude data of each of the deleted reflected signals is taken as the row number of the corresponding acquisition matrix. The last one of the amplitude data of the reflected signals, which is greater than a first preset threshold (which can be set according to actual needs), may be used as a termination signal, so that when the deletion processing is performed, the data after the last one of the amplitude data greater than the first preset threshold is deleted.
For example, it is assumed that n frames of reflected signal amplitude data are corresponding to one scan cross section, after the deletion processing is performed, the number of amplitude data of each frame of reflected signal amplitude data is m1 and m2 … … mn, and if the maximum value is m2, the number of rows of the acquisition matrix corresponding to the scan cross section is m2.
When the amplitude data of each frame of reflected signal is sequentially filled into each column of data of the acquisition matrix, the amplitude data of each frame of reflected signal can be sequentially filled into each column of data of the acquisition matrix from left to right according to the acquisition order of the amplitude data of each frame of reflected signal. When filling a frame of reflected signal amplitude data, replacing the amplitude data of the frame of reflected signal amplitude data with the 0 value in the corresponding column data in the acquisition matrix from top to bottom. For example, the i-th frame of reflected signal amplitude data is filled into the i-th column data of the acquisition matrix, the i-th frame of reflected signal amplitude data includes j amplitude data, so that the 1-th amplitude data of the i-th frame of reflected signal amplitude data is replaced with the 0 value of the (1, i) data of the acquisition matrix, the 2-th amplitude data of the i-th frame of reflected signal amplitude data is replaced with the 0 value of the (2, i) data of the acquisition matrix, the 3-th amplitude data of the i-th frame of reflected signal amplitude data is replaced with the 0 value of the (3,i) data of the acquisition matrix, and so on until the j-th amplitude data of the i-th frame of reflected signal amplitude data is replaced with the 0 value of the (j, i) data of the acquisition matrix.
Preferably, the acquisition device 1 performs, when performing screening and zeroing processing on the data of each acquisition matrix:
ascending order sorting is carried out on each numerical value in the same acquisition matrix;
selecting a numerical value of a preset ranking in the acquisition matrix (assuming that the preset ranking is T, the numerical value of the preset ranking is the numerical value of the ranking is T in the acquisition matrix) as a screening threshold of the acquisition matrix according to the sequencing result;
and setting the value smaller than the screening threshold value in the acquisition matrix to zero.
In practice, there may be a large number of small gaps or small voids in the target shale 90, which are not effective shale gas reservoir space, but reflect the ultrasonic beam to form relatively weak reflected signals, which may be received by the phased array ultrasonic probe 101, so that there is a large number of smaller amplitude data in the amplitude data of each frame of reflected signals, which have no beneficial effect on the detection of shale gas reservoir space, but rather greatly increase the data throughput. By the method, the data volume can be reduced while effective characteristic data is reserved, so that the data transmission bandwidth requirement is reduced, the shale gas reservoir space distribution model is built quickly (the upper computer 2 only needs to perform calculation processing according to non-zero values in the acquisition matrix when performing subsequent processing), and the detection result is ensured to have enough accuracy.
The preset ranking can be N, M/4, M is the number of rows of the acquisition matrix, and N is the number of columns of the acquisition matrix; but is not limited thereto, and may be specifically set according to actual needs.
Further, the acquisition device 1 performs, when transmitting the acquisition matrix after the screening and zeroing processing and the corresponding reflection position data to the upper computer 2:
and transmitting the filtered and zeroed acquisition matrix and the reflection position data corresponding to the non-zero values in the acquisition matrix to the upper computer 2.
Namely, the reflection position data corresponding to the original non-zero numerical value and the set numerical value is abandoned, so that the data quantity can be further reduced, the burden of data transmission is reduced, and the data transmission bandwidth requirement is reduced.
In fact, in order to facilitate the upper computer 2 to determine the correspondence between the values in the acquisition matrix and the reflected position data, the reflected position data may be recorded by using a position matrix, where, since the reflected position data is two-dimensional coordinate data (polar coordinate data or coordinates of the ultrasonic probe coordinate system), the position matrix includes a first position matrix and a second position matrix, the number of rows and columns of the first position matrix and the second position matrix are the same as the number of rows and columns of the acquisition matrix and are both initialized to be a zero matrix, the first position matrix is used for recording a first coordinate value of the reflected position data, and the recording position of the first coordinate value in the first position matrix is the same as the recording position of the corresponding amplitude data in the acquisition matrix; the second position matrix is used for recording a second coordinate value of the reflection position data, and the recording position of the second coordinate value in the second position matrix is the same as the recording position of the corresponding amplitude value data in the acquisition matrix. When transmitting the reflection position data to the host computer 2, the acquisition device 1 transmits the reflection position data to the first position matrix and the second position matrix. Therefore, the upper computer 2 can extract numerical values from corresponding positions in the first position matrix and the second position matrix according to the recording positions of the amplitude data in the acquisition matrix to form reflection position data corresponding to the amplitude data.
Wherein the reference two-dimensional model and the reference three-dimensional model are models previously established according to shape sizes of various typical structures of shale. These models are stored in a reference model database of the host computer 2.
In some embodiments, the upper computer 2 performs when determining a reference two-dimensional model matching each scan cross section and a reference three-dimensional model matching the target shale 90, respectively, from the received acquisition matrix and corresponding reflection position data, as a target two-dimensional model and a target three-dimensional model:
acquiring edge contour lines of each scanning cross section according to the received acquisition matrix and the corresponding reflection position data;
determining a reference two-dimensional model matched with each scanning cross section from the reference two-dimensional models stored by the upper computer according to the edge contour line of each scanning cross section, and marking the reference two-dimensional model as a target two-dimensional model;
scaling the target two-dimensional model corresponding to each scanning cross section according to the edge contour line of each scanning cross section;
according to the target two-dimensional model after the scaling treatment corresponding to each scanning cross section, determining a reference three-dimensional model matched with the target shale from the reference three-dimensional models stored by the upper computer 2, and marking the reference three-dimensional model as a target three-dimensional model;
And scaling the target three-dimensional model according to the scaled target two-dimensional model corresponding to each scanning cross section.
The reflection position data corresponding to the last non-zero value in each column of the acquisition matrix are position data of points on the contour line of the corresponding scanning cross section, and the edge contour line of the corresponding scanning cross section can be obtained by extracting the reflection position data corresponding to the last non-zero value in each column of the acquisition matrix and performing curve fitting.
And calculating the similarity between the scanning cross section and each reference two-dimensional model according to the edge contour line of the scanning cross section, and extracting the reference two-dimensional model corresponding to the maximum value of the similarity as a target two-dimensional model of the scanning cross section.
When the scaling treatment is carried out on the target two-dimensional model, the target two-dimensional model can be scaled in equal proportion until the area of the target two-dimensional model is equal to the area of an area surrounded by the edge contour line of the corresponding scanning cross section; or respectively generating the minimum circumscribed rectangle of the target two-dimensional model and the minimum circumscribed rectangle of the edge contour line of the corresponding scanning cross section, and then enabling the target two-dimensional model to be scaled in equal proportion along with the minimum circumscribed rectangle of the target two-dimensional model until the diagonal length of the minimum circumscribed rectangle of the target two-dimensional model is equal to the diagonal length of the minimum circumscribed rectangle of the edge contour line of the corresponding scanning cross section. And scaling the target two-dimensional model to enable the target two-dimensional model to be more matched with the corresponding scanning cross section.
The position of the phased array ultrasonic probe 101 (which is a position in the world coordinate system) may be calculated (for example, calculated according to the moving speed or moving step length of the phased array ultrasonic probe 101 and the extending track of the probe hole 91), or may be measured by a positioning device (such as an inertial positioning device) disposed on the phased array ultrasonic probe 101. The position (world coordinate position) of the phased array ultrasonic probe 101 when each scanning cross section is detected can be obtained, the target two-dimensional models after each scaling treatment are arranged according to the relative position relation between the positions of the phased array ultrasonic probe 101 when each scanning cross section is detected, then a plurality of contour points are extracted from the edge contour lines of the target two-dimensional models after each scaling treatment and used for fitting the surface curved surface of the target shale 90, the similarity between the surface curved surface of the target shale 90 and each reference three-dimensional model is calculated according to the surface curved surface obtained by fitting (the similarity is the prior art and is not described in detail), and the reference three-dimensional model corresponding to the maximum value of the similarity is used as the target three-dimensional model.
When the scaling treatment is carried out on the target three-dimensional model, the target three-dimensional model can be scaled in equal proportion until the volume of the target three-dimensional model is equal to the volume of the area surrounded by the surface curved surface; or respectively generating a minimum bounding box of the target three-dimensional model and a minimum bounding box of the surface curved surface, and then enabling the target three-dimensional model to be scaled in equal proportion along with the minimum bounding box of the target three-dimensional model until the length of the longest diagonal of the minimum bounding box of the target three-dimensional model is equal to the length of the longest diagonal of the minimum bounding box of the surface curved surface. The target three-dimensional model is more closely matched to the target shale 90 by scaling the target three-dimensional model.
In some embodiments, the upper computer 2 performs color filling on the target two-dimensional model corresponding to each scanning cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data, and when obtaining the corresponding two-dimensional color map, performs:
determining the color corresponding to each non-zero value according to the non-zero value of the acquisition matrix, and marking the color as a filling color;
and setting the color of the corresponding position point in the target two-dimensional model as the corresponding filling color according to the reflection position data corresponding to the non-zero numerical value of the acquisition matrix.
The non-zero value of the acquisition matrix is amplitude data of the reflected signal, and represents the size of the reflecting structure, and generally reflects the size of shale gas reservoir space causing reflection. The chromaticity or the shade of the corresponding color can be preset for different amplitude data ranges, and then the chromaticity or the shade of the color corresponding to the non-zero value is determined according to the amplitude data range of the non-zero value of the acquisition matrix and recorded as the filling color. For example, the larger the amplitude data, the more red the color and the smaller the amplitude data, the more yellow the color; alternatively, the filling colors are red, and the larger the amplitude data is, the darker the red color is, and the smaller the amplitude data is, the lighter the red color is, but not limited thereto.
Wherein, the colors of other positions in the target two-dimensional model except the position set as the corresponding filling color are kept as the first preset color (can be set according to actual needs, for example, white, but is not limited thereto).
Further, the upper computer 2 performs, when filling the two-dimensional color map stack into the target three-dimensional model to generate a shale gas reservoir spatial distribution model of the target shale:
the position (world coordinate position) of the phased array ultrasonic probe 101 at the time of scanning and detecting each scanning cross section is acquired, and each two-dimensional color map is arranged in the target three-dimensional model according to the relative positional relationship between the positions of the phased array ultrasonic probe 101 at the time of scanning and detecting each scanning cross section.
When the two-dimensional color patterns are arranged tightly enough (for example, when the distance between the adjacent two-dimensional color patterns is smaller than a preset distance threshold, the arrangement is considered to be tight enough, the preset distance threshold can be set according to actual needs), and a shale gas reservoir space distribution model is obtained after the arrangement is completed; and when the two-dimensional color maps are not arranged tightly enough, interpolating and filling the colors of the target three-dimensional model at the space positions between the two-dimensional color maps. Wherein the color of the region of the target three-dimensional model which is not color filled is kept to be a first preset color.
The shale gas reservoir space distribution model obtained through the mode can intuitively reflect the distribution position and the size of the shale gas reservoir space in the target shale 90 through the internal color distribution. It is beneficial to the inspector to more efficiently and accurately judge whether the shale gas of the target shale 90 has exploitation value.
Referring to fig. 2, the application provides a shale gas reservoir space detection method, which is applied to an acquisition device 1, wherein the acquisition device 1 comprises a phased array ultrasonic probe 101 and is in communication connection with an upper computer 2;
the shale gas reservoir space detection method comprises the following steps:
A1. controlling the phased array ultrasonic probe 101 to move in the detection hole 91 extending along the surface of the target shale 90, and transversely scanning the target shale 90 in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections; the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data;
A2. generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section;
A3. screening and zeroing the data of each acquisition matrix;
A4. the acquisition matrix and the corresponding reflection position data after the screening and zeroing processing are sent to the upper computer 2, so that the upper computer 2 generates a shale gas reservoir space distribution model of the target shale 90 based on a reference two-dimensional model and a reference three-dimensional model of various typical structures of the shale which are established in advance according to the acquisition matrix and the reflection position data.
For the specific procedure of step A1 to step A4, reference may be made to the related process procedure described above.
Referring to fig. 3, the application provides a shale gas reservoir space detection method, which is applied to an upper computer 2, wherein the upper computer 2 is in communication connection with a collecting device 1, and the collecting device 1 comprises a phased array ultrasonic probe 101; the upper computer 2 stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the shale gas reservoir space detection method comprises the following steps:
B1. determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale 90 according to the acquisition matrix and the corresponding reflection position data from the acquisition device 1, and respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model; the acquisition matrix and the corresponding reflection position data are generated by the acquisition device 1 according to ultrasonic detection data of a scanning cross section obtained by transversely scanning the target shale 90 in the process that the phased array ultrasonic probe 101 moves in the detection hole 91 extending along the surface of the target shale 90;
B2. according to the numerical value of the acquisition matrix and the corresponding reflection position data, performing color filling on the target two-dimensional model corresponding to each scanning cross section to obtain a corresponding two-dimensional color map;
B3. The two-dimensional color map stack is filled into the target three-dimensional model, and a shale gas reservoir spatial distribution model of the target shale 90 is generated.
For the specific procedure of step B1 to step B3, reference may be made to the related process procedure described above.
Referring to fig. 4, fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present application, where the electronic device includes: a processor 301 and a memory 302, the processor 301 and the memory 302 being interconnected and in communication with each other by a communication bus 303 and/or other form of connection mechanism (not shown), the memory 302 storing a computer program executable by the processor 301, the processor 301 executing the computer program when the electronic device is running to perform the shale gas reservoir space detection method in any of the alternative implementations of the above embodiments to perform the following functions: controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale 90 in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections; the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data; generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section; screening and zeroing the data of each acquisition matrix; transmitting the acquisition matrix after screening and zeroing treatment and corresponding reflection position data to an upper computer, so that the upper computer generates a shale gas reservoir space distribution model of target shale based on a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale which are established in advance according to the acquisition matrix and the reflection position data; or the following functions are realized: determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to the acquisition matrix from the acquisition device and corresponding reflection position data, and respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model; the acquisition matrix and the corresponding reflection position data are generated by an acquisition device according to ultrasonic detection data of a scanning cross section obtained by transversely scanning the target shale in the process that the phased array ultrasonic probe moves in a detection hole extending along the surface of the target shale; according to the numerical value of the acquisition matrix and the corresponding reflection position data, performing color filling on the target two-dimensional model corresponding to each scanning cross section to obtain a corresponding two-dimensional color map; and filling the two-dimensional color map stack into the target three-dimensional model to generate a shale gas reservoir space distribution model of the target shale.
The embodiment of the application provides a computer readable storage medium, on which a computer program is stored, which when executed by a processor, performs the shale gas reservoir space detection method in any optional implementation of the above embodiment, so as to implement the following functions: controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale 90 in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections; the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data; generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section; screening and zeroing the data of each acquisition matrix; transmitting the acquisition matrix after screening and zeroing treatment and corresponding reflection position data to an upper computer, so that the upper computer generates a shale gas reservoir space distribution model of target shale based on a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale which are established in advance according to the acquisition matrix and the reflection position data; or the following functions are realized: determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to the acquisition matrix from the acquisition device and corresponding reflection position data, and respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model; the acquisition matrix and the corresponding reflection position data are generated by an acquisition device according to ultrasonic detection data of a scanning cross section obtained by transversely scanning the target shale in the process that the phased array ultrasonic probe moves in a detection hole extending along the surface of the target shale; according to the numerical value of the acquisition matrix and the corresponding reflection position data, performing color filling on the target two-dimensional model corresponding to each scanning cross section to obtain a corresponding two-dimensional color map; and filling the two-dimensional color map stack into the target three-dimensional model to generate a shale gas reservoir space distribution model of the target shale.
The computer readable storage medium may be implemented by any type or combination of volatile or non-volatile Memory devices, such as static random access Memory (Static Random Access Memory, SRAM), electrically erasable Programmable Read-Only Memory (EEPROM), erasable Programmable Read-Only Memory (Erasable Programmable Read Only Memory, EPROM), programmable Read-Only Memory (PROM), read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk, or optical disk.
In summary, the shale gas reservoir space detection method and the related equipment thereof provided by the application have the following advantages:
1. in the process of detecting the shale gas reservoir space in shale, ultrasonic detection is carried out through a phased array ultrasonic probe, and screening and zeroing treatment are carried out on a collected signal data set, so that the distribution condition of natural gas in shale can be known faster and better, and the requirements on ultrasonic detection equipment and wireless transmission equipment in terms of hardware are lower;
2. In the conventional method, in the process of detecting the storage space of shale gas, the integral distribution position of shale gas is detected (local distribution condition is not detected or low-accuracy detection is carried out), only one-dimensional detection signals are collected, all detection signals are reserved, and then data post-processing is carried out.
3. The shale type can be directly distinguished according to the detection data when the subsequent ultrasonic signal processing is carried out, and the two-dimensional and three-dimensional model is matched and called, so that the two-dimensional and three-dimensional modeling efficiency of the shale gas reservoir space distribution structure in the shale is greatly improved.
In the embodiments provided in the present application, it should be understood that the disclosed apparatus and method may be implemented in other manners. The above-described apparatus embodiments are merely illustrative, for example, the division of the units is merely a logical function division, and there may be other manners of division in actual implementation, and for example, multiple units or components may be combined or integrated into another system, or some features may be omitted, or not performed. Alternatively, the coupling or direct coupling or communication connection shown or discussed with each other may be through some communication interface, device or unit indirect coupling or communication connection, which may be in electrical, mechanical or other form.
Further, the units described as separate components may or may not be physically separate, and components shown as units may or may not be physical units, may be located in one place, or may be distributed over a plurality of network units. Some or all of the units may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
Furthermore, functional modules in various embodiments of the present application may be integrated together to form a single portion, or each module may exist alone, or two or more modules may be integrated to form a single portion.
In this document, relational terms such as first and second, and the like may be used solely to distinguish one entity or action from another entity or action without necessarily requiring or implying any actual such relationship or order between such entities or actions.
The foregoing is merely exemplary embodiments of the present application and is not intended to limit the scope of the present application, and various modifications and variations may be suggested to one skilled in the art. Any modification, equivalent replacement, improvement, etc. made within the spirit and principles of the present application should be included in the protection scope of the present application.
Claims (10)
1. The shale gas reservoir space detection system is characterized by comprising a collection device and an upper computer, wherein the collection device comprises a phased array ultrasonic probe; the acquisition device is in communication connection with the upper computer, and the upper computer stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the acquisition device is used for controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections, wherein the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data, an acquisition matrix of each scanning cross section is generated according to the reflected signal amplitude data of each scanning cross section, screening and zeroing processing is carried out on the data of each acquisition matrix, and the acquisition matrix and the corresponding reflected position data after the screening and zeroing processing are sent to the upper computer;
the upper computer is used for determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to the received acquisition matrix and the corresponding reflection position data, respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model, carrying out color filling on the target two-dimensional model corresponding to each scanning cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data to obtain corresponding two-dimensional color images, stacking and filling the two-dimensional color images into the target three-dimensional model, and generating a shale gas storage space distribution model of the target shale.
2. The shale gas reservoir space detection system of claim 1, wherein the ultrasonic detection data of each scanning cross section comprises one-dimensional reflected signal amplitude data obtained by multi-frame detection along different detection angles and reflection position data corresponding to each amplitude data in each frame of the reflected signal amplitude data;
the acquisition device performs the following steps when generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section:
taking the frame number of the reflected signal amplitude data of the same scanning cross section as the column number of the corresponding acquisition matrix, and determining the row number of the corresponding acquisition matrix according to the reflected signal amplitude data of the scanning cross section;
generating the acquisition matrix according to the number of rows and the number of columns and initializing the acquisition matrix to be a zero matrix;
and filling the amplitude data of the reflected signals of each frame into the data of each column of the acquisition matrix in sequence.
3. The shale gas reservoir space detection system of claim 2, wherein the acquisition device, when screening and zeroing the data of each acquisition matrix, performs:
Ascending order sorting is carried out on all the numerical values in the same acquisition matrix;
selecting a numerical value of a preset ranking in the acquisition matrix as a screening threshold of the acquisition matrix according to the ranking result;
and setting the numerical value smaller than the screening threshold value in the acquisition matrix to zero.
4. The shale gas reservoir space detection system of claim 2, wherein the acquisition device, when transmitting the acquisition matrix and corresponding reflected position data after screening and zeroing processing to the host computer, performs:
and transmitting the acquisition matrix after screening and zero setting processing and reflection position data corresponding to each non-zero value in the acquisition matrix to the upper computer.
5. The shale gas reservoir space detection system of claim 2, wherein the upper computer performs, when determining a reference two-dimensional model matching each of the scan cross sections and a reference three-dimensional model matching the target shale based on the received acquisition matrix and the corresponding reflection position data, respectively, as a target two-dimensional model and a target three-dimensional model:
acquiring edge contour lines of the scanning cross sections according to the received acquisition matrix and the corresponding reflection position data;
Determining a reference two-dimensional model matched with each scanning cross section from the reference two-dimensional models stored by the upper computer according to the edge contour line of each scanning cross section, and marking the reference two-dimensional model as a target two-dimensional model;
scaling the target two-dimensional model corresponding to each scanning cross section according to the edge contour line of each scanning cross section;
determining a reference three-dimensional model matched with the target shale from the reference three-dimensional models stored by the upper computer according to the target two-dimensional model subjected to scaling processing corresponding to each scanning cross section, and marking the reference three-dimensional model as a target three-dimensional model;
and scaling the target three-dimensional model according to the scaled target two-dimensional model corresponding to each scanning cross section.
6. The shale gas reservoir space detection system according to claim 2, wherein the upper computer performs color filling on the target two-dimensional model corresponding to each scanning cross section according to the numerical value of the acquisition matrix and the corresponding reflection position data to obtain a corresponding two-dimensional color map, and performs:
determining the color corresponding to each non-zero value according to the non-zero value of the acquisition matrix, and marking the color as a filling color;
And setting the color of the corresponding position point in the target two-dimensional model as the corresponding filling color according to the reflection position data corresponding to the non-zero numerical value of the acquisition matrix.
7. The shale gas reservoir space detection method is characterized by being applied to an acquisition device, wherein the acquisition device comprises a phased array ultrasonic probe and is in communication connection with an upper computer;
the shale gas reservoir space detection method comprises the following steps:
A1. controlling the phased array ultrasonic probe to move in a detection hole extending along the surface of the target shale, and transversely scanning the target shale in the moving process to acquire ultrasonic detection data of a plurality of scanning cross sections; the ultrasonic detection data comprise reflected signal amplitude data and corresponding reflected position data;
A2. generating an acquisition matrix of each scanning cross section according to the reflected signal amplitude data of each scanning cross section;
A3. screening and zeroing the data of each acquisition matrix;
A4. and sending the acquisition matrix after screening and zeroing treatment and corresponding reflection position data to the upper computer, so that the upper computer generates a shale gas reservoir space distribution model of the target shale based on a reference two-dimensional model and a reference three-dimensional model of various typical structures of the shale, which are established in advance, according to the acquisition matrix and the reflection position data.
8. The shale gas reservoir space detection method is characterized by being applied to an upper computer, wherein the upper computer is in communication connection with a collection device, and the collection device comprises a phased array ultrasonic probe; the upper computer stores a reference two-dimensional model and a reference three-dimensional model of various typical structures of shale;
the shale gas reservoir space detection method comprises the following steps:
B1. determining a reference two-dimensional model matched with each scanning cross section and a reference three-dimensional model matched with the target shale according to an acquisition matrix and corresponding reflection position data from the acquisition device, and respectively marking the reference two-dimensional model and the reference three-dimensional model as a target two-dimensional model and a target three-dimensional model; the acquisition matrix and the corresponding reflection position data are generated by the acquisition device according to ultrasonic detection data of a scanning cross section obtained by transversely scanning the target shale in the process that the phased array ultrasonic probe moves in a detection hole extending along the surface of the target shale;
B2. according to the numerical value of the acquisition matrix and the corresponding reflection position data, performing color filling on the target two-dimensional model corresponding to each scanning cross section to obtain a corresponding two-dimensional color map;
B3. and filling the two-dimensional color map stack into the target three-dimensional model to generate a shale gas reservoir space distribution model of the target shale.
9. An electronic device comprising a processor and a memory, the memory storing a computer program executable by the processor, when executing the computer program, running the steps in the shale gas reservoir space detection method of any of claims 7-8.
10. A computer readable storage medium having stored thereon a computer program, which when executed by a processor performs the steps in the shale gas reservoir space detection method as claimed in any of claims 7-8.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311770787.5A CN117451857B (en) | 2023-12-21 | 2023-12-21 | Shale gas reservoir space detection method and related equipment thereof |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202311770787.5A CN117451857B (en) | 2023-12-21 | 2023-12-21 | Shale gas reservoir space detection method and related equipment thereof |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117451857A true CN117451857A (en) | 2024-01-26 |
| CN117451857B CN117451857B (en) | 2024-03-08 |
Family
ID=89584020
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202311770787.5A Active CN117451857B (en) | 2023-12-21 | 2023-12-21 | Shale gas reservoir space detection method and related equipment thereof |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117451857B (en) |
Family Cites Families (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106846458A (en) * | 2016-12-15 | 2017-06-13 | 哈尔滨工业大学（威海） | Stereoscopic ultrasonic model building method and device based on 3D printing |
| CN109490165A (en) * | 2018-11-02 | 2019-03-19 | 中国石油化工股份有限公司 | The method for characterizing the non-fabric selectivity reservoir space of carbonate rock |
| CN111608654A (en) * | 2020-06-05 | 2020-09-01 | 中国石油大学(华东) | Shale reservoir natural frequency measuring method and device |
| CN112627809A (en) * | 2020-12-18 | 2021-04-09 | 中国矿业大学(北京) | Ultrasonic imaging logging device and coal seam crack prediction method |
| CN115929282A (en) * | 2021-09-22 | 2023-04-07 | 中国石油化工股份有限公司 | Logging system and method based on orientation while drilling gamma and ultrasonic imaging |
| CN114897767A (en) * | 2022-03-15 | 2022-08-12 | 中海油能源发展股份有限公司 | Multi-scale representation and reservoir classification method for reservoir space of compact mixed-reservoir rock reservoir |
| CN115078422B (en) * | 2022-06-13 | 2024-07-26 | 中国地质调查局油气资源调查中心 | Method for establishing three-dimensional mineral model |
| CN116906031A (en) * | 2023-07-07 | 2023-10-20 | 中铁第四勘察设计院集团有限公司 | Three-dimensional high-precision detection method for surrounding rock of drilling hole |
| CN116976519A (en) * | 2023-08-03 | 2023-10-31 | 中国地质大学（北京） | Shale oil reservoir single well recoverable reserve prediction method and system |
-
2023
- 2023-12-21 CN CN202311770787.5A patent/CN117451857B/en active Active
Also Published As
| Publication number | Publication date |
|---|---|
| CN117451857B (en) | 2024-03-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bornaz et al. | Terrestrial laser scanner data processing | |
| Tang et al. | Quantification of edge loss of laser scanned data at spatial discontinuities | |
| Barbarella et al. | Monitoring of large landslides by Terrestrial Laser Scanning techniques: field data collection and processing | |
| EP2304688B1 (en) | Automated building outline detection | |
| Lindenbergh et al. | Accuracy analysis of the Leica HDS3000 and feasibility of tunnel deformation monitoring | |
| Drews et al. | Validation of fracture data recognition in rock masses by automated plane detection in 3D point clouds | |
| CN104569972B (en) | Plant root system three-dimensional configuration nondestructive testing method | |
| RU2364887C2 (en) | Method for navigation of aircraft by radar images of earth surface with application of digital area models | |
| US11604301B2 (en) | Methods and systems for automated sonic imaging | |
| CN105761308A (en) | Ground LiDAR and image data fused occlusion region building facade reconstruction method | |
| CN108074232B (en) | Voxel segmentation-based airborne LIDAR building detection method | |
| CN108663041A (en) | A kind of method and device for drawing navigation map | |
| CN112068153A (en) | Crown clearance rate estimation method based on foundation laser radar point cloud | |
| CN116310177B (en) | Soil map collaborative drawing method under three-dimensional geological landform scanning information | |
| NL2020152A (en) | Kirchhoff Beam Migration Method Under Complex Topography | |
| BRPI0413671B1 (en) | lighting monitoring process to make filling decisions | |
| CN112581456B (en) | Image crack real-time extraction method suitable for ultrasonic imaging logging | |
| Singh et al. | Moving from plot-based to hillslope-scale assessments of savanna vegetation structure with long-range terrestrial laser scanning (LR-TLS) | |
| CN1684105A (en) | Automatic registration method for large-scale three-dimensional scene multi-view laser scanning data | |
| Wang et al. | Fine detection technology of rock mass structure based on borehole acousto-optic combined measurement | |
| CN117451857B (en) | Shale gas reservoir space detection method and related equipment thereof | |
| Milde et al. | Building reconstruction using a structural description based on a formal grammar | |
| Calvo et al. | Unlocking the correlation in fluvial outcrops by using a DOM-derived virtual datum: Method description and field tests in the Huesca fluvial fan, Ebro Basin (Spain) | |
| Lin et al. | $ k $-segments-based geometric modeling of VLS scan lines | |
| Wang et al. | Improved SSD Framework for Automatic Subsurface Object Indentification for GPR Data Processing |
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 |