CN113281810B - Volcanic channel identification and division method - Google Patents
Volcanic channel identification and division method Download PDFInfo
- Publication number
- CN113281810B CN113281810B CN202110496212.3A CN202110496212A CN113281810B CN 113281810 B CN113281810 B CN 113281810B CN 202110496212 A CN202110496212 A CN 202110496212A CN 113281810 B CN113281810 B CN 113281810B
- Authority
- CN
- China
- Prior art keywords
- volcanic
- channel
- seismic data
- coherence
- anomaly
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 35
- 238000012545 processing Methods 0.000 claims abstract description 50
- 230000001427 coherent effect Effects 0.000 claims abstract description 43
- 230000002159 abnormal effect Effects 0.000 claims abstract description 38
- 230000005856 abnormality Effects 0.000 claims abstract description 31
- 239000000758 substrate Substances 0.000 claims abstract description 23
- 230000005484 gravity Effects 0.000 claims abstract description 16
- 230000035515 penetration Effects 0.000 claims abstract description 9
- 239000002002 slurry Substances 0.000 claims description 23
- 239000011435 rock Substances 0.000 claims description 22
- 238000004321 preservation Methods 0.000 claims description 6
- 230000004044 response Effects 0.000 claims description 6
- 230000009545 invasion Effects 0.000 claims description 4
- 238000004880 explosion Methods 0.000 claims description 2
- 238000005553 drilling Methods 0.000 abstract description 9
- 230000007246 mechanism Effects 0.000 description 6
- 238000011161 development Methods 0.000 description 5
- CCEKAJIANROZEO-UHFFFAOYSA-N sulfluramid Chemical group CCNS(=O)(=O)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)C(F)(F)F CCEKAJIANROZEO-UHFFFAOYSA-N 0.000 description 5
- 238000010586 diagram Methods 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 230000002349 favourable effect Effects 0.000 description 3
- 238000001914 filtration Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 230000000717 retained effect Effects 0.000 description 3
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000005611 electricity Effects 0.000 description 2
- 230000005389 magnetism Effects 0.000 description 2
- 239000011324 bead Substances 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000012512 characterization method Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 239000002360 explosive Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/30—Analysis
- G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
- G01V3/08—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with magnetic or electric fields produced or modified by objects or geological structures or by detecting devices
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V7/00—Measuring gravitational fields or waves; Gravimetric prospecting or detecting
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Electromagnetism (AREA)
- Acoustics & Sound (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention discloses a volcanic channel identification and division method, which belongs to the field of geophysical seismic interpretation and is characterized by comprising the following steps of: a. processing the high-resolution electrical data; b. identifying a coherence anomaly within a coherence attribute volume; c. vertically identifying high-resistance abnormal strips; d. the seismic data is a volcanic rock-magma channel when the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and correspond to the magnetic anomaly, and the seismic data is a volcanic hydrothermal channel when the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and have no magnetic anomaly; e. the volcanic rock-magma channel with the abnormal magnetic property, electric property and earthquake in the threshold value is superposed with the gravity abnormality, the volcanic rock-magma channel with the abnormal coherence and the abnormal substrate fracture superposition and the plane area larger than ten square kilometers is a penetration type volcanic rock-magma channel, and the rest is a crack type volcanic rock-magma channel. The method can reduce the multi-solution of single seismic data describing the volcanic tunnel, and provides accurate theoretical support for the implementation of a subsequent well drilling sidetrack drilling scheme.
Description
Technical Field
The invention relates to the technical field of geophysical seismic interpretation, in particular to a volcanic channel identification and division method.
Background
The outbreak phase development of the two-fold unified volcanic rock under the Sichuan basin is a newly discovered important successive layer system for oil and gas exploration, the conventional research degree is not high, drilling and core data are lacked, single technical means are different in magnitude, and the difficulty in describing volcanic channels is high, so that the volcanic mechanism and reservoir formation mechanism research is weak, the understanding of reservoir distribution and reservoir formation evolution is insufficient, and the drilling verification of the related mode and experience of the peripheral basin is used for reference, so that the volcanic rock is not suitable for the volcanic rock exploration of the Sichuan basin. Therefore, the combined identification of multiple technical means of gravity electromagnetic vibration, the refinement of volcanic channels, the improvement of the overall research of the basin grade related to volcanic mechanisms and the support for expanding the new breakthrough of volcanic rock exploration are urgently needed to be developed.
The Chinese patent document with the publication number of CN 112114358A and the publication date of 2020, 12 and 22 discloses an underground volcano channel identification method based on three-dimensional seismic data representation, which is characterized by comprising the following steps of:
step 1: performing anisotropic diffusion filtering on the stacked three-dimensional seismic data volume to obtain a seismic data volume after diffusion filtering;
and 2, step: acquiring a similarity data volume based on the seismic data volume after diffusion filtering;
and 3, step 3: carrying out binarization processing on the similarity data volume to obtain a binary attribute data volume;
and 4, step 4: calculating a binary attribute accumulated value below a target layer based on the binary attribute data volume;
and 5: performing gain calculation on the binary attribute accumulated value to obtain a post-gain attribute;
step 6: searching the gained attributes, determining a local most value point set, and verifying an underground volcanic channel prediction result through the local most value point set to obtain an underground volcanic channel coordinate set;
and 7: and obtaining the spatial distribution of the underground volcanic channel point set based on the underground volcanic channel coordinate set and the binary attribute data volume.
The patent document discloses a method for identifying underground volcanic tunnels based on three-dimensional seismic data characterization, although the development position of the volcanic tunnels can be predicted. However, since the seismic data are single, the influence of the resolution of the seismic data is large, the multi-solution of the result cannot be effectively solved, and the recognition accuracy is low.
Disclosure of Invention
The invention provides a volcanic channel identification and division method for overcoming the defects of the prior art, and the method divides different volcanic channel types from different dimensions by organically combining four kinds of data, namely heavy data, magnetic data, electric data and vibration data, thereby greatly improving the identification accuracy; the method can reduce the ambiguity of depicting the volcanic channel by a single seismic data, can clearly determine the eruption type of the two-fold volcanic rock and the mode of the corresponding volcanic mechanism, realizes the development and plane distribution rule of the reservoir favorable for outburst of the volcanic rock on the basis, and can reasonably explain the unexplained reason of the well-drilled reservoir, thereby providing accurate theoretical support for the implementation of a subsequent well-drilling sidetracking scheme.
The invention is realized by the following technical scheme:
a volcanic channel identification and division method is characterized by comprising the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile;
d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 6-35Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In step b, the depth range refers to a range of the threshold value being 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherent anomalies in the coherent attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
In the step c, the vertical identification means that residual magnetic anomaly plane distribution is superposed according to a high-resistance piercing threshold value of 8000-11000 omega-m, and the volcanic channel response is determined when the high-resistance strip is transversely larger than 1km and magnetic anomalies exist.
In the step d, comparing specifically means comparing various earthquake coherence anomalies in the reserved effective depth range with the electrical profile, and if magnetic anomalies exist and the width is greater than 1km and corresponds to a strong coherence anomaly on the earthquake profile, dividing the earthquake coherence anomalies into volcanic explosion magma channels; if the longitudinal depth is below the target layer, dividing the target layer into volcano invasion rock-slurry channels;
in the step e, the rest of the gap type volcanic rock slurry channels are strips, wherein the volcanic rock slurry channels with abnormal magnetic property, electric property and earthquake in the threshold value are superposed with gravity abnormity, the coherence abnormity is superposed with the substrate fracture, and the plane area is less than or equal to ten square kilometers.
The crack type volcanic rock-slurry channel comprises a linear volcanic rock-slurry channel, a point volcanic rock-slurry channel and a beaded volcanic rock-slurry channel.
The beneficial effects of the invention are mainly shown in the following aspects:
1. according to the method, a, fidelity and amplitude-preserving seismic processing is carried out on seismic data based on single-point cutting, and high-resolution electrical data are processed; b. on the basis of seismic data of fidelity amplitude, a coherence attribute body in the depth range from a base to a target layer is calculated, and columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body are identified; c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile; d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel; e. compared with the prior art, different volcanic channel types are divided from different dimensions by organically combining four data of gravity, magnetism, electricity and earthquake, so that the identification accuracy is greatly improved; the method can reduce the ambiguity of single seismic data depicting volcanic channels, can clarify the eruption type of the two-fold volcanic rock and the mode of the corresponding volcanic mechanism, realizes the development and plane distribution rule of the favorable outburst phase volcanic rock reservoir on the basis, and can reasonably explain the reasons of the undeveloped well reservoir, thereby providing accurate theoretical support for the implementation of the subsequent well drilling sidetracking scheme.
2. According to the invention, the volcanic channels are identified and accurately divided, so that a constructive guiding significance is provided for adjustment of the exploration direction and deployment of the drilling well.
3. The method is particularly suitable for identifying the second-folded volcanic rock under the Sichuan basin, can identify that the volcanic rock channels in the West region of Sichuan are mainly in a crack type bead-like shape, and parts of treetop-shaped and punctiform explosive phase reservoirs are distributed along the direction close to south and north of the crack, so that the planar distribution range of the volcanic channels can be well identified, and the longitudinal distribution characteristics of the volcanic channels can be accurately depicted.
Drawings
The invention will be further described in detail with reference to the drawings and the following detailed description:
FIG. 1 is a block flow diagram of the present invention;
FIG. 2 is a block diagram of a process flow for fidelity and amplitude-preserving seismic data in accordance with the present invention;
FIG. 3 is a high resolution time-frequency electromagnetic data processing flow of the present invention;
FIG. 4 is a schematic view of volcanic tunnels of different plan shapes and sizes in accordance with the present invention;
FIG. 5 is a diagram illustrating the calculation of coherence values according to the present invention.
Detailed Description
Example 1
Referring to fig. 1 to 5, a volcanic tunnel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on a high-resolution electrical profile;
d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
a. Performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data; b. on the basis of seismic data of fidelity amplitude, a coherence attribute body in the depth range from a base to a target layer is calculated, and columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body are identified; c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile; d. through comparison, if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and correspond to the magnetic anomaly, the seismic data are divided into volcanic magma channels, and if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and no magnetic anomaly exists, the seismic data are divided into volcanic hydrothermal channels; e. compared with the prior art, different volcanic channel types are divided from different dimensions by organically combining four data of gravity, magnetism, electricity and earthquake, so that the identification accuracy is greatly improved; the method can reduce the ambiguity of single seismic data depicting volcanic channels, can clarify the eruption type of the two-fold volcanic rock and the mode of the corresponding volcanic mechanism, realizes the development and plane distribution rule of the favorable outburst phase volcanic rock reservoir on the basis, and can reasonably explain the reasons of the undeveloped well reservoir, thereby providing accurate theoretical support for the implementation of the subsequent well drilling sidetracking scheme.
Example 2
Referring to fig. 1 to 5, a volcanic channel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on a high-resolution electrical profile;
d. through comparison, if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and correspond to the magnetic anomaly, the seismic data are divided into volcanic magma channels, and if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and no magnetic anomaly exists, the seismic data are divided into volcanic hydrothermal channels;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 6Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In the step b, the depth range refers to a range with a threshold value of 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
Example 3
Referring to fig. 1 to 5, a volcanic channel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of seismic data of fidelity amplitude, a coherence attribute body in the depth range from a base to a target layer is calculated, and columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body are identified;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile;
d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 10Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In the step b, the depth range refers to a range with a threshold value of 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
In the step c, the vertical identification means that residual magnetic anomaly plane distribution is superposed according to a high-resistance piercing threshold value of 8000 omega-m, and the volcanic channel response is determined when the high-resistance strip is larger than 1km transversely and has magnetic anomaly.
The volcanic channels are identified and accurately divided, so that the adjustment of the drilling exploration direction and the deployment is facilitated to provide a constructive guiding significance.
Example 4
Referring to fig. 1 to 5, a volcanic tunnel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile;
d. through comparison, if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and correspond to the magnetic anomaly, the seismic data are divided into volcanic magma channels, and if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and no magnetic anomaly exists, the seismic data are divided into volcanic hydrothermal channels;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 20Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In step b, the depth range refers to a range of the threshold value being 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
In the step c, the vertical identification refers to superposing the residual magnetic anomaly plane distribution according to a high-resistance piercing threshold value of 9000 omega. M, and determining that the high-resistance strip with the transverse direction being greater than 1km and the magnetic anomaly exists as a volcanic channel response.
In the step d, the comparison specifically means that various earthquake coherence anomalies within an effective depth range are retained and compared with the electrical profile, and if magnetic anomalies exist and the width is greater than 1km and the seismic coherence anomalies on the seismic profile correspond to, the seismic coherence anomalies are divided into volcanic eruption magma channels; if the longitudinal depth is below the target layer, dividing the target layer into volcano invasion rock-slurry channels;
example 5
Referring to fig. 1 to 5, a volcanic tunnel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on a high-resolution electrical profile;
d. through comparison, if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and correspond to the magnetic anomaly, the seismic data are divided into volcanic magma channels, and if the coherent anomaly and the high-resistance anomaly strips on the seismic data are consistent and no magnetic anomaly exists, the seismic data are divided into volcanic hydrothermal channels;
e. and the volcanic rock-slurry channels with abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic rock-slurry channels, and the rest are crack type volcanic rock-slurry channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 25Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In step b, the depth range refers to a range of the threshold value being 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
In the step c, the vertical identification means that residual magnetic anomaly plane distribution is superposed according to a high-resistance piercing threshold value of 10000 omega.m, and the volcanic channel response is determined when the high-resistance strip is transversely larger than 1km and magnetic anomalies exist.
In the step d, the comparison specifically means that various earthquake coherence anomalies within an effective depth range are retained and compared with the electrical profile, and if magnetic anomalies exist and the width is greater than 1km and the seismic coherence anomalies on the seismic profile correspond to, the seismic coherence anomalies are divided into volcanic eruption magma channels; if the longitudinal depth is below the target layer, dividing the target layer into volcano invasion rock-slurry channels;
in the step e, the rest of the gap type volcanic rock slurry channels are strips, wherein the volcanic rock slurry channels with abnormal magnetic property, electric property and earthquake in the threshold value are superposed with gravity abnormity, the coherence abnormity is superposed with the substrate fracture, and the plane area is less than or equal to ten square kilometers.
Example 6
Referring to fig. 1 to 5, a volcanic tunnel identification and division method includes the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of seismic data of fidelity amplitude, a coherence attribute body in the depth range from a base to a target layer is calculated, and columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body are identified;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on a high-resolution electrical profile;
d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
In the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 35Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
In step b, the depth range refers to a range of the threshold value being 2500-3500 ms.
In the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
In the step c, the vertical identification means that residual magnetic anomaly plane distribution is superposed according to a high-resistance piercing threshold value of 11000 omega-m, and the volcanic channel response is determined when the high-resistance strip is larger than 1km transversely and has magnetic anomaly.
In the step d, the comparison specifically means that various earthquake coherence anomalies within an effective depth range are retained and compared with the electrical profile, and if magnetic anomalies exist and the width is greater than 1km and the seismic coherence anomalies on the seismic profile correspond to, the seismic coherence anomalies are divided into volcanic eruption magma channels; if the longitudinal depth is below the target layer, dividing the target layer into volcano intrusion rock-magma channels;
in the step e, the rest of the gap type volcanic rock slurry channels are strips, wherein the volcanic rock slurry channels with abnormal magnetic property, electric property and earthquake in the threshold value are superposed with gravity abnormity, the coherence abnormity is superposed with the substrate fracture, and the plane area is less than or equal to ten square kilometers.
The crack type volcanic magma channel comprises a linear volcanic magma channel, a point volcanic magma channel and a bead-like volcanic magma channel.
The method is particularly suitable for identifying the two-fold volcanic rock under the Sichuan basin, can identify that volcanic rock channels in the West Sichuan area are mainly in a slit type bead shape, and parts of treetop-shaped and punctiform reservoirs with high-quality outburst phases are spread in the south-south direction and the north-north direction of the slit, can well identify the plane distribution range of the volcanic channels, and can accurately depict the longitudinal spreading characteristics of the volcanic channels.
In FIG. 5, A (x) i ,y i T) is the desired point, B (x) i ,y i+1 ,t)、C(x i+1 ,y i T) is the adjacent seismic trace of point A in the horizontal and vertical measuring line directions respectively, and A (x) i ,y i T) defining a time window length, tau, centered on the seismic trace in which the point is located x 、τ y Respectively are two paths of B and Cthe coherent delay of time t with track A.
Thus, the C1 correlation coefficient is defined as:
where ρ is x (t,τ x ,x i ,y i ) Is the cross-correlation coefficient of A and B, ρ y (t,τ y ,x i ,y i ) Is the cross-correlation coefficient of A and C; the C1 correlation coefficient is the product of the maximum of the normalized correlation coefficient for a number of target tracks in the x-direction and a number of target tracks in the y-direction, at which time τ is the product of the maximum of the normalized correlation coefficient for high quality data x 、τ y Approximately equal to the apparent time tilt of the x and y axes.
Claims (8)
1. A volcanic channel identification and division method is characterized by comprising the following steps:
a. performing fidelity and amplitude-preserving seismic processing on the seismic data based on the single-point cutting, and processing the high-resolution electrical data;
b. on the basis of the seismic data of fidelity and amplitude preservation, calculating a coherent attribute body in the depth range from a base to a target layer and identifying the columnar, chimney-shaped and linear coherent abnormality in the coherent attribute body;
c. vertically identifying high-resistance abnormal strips from a substrate to a target layer on the high-resolution electrical profile;
d. through comparison, if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and corresponds to the magnetic abnormality, the seismic data is divided into a volcanic magma channel, and if the coherent abnormality on the seismic data is consistent with the high-resistance abnormal strip and does not have the magnetic abnormality, the seismic data is divided into a volcanic hydrothermal channel;
e. and the volcanic magma channels with the abnormal magnetic property, electric property and earthquake in the threshold value are superposed with the gravity abnormity, the strip with the coherence abnormity superposed with the substrate fracture and the plane area larger than ten square kilometers is divided into penetration type volcanic magma channels, and the rest are crack type volcanic magma channels.
2. The volcanic channel identification and division method according to claim 1, wherein the volcanic channel identification and division method comprises the following steps: in the step a, the step of performing fidelity amplitude-preserving seismic processing specifically refers to processing seismic data into seismic data with 30Hz main frequency and 6-35Hz bandwidth, and the step of processing high-resolution electrical data specifically refers to processing the high-resolution electrical data into electrical data with transverse identification power reaching more than 1 km.
3. The volcanic channel identification and division method according to claim 1, wherein the volcanic channel identification and division method comprises the following steps: in the step b, the depth range refers to a range with a threshold value of 2500-3500 ms.
4. The volcanic channel identification and division method according to claim 1, characterized in that: in the step b, identifying the columnar, chimney-shaped and linear coherence anomalies in the coherence attribute body specifically means removing the anomaly coherence reaching a shallower depth by comparing with the depth range.
5. The volcanic channel identification and division method according to claim 1, wherein the volcanic channel identification and division method comprises the following steps: in the step c, the vertical identification means that residual magnetic anomaly plane distribution is superposed according to a high-resistance piercing threshold value of 8000-11000 omega-m, and the volcanic channel response is determined when the high-resistance strip is transversely larger than 1km and magnetic anomalies exist.
6. The volcanic channel identification and division method according to claim 1, characterized in that: in the step d, comparing specifically means comparing various earthquake coherence anomalies in the reserved effective depth range with the electrical profile, and if magnetic anomalies exist and the width is greater than 1km and corresponds to a strong coherence anomaly on the earthquake profile, dividing the earthquake coherence anomalies into volcanic explosion magma channels; if the longitudinal depth is below the target layer, the tunnel is divided into volcanic invasion rock-slurry channels.
7. The volcanic channel identification and division method according to claim 1, wherein the volcanic channel identification and division method comprises the following steps: in the step e, the rest of the gap type volcanic rock slurry channels are strips, wherein the volcanic rock slurry channels with abnormal magnetic property, electric property and earthquake in the threshold value are superposed with gravity abnormity, the coherence abnormity is superposed with the substrate fracture, and the plane area is less than or equal to ten square kilometers.
8. The volcanic channel identification and division method according to claim 1, characterized in that: the crack type volcanic rock-slurry channel comprises a linear volcanic rock-slurry channel, a point volcanic rock-slurry channel and a beaded volcanic rock-slurry channel.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110496212.3A CN113281810B (en) | 2021-05-07 | 2021-05-07 | Volcanic channel identification and division method |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110496212.3A CN113281810B (en) | 2021-05-07 | 2021-05-07 | Volcanic channel identification and division method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113281810A CN113281810A (en) | 2021-08-20 |
| CN113281810B true CN113281810B (en) | 2022-11-11 |
Family
ID=77278108
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110496212.3A Active CN113281810B (en) | 2021-05-07 | 2021-05-07 | Volcanic channel identification and division method |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113281810B (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104007468A (en) * | 2014-05-23 | 2014-08-27 | 中国地质大学（武汉） | Method for depicting volcanic space distribution based on amplitude-variance cube seismic attributes |
| CN112034515A (en) * | 2020-08-14 | 2020-12-04 | 中国海洋石油集团有限公司 | Volcanic channel identification method based on unsupervised neural network |
| CN112114358A (en) * | 2019-06-20 | 2020-12-22 | 中国石油化工股份有限公司 | Underground volcanic channel identification method based on three-dimensional seismic data representation |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US8861309B2 (en) * | 2011-01-31 | 2014-10-14 | Chevron U.S.A. Inc. | Exploitation of self-consistency and differences between volume images and interpreted spatial/volumetric context |
-
2021
- 2021-05-07 CN CN202110496212.3A patent/CN113281810B/en active Active
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104007468A (en) * | 2014-05-23 | 2014-08-27 | 中国地质大学（武汉） | Method for depicting volcanic space distribution based on amplitude-variance cube seismic attributes |
| CN112114358A (en) * | 2019-06-20 | 2020-12-22 | 中国石油化工股份有限公司 | Underground volcanic channel identification method based on three-dimensional seismic data representation |
| CN112034515A (en) * | 2020-08-14 | 2020-12-04 | 中国海洋石油集团有限公司 | Volcanic channel identification method based on unsupervised neural network |
Non-Patent Citations (4)
| Title |
|---|
| 大丰兴化地区火山岩地震识别方法研究;左国平 等;《石油物探》;20110525;第50卷(第03期);252-259 * |
| 渤海海域沙东地区馆陶组火山通道相构型及成因机制;孙希家 等;《特种油气藏》;20170228;第24卷(第01期);53-57 * |
| 火山发育区通道相类型、特征、成因及对油气的控制作用;孙希家 等;《地质论评》;20180715;第64卷(第04期);937-945 * |
| 用重磁电异常信息模式识别石炭系火山岩岩性――以准噶尔盆地陆东地区为例;索孝东 等;《新疆石油地质》;20110601;第32卷(第03期);318-320 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113281810A (en) | 2021-08-20 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11193372B2 (en) | Oil and gas zone effectiveness evaluation method and apparatus | |
| US10203427B2 (en) | Correlation techniques for passive electroseismic and seismoelectric surveying | |
| EP2962135B1 (en) | System and method for detecting a fracture in a rock formation using an electromagnetic source | |
| Khesin et al. | Interpretation of geophysical fields in complicated environments | |
| West et al. | Interactive seismic facies classification using textural attributes and neural networks | |
| GB2527239B (en) | Sensors for passive electroseismic and seismoelectric surveying | |
| US10810331B2 (en) | System for predicting induced seismicity potential resulting from injection of fluids in naturally fractured reservoirs | |
| Sullivan et al. | An integrated approach to characterization and modeling of deep-water reservoirs, Diana field, western Gulf of Mexico | |
| CN109709626B (en) | Structural closed type weak open layered rock heat-storage geothermal field prospecting method | |
| CN111696208B (en) | Geological-geophysical three-dimensional modeling method based on multi-data fusion | |
| CN113281810B (en) | Volcanic channel identification and division method | |
| Zhenwu et al. | Practices and expectation of high-density seismic exploration technology in CNPC | |
| CN114740538A (en) | Skarn type iron-rich ore deep exploration method and system based on multivariate geophysical | |
| CN104330824A (en) | Oil layer identification method by utilizing energy relative change rate | |
| CN109425890B (en) | Carbonate karst cave reservoir development scale earthquake identification method and system | |
| Pearce et al. | Interaction between hydrothermal fluids and fault systems in the in the Southern Andes revealed by magnetotelluric and seismic data | |
| CN111965724B (en) | Stratum fracture-cavity type identification method and device | |
| EP3943985A1 (en) | System and method for mapping and monitoring reservoirs by electromagnetic crosswell and optimizing production | |
| Xiao et al. | Application of scattering image wavelet transform in cave recognition: A case study on a bedrock buried hill reservoir in Bongor Basin, Chad | |
| Kusumah et al. | Horizontal derivative from gravity data as a tool for drilling target guide in Wayang Windu Geothermal Field, Indonesia | |
| Wang et al. | 3D real-time imaging for electromagnetic fracturing monitoring based on deep learning | |
| Qiu et al. | Application of stochastic seismic invertion of simulated anneal method in reservoir prediction of Xing'anling Group in Bei14 Block | |
| CN117270037A (en) | Oil-gas migration prediction method based on seepage structure fine description | |
| CN115016017A (en) | Shale bed series identification method and device | |
| CN117950021A (en) | Ultra-deep carbonate broken solution reservoir earthquake phase control inversion method and device |
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 |