CN115903032A - Underground fluid storage space and migration channel detection method and device - Google Patents

Underground fluid storage space and migration channel detection method and device Download PDF

Info

Publication number
CN115903032A
CN115903032A CN202211307019.1A CN202211307019A CN115903032A CN 115903032 A CN115903032 A CN 115903032A CN 202211307019 A CN202211307019 A CN 202211307019A CN 115903032 A CN115903032 A CN 115903032A
Authority
CN
China
Prior art keywords
matrix
data
storage space
fluid storage
migration
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202211307019.1A
Other languages
Chinese (zh)
Other versions
CN115903032B (en
Inventor
向阳
彭苏萍
林朋
李闯建
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Mining and Technology Beijing CUMTB
Original Assignee
China University of Mining and Technology Beijing CUMTB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Mining and Technology Beijing CUMTB filed Critical China University of Mining and Technology Beijing CUMTB
Priority to CN202211307019.1A priority Critical patent/CN115903032B/en
Publication of CN115903032A publication Critical patent/CN115903032A/en
Application granted granted Critical
Publication of CN115903032B publication Critical patent/CN115903032B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

The invention belongs to the technical field of detection, and relates to a method and a device for detecting an underground fluid storage space and a migration channel, wherein common offset seismic wave field data of a region to be processed are collected firstly; processing the common offset seismic wave field data to obtain scattered wave data; then, performing homing imaging on scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, and realizing the detection of the underground fluid storage space and the migration channel; the method takes the kinematic characteristic difference of the wave field and the dynamic characteristic of the wave field into consideration, so that the wave field information can be extracted more quickly and conveniently, the matrix decomposition technology is applied, the scattered wave data of discontinuous geologic bodies can be accurately acquired, the calculation efficiency and the accuracy are high, the method is beneficial to positioning underground abnormal geologic bodies, and the method is convenient for geothermal heating, mineral resources and CO 2 And items such as sealing and the like are safely and effectively promoted.

Description

Underground fluid storage space and migration channel detection method and device
The technical field is as follows:
the invention belongs to the technical field of detection, and relates to a method and a device for detecting an underground fluid storage space and a migration channel.
Background art:
underground discontinuous geological anomalous bodies have important influence on the development and utilization of geothermy, mineral resources and the like. The discontinuous geological abnormal body can be used as a geothermal storage space heat conduction channel for heat conduction, and is of great importance for geothermal resource utilization and well drilling position selection. The underground discontinuous fractures serving as storage spaces and migration channels of mineral resources such as coal bed gas, petroleum and natural gas play a vital role in selecting drilling positions and exploring, developing and utilizing oil and gas resources. The scattering data is the wave field response of the underground geological abnormal body, can be used for accurately positioning the non-uniform discontinuous body and provides stronger underground space illumination, and the premise of using the scattering data is to extract the scattering data from the wave field data and then carry out high-precision imaging to accurately position and image the geological abnormal body.
In the existing method, a plane wave decomposition method is usually adopted to obtain scattering data, but the method only considers the kinematic characteristic difference of a wave field, ignores the dynamic characteristic of the wave field and is not beneficial to the extraction of wave field information.
The invention content is as follows:
the invention aims to overcome the defects in the prior art, and provides a method and a device for detecting an underground fluid storage space and a migration channel.
In order to achieve the purpose, the specific process of detecting the underground fluid storage space and the migration channel comprises the following steps:
(1) Acquiring common offset distance seismic wave field data of a region to be processed;
(2) Processing the common offset seismic wave field data to obtain scattered wave data;
(3) And performing homing imaging on the scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, so as to realize detection of the underground fluid storage space and the migration channel.
As a further technical scheme of the invention, the specific process of acquiring the scattered wave data in the step (2) is as follows:
(21) Constructing frequency space domain seismic data by utilizing Fourier transform on the common offset seismic wave data acquired in the step (1), and converting the frequency space domain seismic data into a Hankel matrix;
(22) Selecting a plurality of rows and columns from a Hankel matrix according to a matrix row-column contribution value weight function to construct a decomposition submatrix;
(23) And approximating the low-rank information of the Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component, and acquiring scattered wave data.
As a further technical scheme of the invention, the Hankle matrix obtained in the step (21) is as follows:
Figure BDA0003906240050000021
/>
wherein d is i For the ith channel of data at a certain frequency, m and n are the rows and columns of the matrix.
As a further technical solution of the present invention, the matrix row and column contribution weight functions in step (22) are respectively:
Figure BDA0003906240050000022
Figure BDA0003906240050000023
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,
Figure BDA0003906240050000024
for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
Figure BDA0003906240050000025
Figure BDA0003906240050000026
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row and column contribution weight function.
As a further technical solution of the present invention, the approximate low rank component shown in step (23)
Figure BDA0003906240050000027
Scattered wave data->
Figure BDA0003906240050000028
As a further technical scheme of the invention, the shift algorithm adopted in the step (3) is a Kirchoff shift algorithm.
Compared with the prior art, the method simultaneously considers the kinematic characteristic difference of the wave field and the dynamic characteristic of the wave field, so that the wave field information is extracted more quickly and conveniently, the matrix decomposition technology is used, the scattered wave data of discontinuous geologic bodies can be accurately obtained, the calculation efficiency and the accuracy are high, the method is favorable for positioning underground abnormal geologic bodies, and is convenient for geothermal heat, mineral resources and CO 2 And items such as sealing and the like are safely and effectively promoted.
Description of the drawings:
FIG. 1 is a block diagram of the working process of the present invention.
Fig. 2 is a block diagram of a process for acquiring scattered wave data according to the present invention.
FIG. 3 is a block diagram of the underground fluid storage space and migration channel detection apparatus according to the present invention.
FIG. 4 is a raw seismic stack data record of common-offset seismic wavefield data of example 1 of the present invention.
FIG. 5 is a graph of scattered wave data extracted in embodiment 1 of the present invention.
FIG. 6 is a diagram showing the result of the scattered wave data imaging in example 1 of the present invention.
The specific implementation mode is as follows:
to make the objects, technical solutions and advantages of the embodiments of the present invention clearer and more complete, the technical solutions of the present invention will be described below with reference to the accompanying drawings, it is obvious that the described embodiments are some, but not all embodiments of the present invention, and all other embodiments obtained by a person skilled in the art without making creative efforts based on the embodiments of the present invention belong to the protection scope of the present invention.
Example 1:
as shown in fig. 1 and fig. 2, the present embodiment provides a method for detecting an underground fluid storage space and a migration channel, which includes the following specific processes:
(1) Determining a region to be processed, and acquiring common offset seismic wave field data of the region to be processed, as shown in figure 1;
(2) Constructing frequency-space domain seismic data by utilizing Fourier transform according to the common offset seismic wave field data acquired in the step (1), and converting the frequency-space domain seismic data into a Hankel matrix H:
Figure BDA0003906240050000031
wherein d is i The ith channel of data under a certain frequency, and m and n are rows and columns of a matrix;
(3) Selecting a plurality of rows and columns from a Hankel matrix to construct a decomposition submatrix according to a matrix row and column contribution value weight function, wherein the matrix row and column contribution weight function is respectively as follows:
Figure BDA0003906240050000041
Figure BDA0003906240050000042
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,
Figure BDA0003906240050000043
for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
Figure BDA0003906240050000044
Figure BDA0003906240050000045
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method comprises the following steps of obtaining three matrixes A, U and B through matrix row and column contribution weight functions as a generalized inverse matrix;
(4) According to the submatrix obtained in the step (3), utilizing a matrix decomposition technology to approximate Hankle matrix low-rank information to obtain approximate low-rank components
Figure BDA0003906240050000046
Thereby acquiring scattered wave data->
Figure BDA0003906240050000047
As shown in fig. 5;
(5) And (3) performing homing imaging on the scattered wave data obtained in the step (4) by using a Kirchhoff migration algorithm, and obtaining a fine positioning result of the underground fluid storage space and the migration channel as shown in fig. 6 to realize detection of the underground fluid storage space and the migration channel.
Example 2:
as shown in fig. 3, the present embodiment provides an underground fluid storage space and migration channel detection apparatus, including:
the data acquisition module is used for acquiring common offset seismic wave field data of the area to be detected;
the data processing module is used for processing the common offset seismic wave field data to obtain scattered wave data;
and the imaging module is used for carrying out offset imaging on the scattered wave data to obtain a fine detection result of the underground fluid storage space and the migration channel.
Specifically, the data processing module comprises a Hankel matrix conversion unit, a sub-matrix decomposition unit and a scattered wave data acquisition unit, the Hankel matrix conversion unit constructs frequency space domain seismic data from the collected common offset seismic wave field data by utilizing Fourier transform and converts the frequency space domain seismic data into a Hankel matrix H, and the sub-matrix decomposition unit selects a plurality of rows and columns from the Hankel matrix to construct a decomposition sub-matrix and construct a target function by utilizing a matrix row-column contribution value weight function; the scattered wave data acquisition unit approximates the low-rank information of a Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component, so that the scattered wave data is acquired.
More specifically, the Hankel matrix H obtained by the Hankel matrix conversion unit is:
Figure BDA0003906240050000051
wherein d is i The ith channel of data under a certain frequency, and m and n are rows and columns of a matrix;
matrix row and column contribution weight functions adopted by the sub-matrix decomposition unit are respectively as follows:
Figure BDA0003906240050000052
Figure BDA0003906240050000053
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,
Figure BDA0003906240050000054
for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
Figure BDA0003906240050000055
Figure BDA0003906240050000056
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and the matrix U is formed by intersecting elements of the A and the B, namely:
U=A + HB +
wherein, (.) + The matrix is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row-column contribution weight function;
approximate low rank component obtained by scattered wave acquisition unit
Figure BDA0003906240050000057
Thereby obtaining scattered wave data
Figure BDA0003906240050000058
In the embodiment, the common offset seismic wave field data of the area to be processed is obtained through a data acquisition module; the data processing module constructs frequency-space domain seismic data by utilizing Fourier transform and converts the frequency-space domain seismic data into a Hankel matrix; selecting a plurality of rows and columns from the Hankel matrix according to the matrix row and column contribution value weight function to construct a decomposition submatrix; then, approximating Hankle matrix low-rank information by using a matrix decomposition technology to obtain scattering wave field data; and finally, the imaging module performs homing imaging on the scattered wave data through a migration algorithm to obtain detection results of the underground fluid storage space and the migration channel.
The computer program product of the method and the device for finely detecting the underground fluid storage space and the migration channel provided by the embodiment of the invention comprises a computer readable storage medium storing program codes, wherein instructions included in the program codes can be used for executing the method described in the previous embodiment, specific implementation can refer to the method embodiment, and detailed description is omitted, and algorithms not described in detail are all general technologies in the field.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk, or an optical disk, and various media capable of storing program codes.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: those skilled in the art can still make modifications or changes to the embodiments described in the foregoing embodiments, or make equivalent substitutions for some features, within the scope of the disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (6)

1. A method for detecting underground fluid storage space and migration channels is characterized by comprising the following specific processes:
(1) Acquiring common offset distance seismic wave field data of a region to be processed;
(2) Processing the common offset seismic wave field data to obtain scattered wave data;
(3) And performing homing imaging on the scattered wave data through a migration algorithm to obtain a fine positioning result of the underground fluid storage space and the migration channel, so as to realize detection of the underground fluid storage space and the migration channel.
2. The method for detecting the underground fluid storage space and the migration channel according to claim 1, wherein the step (2) of acquiring the scattered wave data comprises the following specific steps:
(21) Constructing frequency space domain seismic data by utilizing Fourier transform on the common offset seismic wave data acquired in the step (1), and converting the frequency space domain seismic data into a Hankel matrix;
(22) Selecting a plurality of rows and columns from the Hankel matrix according to a matrix row and column contribution value weight function to construct a decomposition submatrix;
(23) And approximating the low-rank information of the Hankle matrix by using a matrix decomposition technology to obtain an approximate low-rank component and acquiring scattered wave data.
3. The method for detecting the underground fluid storage space and the migration channel according to claim 2, wherein the Hankle matrix obtained in the step (21) is:
Figure FDA0003906240040000011
wherein d is i For the ith channel of data at a certain frequency, m and n are the rows and columns of the matrix.
4. A method according to claim 3, wherein the matrix row and column contribution weighting functions of step (22) are:
Figure FDA0003906240040000012
Figure FDA0003906240040000013
wherein i and j are the numbers of indexes in rows and columns of the Hankle matrix, x represents the number of columns, y represents the number of rows,
Figure FDA0003906240040000016
for the contribution value function of the component, the weight value of each row and column can be circularly calculated by the formula, and the column vector with larger contribution weight is selected to construct decomposed matrixes A and B based on the weight function, wherein the matrix A is composed of column vectors, the matrix B is composed of row vectors, and the screening principle of the column vectors and the row vectors is as follows:
Figure FDA0003906240040000014
Figure FDA0003906240040000015
according to the formula, c column vectors and r row vectors with larger weights are respectively selected from each row vector and column vector to construct matrixes A and B, and a matrix U is formed by intersecting elements of A and B, namely:
U=A + HB +
wherein, (.) + The method is a generalized inverse matrix, and three matrixes A, U and B are obtained through a matrix row and column contribution weight function.
5. The method of claim 4, wherein the approximate low rank component of step (23) is
Figure FDA0003906240040000021
Scattered wave data->
Figure FDA0003906240040000022
/>
6. The method for detecting underground fluid storage space and migration passage according to claim 1, wherein the migration algorithm used in step (3) is Kirchhoff migration algorithm.
CN202211307019.1A 2022-10-25 2022-10-25 Underground fluid storage space and migration channel detection method and device Active CN115903032B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211307019.1A CN115903032B (en) 2022-10-25 2022-10-25 Underground fluid storage space and migration channel detection method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211307019.1A CN115903032B (en) 2022-10-25 2022-10-25 Underground fluid storage space and migration channel detection method and device

Publications (2)

Publication Number Publication Date
CN115903032A true CN115903032A (en) 2023-04-04
CN115903032B CN115903032B (en) 2023-12-15

Family

ID=86470065

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211307019.1A Active CN115903032B (en) 2022-10-25 2022-10-25 Underground fluid storage space and migration channel detection method and device

Country Status (1)

Country Link
CN (1) CN115903032B (en)

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018129844A1 (en) * 2017-01-10 2018-07-19 中国科学院地质与地球物理研究所 Seismic diffracted wave separation method and device
CN111999770A (en) * 2020-09-03 2020-11-27 中国地质科学院 Precise beam offset imaging method and system for converting TTI medium into PS wave
CN112684504A (en) * 2020-12-03 2021-04-20 中国矿业大学 Urban underground cavity rapid detection method based on total scattering model
CN113945967A (en) * 2021-10-14 2022-01-18 中国矿业大学(北京) Diffraction wave separation method and device
CN113945969A (en) * 2021-10-19 2022-01-18 中国矿业大学(北京) Discontinuous geologic body detection method and device and electronic equipment
CN114924315A (en) * 2022-05-10 2022-08-19 北京世纪金道石油技术开发有限公司 Polymorphic seismic exploration method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2018129844A1 (en) * 2017-01-10 2018-07-19 中国科学院地质与地球物理研究所 Seismic diffracted wave separation method and device
CN111999770A (en) * 2020-09-03 2020-11-27 中国地质科学院 Precise beam offset imaging method and system for converting TTI medium into PS wave
CN112684504A (en) * 2020-12-03 2021-04-20 中国矿业大学 Urban underground cavity rapid detection method based on total scattering model
CN113945967A (en) * 2021-10-14 2022-01-18 中国矿业大学(北京) Diffraction wave separation method and device
CN113945969A (en) * 2021-10-19 2022-01-18 中国矿业大学(北京) Discontinuous geologic body detection method and device and electronic equipment
CN114924315A (en) * 2022-05-10 2022-08-19 北京世纪金道石油技术开发有限公司 Polymorphic seismic exploration method and system

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
史恒;刘洪;高星;丁仁伟;王伟;: "散射波叠前时间偏移成像", 石油地球物理勘探, no. 04 *
李纯阳,周官群,章俊: "散射波成像技术在地铁空洞探测中的应用", 安徽理工大学学报( 自然科学版), vol. 38, no. 3, pages 40 - 44 *
王瑞贞;白旭明;王金宽;赵利庆;韩力;邱文平;: "朝克乌拉凹陷火成岩覆盖区地震勘探方法研究", 石油物探, no. 03 *

Also Published As

Publication number Publication date
CN115903032B (en) 2023-12-15

Similar Documents

Publication Publication Date Title
Simmons et al. LLNL‐G3Dv3: Global P wave tomography model for improved regional and teleseismic travel time prediction
Afonso et al. 3‐D multi‐observable probabilistic inversion for the compositional and thermal structure of the lithosphere and upper mantle. II: General methodology and resolution analysis
Xiao et al. Siamese earthquake transformer: A pair‐input deep‐learning model for earthquake detection and phase picking on a seismic array
Yang et al. Dual-mode FPGA implementation of target and anomaly detection algorithms for real-time hyperspectral imaging
Hu et al. Using deep learning to derive shear‐wave velocity models from surface‐wave dispersion data
Zhang et al. Extracting dispersion curves from ambient noise correlations using deep learning
US10436925B2 (en) Method for processing marine seismic data
Zhang et al. Moho depth variations from receiver function imaging in the northeastern North China Craton and its tectonic implications
Yuan et al. A probabilistic shear wave velocity model of the crust in the central West Australian craton constrained by transdimensional inversion of ambient noise dispersion
Kumar et al. Dynamic time warping as an alternative to windowed cross correlation in seismological applications
CN112114362A (en) Method for reconstructing space-time field of underground shallow layer explosion
Sopher Converting scanned images of seismic reflection data into SEG-Y format
Chevrot et al. Optimized discrete wavelet transforms in the cubed sphere with the lifting scheme—implications for global finite-frequency tomography
Li Recent advances in earthquake monitoring II: Emergence of next-generation intelligent systems
Wang et al. Denoising with weak signal preservation by group-sparsity transform learning
Guevara et al. A data-driven workflow for predicting horizontal well production using vertical well logs
López et al. Spectral gap-based seismic survey design
Ojeda G et al. Spatio‐temporal entropy analysis of the magnetic field to help magnetic cloud characterization
AlSalmi et al. Automated seismic semantic segmentation using attention U-Net
CN115903032A (en) Underground fluid storage space and migration channel detection method and device
Gu et al. A modified excitation amplitude imaging condition for prestack reverse time migration
Liu et al. Line survey joint denoising via low-rank minimization
Song et al. Domain adaptation in automatic picking of phase velocity dispersions based on deep learning
CN114185095B (en) Method for suppressing multiple waves of three-dimensional plane wave domain seismic data
CN113484906B (en) High-resolution energy field reconstruction method based on low-frequency energy spectrum data driving

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