CN111257926B - Method for predicting ancient valley uranium reservoir by using old seismic data - Google Patents
Method for predicting ancient valley uranium reservoir by using old seismic data Download PDFInfo
- Publication number
- CN111257926B CN111257926B CN201811465517.2A CN201811465517A CN111257926B CN 111257926 B CN111257926 B CN 111257926B CN 201811465517 A CN201811465517 A CN 201811465517A CN 111257926 B CN111257926 B CN 111257926B
- Authority
- CN
- China
- Prior art keywords
- uranium
- seismic
- valley
- reservoir
- ancient
- 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
- 229910052770 Uranium Inorganic materials 0.000 title claims abstract description 140
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 title claims abstract description 139
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000005553 drilling Methods 0.000 claims abstract description 25
- 239000011435 rock Substances 0.000 claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 16
- 239000003208 petroleum Substances 0.000 claims abstract description 13
- 238000012795 verification Methods 0.000 claims abstract description 7
- 238000004321 preservation Methods 0.000 claims abstract description 6
- 239000004576 sand Substances 0.000 claims description 32
- 230000003647 oxidation Effects 0.000 claims description 27
- 238000007254 oxidation reaction Methods 0.000 claims description 27
- 239000010410 layer Substances 0.000 claims description 23
- 230000033558 biomineral tissue development Effects 0.000 claims description 19
- 238000011161 development Methods 0.000 claims description 18
- 230000018109 developmental process Effects 0.000 claims description 18
- 238000009826 distribution Methods 0.000 claims description 16
- 239000011229 interlayer Substances 0.000 claims description 13
- 230000002349 favourable effect Effects 0.000 claims description 9
- 208000035126 Facies Diseases 0.000 claims description 8
- 238000003384 imaging method Methods 0.000 claims description 7
- 238000001914 filtration Methods 0.000 claims description 6
- 238000010276 construction Methods 0.000 claims description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- 238000007596 consolidation process Methods 0.000 claims description 4
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- NIFIFKQPDTWWGU-UHFFFAOYSA-N pyrite Chemical compound [Fe+2].[S-][S-] NIFIFKQPDTWWGU-UHFFFAOYSA-N 0.000 claims description 4
- 229910052683 pyrite Inorganic materials 0.000 claims description 4
- 239000011028 pyrite Substances 0.000 claims description 4
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 239000000428 dust Substances 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 230000009467 reduction Effects 0.000 claims description 3
- 230000001629 suppression Effects 0.000 claims description 3
- 238000010183 spectrum analysis Methods 0.000 abstract description 3
- 230000015572 biosynthetic process Effects 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 description 3
- 230000008021 deposition Effects 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 229910052760 oxygen Inorganic materials 0.000 description 3
- 239000001301 oxygen Substances 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 238000009933 burial Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 235000017399 Caesalpinia tinctoria Nutrition 0.000 description 1
- 206010054949 Metaplasia Diseases 0.000 description 1
- 241000388430 Tara Species 0.000 description 1
- 230000009471 action Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 239000007789 gas Substances 0.000 description 1
- 230000008595 infiltration Effects 0.000 description 1
- 238000001764 infiltration Methods 0.000 description 1
- 239000012633 leachable Substances 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 230000015689 metaplastic ossification Effects 0.000 description 1
- 238000005192 partition Methods 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 239000013049 sediment Substances 0.000 description 1
- 238000004062 sedimentation Methods 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 238000003892 spreading Methods 0.000 description 1
- 238000003860 storage Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000013519 translation 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/282—Application of seismic models, synthetic seismograms
-
- 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
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/364—Seismic filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/16—Survey configurations
- G01V2210/163—Cross-well
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/30—Noise handling
- G01V2210/32—Noise reduction
- G01V2210/324—Filtering
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/622—Velocity, density or impedance
- G01V2210/6226—Impedance
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention belongs to the technical field of uranium ore exploration, and particularly relates to a method for predicting an ancient valley uranium reservoir by using old seismic data. At present, drilling is used as a main ore exploration means and method, but the investment cost is high, the period is long, and the prediction of the ancient valley uranium reservoir by the traditional method is difficult. The invention mainly comprises the following steps: the method comprises the following steps: determining the uranium reservoir attribute by using known geological drilling data of uranium ores through a structural criterion, a uranium source criterion, a lithology-lithofacies criterion and a rock geochemistry criterion; step two: performing fidelity and amplitude preservation processing on forehand seismic data by combining a petroleum two-dimensional seismic profile, highlighting shallow seismic reflection information, analyzing ancient valley seismic phases of the seismic profile of a known uranium mine hole, performing wave impedance inversion according to known well wave impedance and seismic wave spectrum analysis, realizing shallow interpretation and structure recognition of the petroleum seismic profile, and determining the uranium reservoir attribute; step three: and (5) performing combination verification. The method can be used for predicting the uranium reservoir in the valley of the ancient river.
Description
Technical Field
The invention belongs to the technical field of uranium ore exploration, and particularly relates to a method for predicting an ancient valley uranium reservoir by using old seismic data.
Background
The sedimentary basin is often a 'treasure collecting basin' rich in various energy resources, not only contains abundant coal and oil gas resources, but also has a large amount of uranium ore resources. The northern regional leachable sandstone type uranium ore is the key point of the uranium ore prospecting work in China, the uranium ore mainly takes a deeply buried blind ore body and basically has no outcrop, the main prospecting means and method is drilling at present, but the investment cost is high, the period is long, and the prediction of the uranium reservoir in the valley of the ancient river by the traditional method is difficult.
The method utilizes reflection characteristics on old seismic sections, wave impedance and seismic wave spectrum analysis and well-seismic combination to carry out wave impedance inversion in a large area without cores or outcrops. The method has the advantages of analyzing basin evolution development and locking a target layer, can accurately determine the distribution range of the uranium reservoir, and the change characteristics of boundary lines, burial depth, lithology and the like of the uranium reservoir, and comprehensively analyzes uranium mineralization conditions by combining the characteristics of reflecting the metaplasia and uranium mineralization of the drilled uranium ore to predict the uranium reservoir in the valley of the ancient river.
Disclosure of Invention
The invention provides a method for predicting an ancient valley uranium reservoir by using old seismic data, which can overcome the defects (only using uranium ore geological drilling data) of the prior art, correct a two-dimensional petroleum seismic profile by using known uranium ore geological drilling, greatly improve the working efficiency of uranium ore prospecting and improve the partition precision of uranium mineralization scenic regions. The prior petroleum seismic section is heavily focused on deep interpretation, has large detection depth, obviously shows the fold structure characteristics of the underlayer of the Sehan group and the base characteristics of a research area on the section, establishes an isochronous framework for sedimentary facies research to search reservoirs and sand bodies which are possibly mineralized, reduces the target range, but lacks intuitive core data; drilling holes of uranium ores are mostly single points, the hole depth is shallow, the quantity is small, the holes are not uniformly distributed, uranium ore geological data can directly display uranium mineralization information and a uranium mineralization geological environment, and the exploration of unknown areas is insufficient.
The technical scheme for realizing the invention is as follows:
a method for predicting an ancient valley uranium reservoir by using old seismic data mainly comprises the following steps: the method comprises the following steps: determining the uranium reservoir attribute by using known geological drilling data of uranium ores through a structural criterion, a uranium source criterion, a lithology-lithofacies criterion and a rock geochemistry criterion;
step two: performing fidelity and amplitude preservation processing on forehand seismic data by combining a petroleum two-dimensional seismic profile, highlighting shallow seismic reflection information, analyzing ancient valley seismic phases of the seismic profile of a known uranium mine hole, performing wave impedance inversion according to known well wave impedance and seismic wave spectrum analysis, realizing shallow interpretation and structure recognition of the petroleum seismic profile, and determining the uranium reservoir attribute;
step three: performing combined verification on the uranium reservoir attribute obtained in the step one and the uranium reservoir attribute obtained in the step two, and applying the combined verification to prediction of a uranium reservoir;
the third step comprises the following specific steps: forecasting and defining a uranium mineralization favorable zone in a uranium mine low working degree area by combining an earthquake phase according to an ancient valley of known well region earthquake inversion of drilling river channel sand; firstly, finding out wave impedance data of a drilled ancient river channel well, analyzing the wave impedance data in a crossed manner, determining a threshold range of the wave impedance data, extracting wave impedance data of a development layer section corresponding to the ancient river channel in an inversion wave impedance body according to the determined threshold, and finally determining an ancient river valley distribution range through plane or three-dimensional display; and (3) obtaining a time-depth relation pair after the synthetic record is calibrated, calibrating the specific time depth of the uranium-bearing ancient river sand body drilled on the time section, analyzing the corresponding seismic waveform, and searching the seismic distribution range with similar waveform characteristics by utilizing attributes such as seismic clustering analysis and the like, namely the possible long-range target area.
The construction criterion is as follows: the most favorable structure of sandstone-type uranium mineralization is on a large stable structural slope belt, the sedimentary stratum has certain lift and valley undercut forms, and the later structural evolution has inheritance to the sedimentary period.
The uranium source criterion is as follows: and the material source region rock stratum/body uranium is high in abundance and strong in quasi-plain-translation, or the target layer and the upper and lower layers of the target layer are uranium-rich layers, and uranium is obviously migrated.
The lithology-lithofacies criterion is as follows: the sand bodies such as braided rivers, sectorial delta and the like with large development scale and good connectivity have stable 'mud-sand-mud' reservoir structures, the sand bodies are rich in reducing media such as carbon dust, pyrite and the like, and the consolidation degree of rocks is relatively low.
The geochemical rule of the rock is as follows: the geochemical zonation of the rock from oxidation to reduction exists, the oxidation types are mainly interlayer oxidation and diving-interlayer oxidation, the development of the interlayer oxidation zone and the diving-interlayer oxidation zone have certain scales, the distribution boundary of an oxidation zone front line or gray (residual) sand body can be presumed on a plane, and the diagenesis is enriched near the oxidation zone front line or in the gray residual sand body; according to the known well oxide sand body development condition, lithology-lithofacies analysis is combined, the oxide sand body distribution direction and scale are identified, the development position and scale of the oxide zone front edge are predicted, and the ancient valley sandstone type uranium mineralization scenic region is defined.
The seismic data fidelity and amplitude-preserving processing comprises the steps of firstly carrying out frequency domain filtering, setting a band-pass filtering window, and suppressing low-frequency surface waves and other interference signals through effective signals; meanwhile, amplitude preservation processing is carried out on the seismic data through noise suppression, amplitude compensation, wavelet processing and migration imaging means.
The method comprises the steps of analyzing ancient valley seismic facies of seismic sections of known uranium mine holes, drilling holes in the uranium mines to obtain the depth of a high-gamma sandstone section, then extracting well side-channel wavelets, calibrating the high-gamma sandstone in detail, carrying out speed analysis and manufacturing a time-depth conversion ruler.
The wave impedance inversion according to the known well wave impedance and seismic wave-front analysis means that firstly, the speed of the underground rock stratum is obtained by utilizing the acoustic wave time difference data of the known well drilling, the speed is multiplied by the density data to obtain well wave impedance data AI, and the reflection coefficient R (t) is extracted according to the well wave impedance data. The formula is as follows:
in the formula, R (t) i ) Reflection coefficients for the ith and (i + 1) th reflection layers on the well, AI (t) i ) The wave impedance of the ith reflection layer on the well, and t is the time depth; then extracting seismic wavelets through well-side seismic data channels, and establishing a geological model through seismic structure interpretation; and introducing a drilled well combination model for interpolation to obtain reflection coefficients of all layers, and deconvoluting the reflection coefficients and wavelets to obtain a wave impedance inversion volume. And solving reservoir parameters of the unknown area according to the correlation between the known well geological information and the wave impedance inversion body.
The invention has the following effects: the invention provides a new uranium ore prospecting method, aiming at a low working degree area of a uranium ore, firstly, by means of an old two-dimensional petroleum seismic profile, shallow seismic data secondary processing is carried out, stratum identification and comparison are carried out by combining known geological drilling holes of the uranium ore, an isochronous stratum trellis is built, the distribution range and the burial depth of a uranium reservoir in an unknown area are presumed, the uranium ore prospecting working efficiency is greatly improved, drilling verification can be directly and accurately carried out, the investment drilling cost is saved, the blind ore prospecting precision is improved, and a level I distant area favorable for uranium ore formation is particularly trapped in a certain depression of a two-link basin.
Drawings
FIG. 1 is a flow chart of an embodiment;
FIG. 2 is a comparison of new and old results of shallow interpretation of two-dimensional seismic profiles in a certain work area;
FIG. 3 is an X-well synthetic record calibration;
FIG. 4 is a top flat bottom convex valley seismic facies identified by well seismic combination in a certain work area;
FIG. 5 is a schematic geological profile of a work area;
FIG. 6 is a diagram of a historic valley seismic section predicted by combining seismic phase and wave impedance inversion in a work area;
in the figure: 1-ancient line + new line; 2-Saohana group upper segment; 3-Saohan tara group lower segment; 4-yellow oxidized sandstone; 5-grey reduced sandstone; 6-angle non-integral boundary line; 7-the rock section boundary line; 8-front edge of oxidation zone; 9-surface boundary line; 10-sandstone type industrial uranium ore body; 11-uranium ore body
Detailed Description
The method for predicting the ancient valley type uranium reservoir by using old seismic data according to the present invention is further described with reference to the accompanying drawings and specific embodiments, and the flow of the specific embodiment is as shown in fig. 1.
1. Determining the uranium reservoir attribute by using known geological drilling data of uranium ores and through technical methods such as a structural criterion, a uranium source criterion, a lithology-lithofacies criterion, a rock geochemistry criterion and the like;
2. and combining a petroleum two-dimensional seismic profile, performing fidelity amplitude-preserving processing on foreigner seismic data, highlighting shallow seismic reflection information, analyzing the ancient valley seismic phase of the seismic profile of the known uranium mine hole, performing wave impedance inversion according to known well wave impedance and seismic wave-front analysis, realizing shallow interpretation and structure recognition of the petroleum seismic profile, and determining the uranium reservoir property. Finally, combining and verifying the two methods to realize a new technical method for accurately predicting the sandstone-type uranium reservoir in the ancient river valley.
The uranium reservoir is as follows: the space in the sedimentation basin for uranium mineralization fluid transportation and uranium ore storage is the most basic condition for uranium ore formation.
The ancient valley type sandstone uranium reservoir is as follows: the ancient river valley refers to a strip valley formed by the tectonic action in a sedimentary basin, and is filled with alluvial substance and lake sediments in the later period, and a uranium deposit produced in sandstone of the ancient river valley is called as an ancient river valley sandstone uranium deposit.
The construction criterion is as follows: the most favorable structure of the sandstone-type uranium mineralization is on a large stable structural slope zone, a stratum after deposition is inclined to a certain degree and has a valley undercut form, later structural evolution has inheritance to the deposition period, favorable structural conditions are one of necessary conditions of the sandstone-type uranium mineralization, and the structure is favorable for long-term infiltration of oxygen-containing water containing uranium, so that the ancient valley sandstone-type uranium deposit is formed. Shallow interpretation of seismic sections enables accurate identification of the developmental conditions and extent of such formations.
The uranium source criterion is as follows: the physical source region rock stratum (body) has high uranium abundance and strong quasi-plain primalization, or the target layer and the upper and lower layers thereof are uranium-rich layers, uranium is obviously migrated, the standard is one of basic conditions for forming the valley type uranium deposit, and abundant uranium sources can permeate into a deposition basin along a construction slope along with oxygen-containing water containing uranium, and abundant uranium sources can be provided for the valley sandstone type uranium reservoir.
The lithology-lithofacies criterion is as follows: the sand bodies such as braided rivers, sectorial delta and the like with large development scale and good connectivity have stable 'mud-sand-mud' reservoir structures, are rich in reducing media such as carbon dust, pyrite and the like, and have relatively low rock consolidation degree. The lithology-lithofacies development condition controls the development of the ancient valley sandstone reservoir, the ancient valley sandstone reservoir is a migration channel of uranium-containing oxygen-containing water carrying a uranium source, and the spreading form and scale of sand bodies can be predicted according to a shallow seismic profile and known well energy.
The geochemical rule of the rock is as follows: the geochemical zonation of the rock from oxidation to reduction exists, the oxidation types are mainly interlayer oxidation and diving-interlayer oxidation, the development of the interlayer oxidation zone and the diving-interlayer oxidation zone have certain scales, the distribution boundary of an oxidation zone front line or gray (residual) sand body can be presumed on a plane, and the mineralization is enriched near the oxidation zone front line or in the gray residual sand body. According to the known well oxide sand body development condition, lithology-lithofacies analysis is combined, the oxide sand body distribution direction and scale are identified, the development position and scale of the oxide zone front edge are predicted, and the ancient valley sandstone type uranium mineralization scenic region is defined.
The seismic data fidelity and amplitude-preserving processing comprises the steps of firstly carrying out frequency domain filtering, setting a band-pass filtering window, and suppressing low-frequency surface waves and other interference signals through effective signals; meanwhile, amplitude preservation processing is carried out on the seismic data through noise suppression, amplitude compensation, wavelet processing, migration imaging and other means. Highlighting shallow seismic reflection information through seismic data fidelity amplitude preservation processing; the data resolution and the imaging effect are improved, and the imaging precision is improved; the forecasting of the sandstone-type uranium reservoir in the valley and the identification of a fracture structure can be met;
and analyzing the ancient valley seismic facies of the seismic section of the known uranium mine hole. Drilling holes in a uranium mine to obtain the depth of a high-gamma sandstone section, then extracting well bypass wavelets, calibrating the high-gamma sandstone in detail, carrying out speed analysis, and manufacturing a time-depth conversion ruler;
and carrying out wave impedance inversion according to the known well wave impedance and seismic wave-front analysis. Firstly, obtaining the underground rock formation speed by using the known well drilling sound wave time difference data, multiplying the underground rock formation speed by the density data to obtain well wave impedance data AI, and extracting a reflection coefficient R (t) according to the well wave impedance data. The formula is as follows:
in the formula, R (t) i ) Reflection coefficients for the ith and (i + 1) th reflection layers on the well, AI (t) i ) Is the i-th well
The wave impedance of each reflective layer, t, is the time depth.
Then, seismic wavelets are extracted through the well-side seismic data traces, and a geological model is built through seismic tectonic interpretation. And introducing a drilled well combination model for interpolation to obtain reflection coefficients of all layers, and deconvoluting the reflection coefficients and wavelets to obtain a wave impedance inversion volume. And solving reservoir parameters of the unknown area according to the correlation between the known well geological information and the wave impedance inversion body.
And obtaining a data volume capable of reflecting the plane distribution range and the space distribution characteristics of the sand body according to the inversion.
The two methods are combined for verification, namely, according to the ancient valley of the well region seismic inversion known in the sand body of the drilling river channel, the beneficial uranium mineralization zone is predicted and defined in the low working degree area of uranium ore by combining the seismic phase. The method comprises the steps of firstly finding out wave impedance data of a well drilled with an ancient river channel, analyzing the wave impedance data through intersection, determining a threshold range of the wave impedance data, extracting wave impedance data of a development interval corresponding to the ancient river channel in an inversion wave impedance body according to the determined threshold, and finally determining the distribution range of the ancient river valley through plane or three-dimensional display. And (3) obtaining a time-depth relation pair after the synthetic record is calibrated, calibrating the specific time depth of the uranium-bearing ancient river sand body drilled on the time section, analyzing the corresponding seismic waveform, and searching the seismic distribution range with similar waveform characteristics by utilizing attributes such as seismic clustering analysis and the like, namely the possible long-range target area.
Fig. 2 identifies a large-scale and good-connectivity river sand body through drilling of a material uranium ore and combining a construction criterion, a uranium source criterion, a lithology-lithology criterion and a rock geochemistry criterion, wherein the river sand body has a stable mud-sand-mud structure, the sand body is rich in carbon chips, pyrite and other reducing media, and the rock consolidation degree is relatively low; the uranium reservoir has characteristics of a diving-interlayer oxidation zone, yellow is an after-generated oxidized sandstone, gray is a primary reduced sandstone, a uranium ore body is produced in the gray sandstone at the lower part of a front line of the oxidation zone, and the attribute of the uranium reservoir is finally determined.
Through fidelity and amplitude-preserving processing of a petroleum two-dimensional seismic section, shallow seismic reflection information is highlighted, data resolution and imaging effect are improved by combining uranium ore geological drilling, imaging precision is improved, and the newly processed seismic section shallow river valley form is more obvious and has a uranium reservoir structure of mud-sand-mud as shown in comparison in FIG. 3;
as shown in FIG. 4, the research on the correlation between synthetic records and seismic phases shows that the correlation is good in a yellow region, a high-gamma sandstone section has a longitudinal wave impedance box shape and is low in numerical value, the uranium content of a target layer sand body is high, the natural gamma value is high and shows a plurality of peak values, a resistivity curve is large in a box shape, and the characteristic of a river channel sand body is shown overall.
In the figure 5, a No. 1 uranium ore geological drilling hole is cast on a corresponding petroleum earthquake section, synthetic record calibration is carried out, an ancient valley type tectonic seismic facies typical of a exploration target stratum top flat bottom convex is identified through well-seismic combination, and the seismic facies is a potential uranium reservoir and comprehensively determines the uranium reservoir attribute.
And FIG. 6, wave impedance inversion is carried out according to drilling and wave impedance and seismic wave-front analysis of the known uranium ore, and the known ancient valley uranium ore deposit is obviously identified to be located in an ancient valley seismic phase region of seismic inversion, so that the uranium reservoir attributes are verified mutually, and a distant scenic region is predicted and defined in combination with the seismic phase in the north (low-working-degree region) of the known uranium ore deposit.
As shown in figures 2-6, the invention is applied to a level I distant view area with a depression in a Liqun basin and favorable uranium mineralization.
Claims (8)
1. A method for predicting an ancient valley uranium reservoir by using old seismic data is characterized by comprising the following steps: the method comprises the following steps: determining the uranium reservoir attribute by using known geological drilling data of uranium ores through a structural criterion, a uranium source criterion, a lithology-lithofacies criterion and a rock geochemistry criterion;
step two: combining a petroleum two-dimensional seismic profile, performing fidelity amplitude-preserving processing on foreigner seismic data, highlighting shallow seismic reflection information, analyzing an ancient valley seismic phase of the seismic profile of a known uranium mine hole, performing wave impedance inversion according to known well wave impedance and seismic wave-front analysis, realizing shallow interpretation and structure recognition of the petroleum seismic profile, and determining the uranium reservoir attribute;
step three: performing combined verification on the uranium reservoir attribute obtained in the step one and the uranium reservoir attribute obtained in the step two, and applying the combined verification to prediction of a uranium reservoir;
the third step comprises the following specific steps: predicting and defining a uranium mineralization favorable zone by combining seismic facies in a low working degree area of uranium ores according to an ancient valley of known well region seismic inversion of drilling river channel sand; firstly, finding out wave impedance data of a well drilled with an ancient river channel, analyzing the wave impedance data through intersection, determining a threshold range of the wave impedance data, extracting wave impedance data of a development interval corresponding to the ancient river channel in an inversion wave impedance body according to the determined threshold, and finally determining an ancient river valley distribution range through plane or three-dimensional display; and (3) obtaining a time-depth relation pair after the synthetic record is calibrated, calibrating the specific time depth of the uranium-bearing ancient river sand body drilled on a time section, analyzing the corresponding seismic waveform, and searching the seismic distribution range with similar waveform characteristics by using the seismic clustering analysis attribute, namely the possible distant view target area.
2. The method for predicting the uranium reservoir in the ancient valley by using the old seismic data as set forth in claim 1, wherein: the construction criterion is as follows: the most favorable structure of sandstone-type uranium mineralization is on a large stable structural slope belt, the sedimentary stratum has certain lift and has a valley undercut form, and the later structural evolution has inheritance to the sedimentary period.
3. The method for predicting the uranium reservoir in the valley of the ancient river according to the old seismic data as claimed in claim 1, wherein: the uranium source criterion is as follows: the abundance of the uranium in the rock strata/body of the source region is high, the quasi-plain is strong, or the target layer and the upper and lower layers of the target layer are uranium-rich layers, and the uranium emigration is obvious.
4. The method for predicting the uranium reservoir in the valley of the ancient river according to the old seismic data as claimed in claim 1, wherein: the lithology-lithofacies criterion is as follows: the braided river sand body and the fan delta sand body which have large development scale and good connectivity have a stable mud-sand-mud reservoir structure, the sand body is rich in carbon dust and pyrite reducing medium, and the consolidation degree of the rock is relatively low.
5. The method for predicting the uranium reservoir in the valley of the ancient river according to the old seismic data as claimed in claim 1, wherein: the geochemical rule of the rock is as follows: the geochemical zonation of the rock from oxidation to reduction exists, the oxidation types are mainly interlayer oxidation and diving-interlayer oxidation, the development of the interlayer oxidation zone and the diving-interlayer oxidation zone have certain scales, the distribution boundary of an oxidation zone front line or gray sand body is presumed on a plane, and the mineralization is enriched near the oxidation zone front line or in the gray residual sand body; according to the development condition of the oxide sand bodies in the known well, lithologic-lithofacies analysis is combined, the distribution direction and scale of the oxide sand bodies are identified, the development position and scale of the front edge of the oxide zone are predicted, and the gulf valley sandstone type uranium mineralization distant scenic spot is defined.
6. The method for predicting the uranium reservoir in the valley of the ancient river according to the old seismic data as claimed in claim 1, wherein: the fidelity and amplitude-preserving processing of the seismic data comprises the steps of firstly carrying out frequency domain filtering, setting a band-pass filtering window, and suppressing low-frequency surface waves and other interference signals through effective signals; meanwhile, amplitude preservation processing is carried out on the seismic data through noise suppression, amplitude compensation, wavelet processing and migration imaging means.
7. The method for predicting the uranium reservoir in the ancient valley by using the old seismic data as set forth in claim 1, wherein: the method comprises the steps of analyzing ancient valley seismic facies of seismic sections of known uranium mine holes, drilling holes in the uranium mines to obtain the depth of a high-gamma sandstone section, then extracting well side-channel wavelets, calibrating the high-gamma sandstone in detail, carrying out speed analysis and manufacturing a time-depth conversion ruler.
8. The method for predicting the uranium reservoir in the valley of the ancient river according to the old seismic data as claimed in claim 1, wherein: the wave impedance inversion is carried out according to the known well wave impedance and seismic wave-front analysis, namely, firstly, the speed of an underground rock stratum is obtained by utilizing the known well drilling sound wave time difference data, the speed is multiplied by density data to obtain well wave impedance data AI, and a reflection coefficient R (t) is extracted according to the well wave impedance data; the formula is as follows:
in the formula, R (t) i ) Reflection coefficients for the ith and (i + 1) th reflection layers on the well, AI (t) i ) The wave impedance of the ith reflection layer on the well, and t is the time depth; then extracting seismic wavelets through well-side seismic data channels, and establishing a geological model through seismic structure interpretation; introducing a drilled well combination model for interpolation to obtain reflection coefficients of all layers, and performing deconvolution on the reflection coefficients and wavelets to obtain a wave impedance inversion body; and solving reservoir parameters of the unknown area according to the correlation between the known well geological information and the wave impedance inversion body.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811465517.2A CN111257926B (en) | 2018-12-03 | 2018-12-03 | Method for predicting ancient valley uranium reservoir by using old seismic data |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811465517.2A CN111257926B (en) | 2018-12-03 | 2018-12-03 | Method for predicting ancient valley uranium reservoir by using old seismic data |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN111257926A CN111257926A (en) | 2020-06-09 |
| CN111257926B true CN111257926B (en) | 2022-07-26 |
Family
ID=70946937
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811465517.2A Active CN111257926B (en) | 2018-12-03 | 2018-12-03 | Method for predicting ancient valley uranium reservoir by using old seismic data |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN111257926B (en) |
Families Citing this family (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN112764122B (en) * | 2019-10-21 | 2023-12-12 | 中国石油化工股份有限公司 | Method for confining uranium reservoir of ancient river channel sandstone |
| CN113806674B (en) * | 2020-06-11 | 2024-10-11 | 中国石油化工股份有限公司 | Quantification method and device for longitudinal scale of ancient river channel, electronic equipment and storage medium |
| CN111983721A (en) * | 2020-08-26 | 2020-11-24 | 核工业北京地质研究院 | Sandstone uranium ore mud-sand-mud geological structure identification method and system |
| CN111983689B (en) * | 2020-08-27 | 2023-06-02 | 青海省第四地质勘查院 | Near-source electromagnetic seismoelectric joint GR inversion method |
| CN114139328B (en) * | 2020-09-03 | 2022-11-04 | 中国石油天然气股份有限公司 | Prediction method for favorable ore-forming zone of sandstone-type uranium ore in interlayer oxidation zone |
| CN112580190A (en) * | 2020-11-20 | 2021-03-30 | 核工业二〇八大队 | Volcanic rock type uranium ore attack depth blind finding exploration method |
| CN113514886B (en) * | 2021-07-22 | 2021-12-10 | 核工业北京地质研究院 | Geological-seismic three-dimensional prediction method for beneficial part of sandstone-type uranium deposit mineralization |
| CN113933911B (en) * | 2021-08-30 | 2024-06-11 | 核工业北京地质研究院 | Method for distinguishing metaplasia sandstone from original sedimentary red sandstone |
| CN114236638B (en) * | 2021-11-04 | 2024-08-16 | 核工业航测遥感中心 | Quick paleo-river channel setting method based on magnetic weight information |
| CN114185105A (en) * | 2021-11-15 | 2022-03-15 | 核工业二〇八大队 | Method for delineating favorable area of sandstone-type uranium mineralization |
| CN114114458B (en) * | 2021-11-17 | 2024-01-12 | 核工业北京地质研究院 | Sandstone type uranium ore deep blind ore body prediction method under thick coverage overburden background |
| CN116381794B (en) * | 2023-05-31 | 2023-07-28 | 东北石油大学三亚海洋油气研究院 | Method for delineating sandstone-type uranium deposit beneficial zone by utilizing oil and gas transportation and aggregation rule |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104237964A (en) * | 2014-08-14 | 2014-12-24 | 核工业北京地质研究院 | Paleochannel type sandstone-type uranium ore deposit positioning method |
| CN106501871A (en) * | 2016-09-14 | 2017-03-15 | 核工业二〇八大队 | A kind of paleo-channel type uranium ore detection method |
| WO2017139271A2 (en) * | 2016-02-08 | 2017-08-17 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| CN107576982A (en) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect Comprehensive Seismic Prediction method |
| EP3376265A1 (en) * | 2017-03-14 | 2018-09-19 | CGG Services SAS | System and method for estimating the spatial distribution of an earth resource |
| CN108897041A (en) * | 2018-08-16 | 2018-11-27 | 中国石油天然气股份有限公司 | Prediction method and device for uranium ore enrichment area |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4081675A (en) * | 1976-11-08 | 1978-03-28 | Phillips Petroleum Company | Geophysical and geochemical exploration |
| US9316757B2 (en) * | 2011-12-06 | 2016-04-19 | Exxonmobil Upstream Research Company | Removal of fracture-induced anisotropy from converted-wave seismic amplitudes |
| US10101493B2 (en) * | 2016-08-26 | 2018-10-16 | Schlumberger Technology Corporation | Method for correcting natural gamma ray logging measurements |
| CN107367770A (en) * | 2017-06-12 | 2017-11-21 | 核工业北京地质研究院 | The technique for delineating of the favourable section of sandrock-type uranium-bearing Basin Metallogenic is soaked a kind ofly |
-
2018
- 2018-12-03 CN CN201811465517.2A patent/CN111257926B/en active Active
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104237964A (en) * | 2014-08-14 | 2014-12-24 | 核工业北京地质研究院 | Paleochannel type sandstone-type uranium ore deposit positioning method |
| WO2017139271A2 (en) * | 2016-02-08 | 2017-08-17 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| CN106501871A (en) * | 2016-09-14 | 2017-03-15 | 核工业二〇八大队 | A kind of paleo-channel type uranium ore detection method |
| EP3376265A1 (en) * | 2017-03-14 | 2018-09-19 | CGG Services SAS | System and method for estimating the spatial distribution of an earth resource |
| CN107576982A (en) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect Comprehensive Seismic Prediction method |
| CN108897041A (en) * | 2018-08-16 | 2018-11-27 | 中国石油天然气股份有限公司 | Prediction method and device for uranium ore enrichment area |
Non-Patent Citations (1)
| Title |
|---|
| "二连盆地中东部含铀古河谷构...建造及典型矿床成矿模式研究";刘波 等;《矿床地质》;20170228;第36卷(第1期);第126-142页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN111257926A (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN111257926B (en) | Method for predicting ancient valley uranium reservoir by using old seismic data | |
| CN101158724B (en) | Reservoir thickness prediction method based on dipolar wavelet | |
| CN106842307B (en) | Reservoir stratum fine prediction method based on waveform classification and retrieval under forward constraint | |
| CN112764122B (en) | Method for confining uranium reservoir of ancient river channel sandstone | |
| CN108802812A (en) | Well-seismic fusion stratum lithology inversion method | |
| CN109541685B (en) | River channel sand body identification method | |
| CN104047597A (en) | Fat gas mud shale stratum well log standardizing method | |
| CN103529475A (en) | Method for identifying and interpreting carbonate rock ancient karst reservoir layer three-dimensional structure | |
| CN105445800A (en) | Thick sand body top differentiation lithologic reservoir identification method | |
| KR101853334B1 (en) | Method for predicting exploitation site of unconventional resources applied to sequence stratigraphic analysis | |
| CN103345001A (en) | Ancient lake depth measuring method | |
| CN105842733A (en) | Shale reservoir earthquake identification method | |
| CN114910964B (en) | Prediction method for gravel rock mass dessert area on steep slope of fractured lake basin | |
| Iltaf et al. | Facies and petrophysical modeling of Triassic Chang 6 tight sandstone reservoir, Heshui oil field, Ordos basin, China | |
| CN112505754B (en) | Method for collaborative partitioning sedimentary microfacies by well-seismic based on high-precision sequence grid model | |
| CN111101935B (en) | Oil shale prediction method under few-well condition | |
| Wibowo et al. | Thinned coal distribution modeling based on integrated geological and geophysical data: Case study CBM resources in Central Palembang Sub-Basin | |
| CN115327627A (en) | Multi-information fusion tight sandstone gas diversion river channel characterization method and device | |
| CN109441438A (en) | A kind of structural complex lithologic boundary characterizing method | |
| CN113945993B (en) | Beach reservoir prediction method and device | |
| CN110795513B (en) | Method for predicting distribution of river facies source storage ectopic type compact oil gas dessert area | |
| CN112147676A (en) | Method for predicting thickness of coal bed and gangue | |
| Huang et al. | Subtle trap recognition based on seismic sedimentology—A case study from Shengli Oilfield | |
| Nwankwo et al. | Seismic Refraction Investigation for Groundwater Potential in Parts of Rivers State, Nigeria. | |
| CN117665933B (en) | Deposition microphase characterization method for well pattern area under low resolution seismic data limit |
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 |