CN108897041A - A kind of prediction technique and device of uranium ore enrichment region - Google Patents
A kind of prediction technique and device of uranium ore enrichment region Download PDFInfo
- Publication number
- CN108897041A CN108897041A CN201810932273.8A CN201810932273A CN108897041A CN 108897041 A CN108897041 A CN 108897041A CN 201810932273 A CN201810932273 A CN 201810932273A CN 108897041 A CN108897041 A CN 108897041A
- Authority
- CN
- China
- Prior art keywords
- research area
- porosity
- data
- shallow
- prediction
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
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. analysis, for interpretation, for correction
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/62—Physical property of subsurface
- G01V2210/624—Reservoir parameters
Abstract
The present invention provides the prediction techniques and device of a kind of uranium ore enrichment region.The prediction technique includes the following steps：Shallow seismic data processing is carried out to research area；Structure interpretation is carried out to the three-dimensional data for completing the processing of shallow-layer data, by faults comparison, obtains construction and the fracture in research area；According to well completion data, the statistics porosity in research area is obtained；According to the construction in research area and fracture, inverting is carried out to the porosity in research area, and compare with the statistics porosity in research area, obtain the uranium ore enrichment region in research area.The present invention also provides a kind of prediction meanss of uranium ore enrichment region, which includes：Seism processing module；Structure interpretation module；Obtain porosity module；Prediction module.Prediction technique and device of the invention can be effectively predicted sandstone-type uranium mineralization with respect enrichment region.
Description
Technical field
The present invention relates to a kind of method and apparatus using seismic technology prediction uranium ore enrichment region, belong to physical prospecting engineering neck
Domain.
Background technique
The enrichment process of sandstone-type uranium mineralization with respect is due to carrying U in fact6+Oxygen-containing underground water migrate forward along reservoir, migrating
During due to the effect U by reducing agent6+It is reduced to U4+To precipitating enrichment of ore-forming, this chemical reaction process by
To the strict control of reservoir heterogeneity, when the porosity of reservoir is very big, oxygen-containing underground water due to flow velocity is too fast cannot be abundant
Reduction occurs with reducing agent, and when porosity is too small, then the concordant flowing for blocking oxygen-containing underground water is more unfavorable for uranium
The formation of ore body.Therefore, the enrichment of prediction sandstone-type uranium mineralization with respect becomes a kind of needs.
It is traditional can ground leaching molding sand lithotype uranium exploration method mainly disposed in strict accordance with well pattern, explore successfully
Rate is relatively low, and by taking certain domestic sandrock-type uranium deposit as an example, the exploration success ratio of encrypted area is only 4 one-tenth or so, wastes a large amount of
Workload and fund.
Therefore, it is necessary to provide it is a kind of be suitble into the reservoir of mine porosity and predict to having, to predict uranium ore richness
The method of Ji Qu.
Summary of the invention
In order to solve the above-mentioned technical problem, the purpose of the present invention is to provide one kind can to sandstone-type uranium mineralization with respect enrichment region into
The method that row is effectively predicted.
In order to achieve the above technical purposes, the present invention provides a kind of prediction technique of uranium ore enrichment region, the prediction techniques
Include the following steps：
Step 1：Shallow seismic data processing is carried out to research area；
Step 2：Structure interpretation is carried out to the three-dimensional data for completing shallow seismic data processing, by faults comparison, is obtained
Study construction and the fracture in area；
Step 3：According to well completion data, the statistics porosity ranges in research area are obtained；
Step 4：According to the construction in research area and fracture, inverting, and the system with research area are carried out to the porosity in research area
Meter porosity ranges compare, and obtain the uranium ore enrichment region in research area.
The prediction technique of uranium ore enrichment region of the invention is by the statistics porosity ranges of formation of sandstone-type uranium deposits band and hole
Porosity inversion technique combines, and realizes and can use three dimensional seismic data in the case where no drilling well to Beneficial Ore-forming band
Enrichment region is accurately predicted.
The present invention also provides a kind of prediction meanss of uranium ore enrichment region, which includes：
Seism processing module, for carrying out shallow seismic data processing to research area；
Structure interpretation module passes through layer for carrying out structure interpretation to the three-dimensional data for completing shallow seismic data processing
Position comparison obtains construction and the fracture in research area；
Porosity module is obtained, for obtaining the statistics porosity ranges in research area according to well completion data；
Prediction module, for according to research area construction and fracture, to research area porosity carry out inverting, and with research
The statistics porosity ranges in area compare, and obtain the uranium ore enrichment region in research area.
The prediction technique and device of uranium ore enrichment region of the invention, solve sandstone-type uranium mineralization with respect by seismic processing technique
The problems such as shallow earthquake data SNR is high, surface wave sound wave and linear disturbance are developed；It is implemented by seismic interpretation technique shallow
The structure development situation of layer；Counted and studied by the Heterogeneous Characteristics to metallogenic belt uranium reservoir, determine be conducive to uranium at
The porosity ranges of mine are effectively predicted using enrichment region of the porosity inversion to metallogenic belt.
The prediction technique and device of uranium ore enrichment region of the invention combine can ground leaching type sandrock-type uranium deposit the characteristics of, by stone
Seismic technology in exploration activity has sufficiently been applied in the exploration of sandstone-type uranium mineralization with respect, improves exploration efficiency, reduce exploration at
This.In terms of shallow seismic data processing, since the reservoir lithology of sandstone-type uranium mineralization with respect is more single, predominantly sand-mud interbed,
For this feature, by conventional petroleum explore in seismic processing technique targetedly simplified and improved, shorten place
The reason period reduces runing time；In terms of porosity inversion, the heterogeneity for carrying out reservoir to all prospect pits in metallogenic belt is united
Meter sums up the aeolotropic characteristics and rule in research area, the accuracy rate of inversion result can be improved.
Detailed description of the invention
Fig. 1 is the structural schematic diagram of an one of the specific embodiment of the invention prediction meanss of uranium ore enrichment region.
Fig. 2 is the flow chart of one of another specific embodiment of present invention prediction technique of uranium ore enrichment region.
Fig. 3 is that GR-POR in another specific embodiment of the present invention crosses figure.
Fig. 4 is the porosity inversion plan view in another specific embodiment of the present invention.
Specific embodiment
In order to which technical characteristic of the invention, purpose and beneficial effect are more clearly understood, now to skill of the invention
Art scheme carries out described further below, but should not be understood as that limiting the scope of the invention.
As shown in Figure 1, providing a kind of prediction meanss of uranium ore enrichment region in a specific embodiment of the invention, this is pre-
Surveying device includes：
Seism processing module, for carrying out shallow seismic data processing to research area；
Structure interpretation module passes through layer for carrying out structure interpretation to the three-dimensional data for completing shallow seismic data processing
Position comparison obtains construction and the fracture in research area；
Porosity module is obtained, for obtaining the statistics porosity ranges in research area according to well completion data；
Prediction module, for according to research area construction and fracture, to research area porosity carry out inverting, and with research
The statistics porosity ranges in area compare, and obtain the uranium ore enrichment region in research area.
Sandstone-type uranium mineralization with respect it is shallower at mine buried depth, it is more smart to carry out that specific aim processing is carried out to the seismic data of shallow-layer
Thin reliable seismic interpretation.Specifically, in seism processing module, when carrying out shallow seismic data processing to research area,
Signal-to-noise ratio is mainly improved by the following means：
Combining excision scanning superposition rationally to retain shallow-layer along direction in space application main road collection seismic data subhead graticule has
Information is imitated, signal-to-noise ratio technology is improved by binning homogenization, prestack and improves shallow-layer data signal-to-noise ratio, frequency technology is opened up using poststack and improves
The resolution ratio of data.
More specifically, subhead graticule combines excision scanning superimposing technique along direction in space application main road collection, rationally retain shallow
Layer effective information.The difficult point that shallow earthquake data moves school excision is to have cut greatly completely, but shallow-layer just disappears；Dynamic school is cut small
Stretching directly affects Overlay.Therefore, effective shallow-layer information in order to obtain scans uniform velocity analysis method in application speed
Under conditions of Accurate Analysis speed, using subhead graticule along direction in space big trace gather define cut off on the basis of, carry out excision sweep
Retouch superposition.It must accomplish in processing：First to greatest extent with small excision；Second cuts off nmo stretching completely.This requires
Accomplish finely to cut off and fine velocity analysis successive ignition in processing.
More specifically, improving the signal-to-noise ratio of shallow-layer data by binning homogenization technology.Shallow-layer data degree of covering is low, believes
It makes an uproar than low.By the degree of covering of shallow-layer data being properly increased, to properly increase shallow-layer in the minimum homogenizing distance parameter of selection
The signal-to-noise ratio of data.Binning homogenization is specifically done in the range of 0-800 meters of geophone offsets, homogenizing distance is 24 meters.It is that face element is equal
Change the superposition of front and back, it can be seen that the signal-to-noise ratio of shallow-layer is increased after homogenizing.
More specifically, improving the signal-to-noise ratio that signal-to-noise ratio technology improves shallow-layer data by prestack.It is protected in both the above method
On the basis of staying shallow-layer information and improving shallow-layer signal-to-noise ratio, the linear disturbance, surface wave interference, arteries and veins being commonly present also are rejected in prestack
Punching interference and powerful random noise.Strong amplitude interference for these abnormal strong amplitudes and part, if do not rejected, except to folded
Outside adduction migration result adversely affects, also there is great side effect to the fidelity processing of amplitude.
For example, the signal-to-noise ratio of shallow-layer data can be improved using high-precision Radon converter technique multiple suppression.This method
Main thought be the primary wave smoothing with after the dynamic correction of wave velocity, multiple wave undercorrection and be a parabola passes through
Radon transformation is completely separable with multiple wave in domain primary wave and is removed, then carries out inverse transformation, reaches the mesh of Multiple attenuation
, this method can preferably multiple suppression.It is especially preferable to shortcut pressing result, and the effect with the anti-alias of high-fidelity
Fruit.
More specifically, opening up the resolution ratio that frequency technology improves data using poststack.By applying prestack inverse Q filtering, earth's surface one
Cause property deconvolution, THE MULTICHANNEL PREDICTIVE DECONVOLUTION, the frequency band of data are effectively widened；Using prestack multiple domain high-fidelity noise elimination technology,
Prestack high-fidelity multiple elimination technology has been suppressed various strong jamming waves, the signal-to-noise ratio of data present in data very well and has been obtained very
It is big to improve.
Specifically, in structure interpretation module, structure interpretation is carried out to the three-dimensional data for completing the processing of shallow-layer data, is passed through
Faults comparison obtains construction and the fracture in research area.
Seismic data signal-to-noise ratio with higher by seism processing resume module, data frequency band are effectively opened up
Width carries out the available fine reliable explanation results of structure interpretation on this basis.
More specifically, it is more than exploration drilling number of the sandstone-type uranium mineralization with respect containing mining area and intensive, drilling data can be made full use of
Carry out composite traces horizon calibration.It the characteristics of according to research area's seismic data, according to seismic reflection configuration, wave group feature, utilizes
The means such as section, slice, block drift, carry out the explanation of layer position；Identification, relevant and water are combined with temporal profile using normal profile
The means such as the extraction of flat slice carry out fault interpretation.To obtain fine seismic data interpretation result.
Specifically, in obtaining porosity module, according to well completion data, the statistics porosity ranges in research area are obtained.Only
Have in the aperture ranges of appropriateness, enough oxygen-containing uranium containing waters could sufficiently be acted on reducing substances, and uranium is made to precipitate, and be formed
Uranium deposit, therefore excessively high porosity and too low porosity are all unfavorable at mine.
More specifically, being to be counted to the porosity of research area's target zone, and sum up reservoir porosity and uranium mineralization
Linear relationship between grade, and then the Appreciation gist as porosity inversion result.For example, in uranium ore block containing mining area, complete well
Abundant information, applied analysis chemically examines means and carries out lacunarity analysis to rock core information, by statistical result it can be seen that uranium ore
The porosity of abundance zone reservoir is mostly in 20%-30%.To the mineral deposit with certain particularity, can be united according to reapective features
Meter.
But rock core information is relatively limited in actual exploration, makes full use of limited rock core information test porosity,
It is compared with log, establishes the correlativity of log and core data, be generalized to the prospect pit for lacking rock core information.
Adopting said method can use well-log information and calculate the reservoir porosity of all prospect pits in research area, thus more acurrate
The porosity ranges for counting uranium ore abundance zone reservoir.
Specifically, in prediction module, according to the construction in research area and fracture, inverting is carried out to the porosity in research area,
And compared with the statistics porosity ranges in research area, obtain the uranium ore enrichment region in research area.
Porosity inversion can accurately depict the porosity distribution characteristics of metallogenic belt reservoir.Uranium deposit enrichment region is general
All there is more prospect pit, drilling data is abundant, and the porosity inversion of reservoir can be carried out using Application of Logging-constrained Inversion.Well logging is about
Beam inverting belongs to earthquake post-stack inversion, obtains Acoustic Impedance Data using convolution model based on being a kind of poststack data by earthquake
Inversion method.
More specifically, the inversion method generates initial reflection coefficient series based on convolution model, according to well-log information
R1(t), post-stack seismic data obtains reflection seimogram x (t), extracted using the wavelet extraction of individual well multiple tracks or statistical wavelet etc.
Method obtains primary earthquake wavelet b1(t), by y1(t)=b1(t)R1(t) artificial synthesized record y is obtained1(t), to x (t) and y1
(t) cross-correlation is done,Wherein RxyFor related coefficient.By modifying y1(t) make RxyAs far as possible
Greatly, that is, primary earthquake wavelet b is modified1(t) and initial reflection coefficient series R1(t).Work as RxyWhen reaching satisfactory value, to density, speed
Degree and reflection coefficient establish initial model using interpolating methods such as anti-square distance, triangular mesh and Ke Lijin, generate initial
Wave Impedance Data Volume.Above procedure is a forward modeling process, it passes through modification log, initial reflection coefficient and initial son
Wave generates reasonable impedance initial value model.On the basis of initial model, to the wave impedance of all interpolations according to conjugate gradient
Method carries out limited times modification in certain variation range, reaches objective function e1(t)=x (t)-y1(t) minimal point generates most
Whole inverting section.Wherein e1(t) degree of agreement of representative model and earthquake record.Above procedure is a refutation process, it
The wave impedance for modifying initial model, its purpose is to final mask and earthquake record are coincide.
Application of Logging-constrained Inversion is organically to combine earthquake with well logging, breaches seismic resolution in traditional sense
Limitation theoretically can be obtained resolution ratio identical with well-log information, be the key technology of fine description lithology.Inverting is required from ground
3 shake wavelet, well-log information and initial model aspects do elaboration.Wavelet is the bridge of building well logging and Earthquakes,
Good wavelet should waveform stabilization, energy is concentrated mainly on the main lobe of wavelet, and side-lobe energy is small and decays rapidly.Well logging money
Material, especially sound wave and density log, are the basic foundations for establishing the basic data and geologic interpretation of initial model, it should be noted that disappear
Unless the influence of geologic(al) factor.Initial model is the fundamental way for reducing final result multi-solution, needs to believe with known geology
Breath constantly comparison, establishes the surge impedance model as close possible to actual formation situation.
Finally, the statistics porosity ranges that area will be studied, the porosity of the 3-D seismics area of coverage obtained with seismic inversion
Predicted value is compared, and the overlapping region of the two is uranium ore enrichment region.
System, device, module or the unit that above-described embodiment illustrates can specifically realize by computer chip or entity,
Or it is realized by the product with certain function.
For convenience of description, it is divided into various modules when description apparatus above with function to describe respectively.Certainly, implementing this
The function of each module can be realized in the same or multiple software and or hardware when application.
As seen through the above description of the embodiments, those skilled in the art can be understood that the application can
It realizes by means of software and necessary general hardware platform.Based on this understanding, the technical solution essence of the application
On in other words the part that contributes to existing technology can be embodied in the form of software products, in a typical configuration
In, calculating equipment includes one or more processors (CPU), input/output interface, network interface and memory.
All the embodiments in this specification are described in a progressive manner, same and similar portion between each embodiment
Dividing may refer to each other, and each embodiment focuses on the differences from other embodiments.Especially for system reality
For applying example, since it is substantially similar to the method embodiment, so being described relatively simple, related place is referring to embodiment of the method
Part explanation.
The present invention can describe in the general context of computer-executable instructions executed by a computer, such as program
Module.Generally, program module includes routines performing specific tasks or implementing specific abstract data types, programs, objects, group
Part, data structure etc..The application can also be practiced in a distributed computing environment, in these distributed computing environments, by
Task is executed by the connected remote processing devices of communication network.In a distributed computing environment, program module can be with
In the local and remote computer storage media including storage equipment.
As shown in Fig. 2, a kind of prediction technique of uranium ore enrichment region is provided in another embodiment of the present invention, it should
Prediction technique includes the following steps：
Step S1：Shallow seismic data processing is carried out to research area；
Step S2：Structure interpretation is carried out to the three-dimensional data for completing the processing of shallow-layer data to be studied by faults comparison
The construction in area and fracture；
Step S3：According to well completion data, the statistics porosity ranges in research area are obtained；
Step S4：According to the construction in research area and fracture, inverting, and the system with research area are carried out to the porosity in research area
Meter porosity ranges compare, and obtain the uranium ore enrichment region in research area.
Sandstone-type uranium mineralization with respect it is shallower at mine buried depth, it is more smart to carry out that specific aim processing is carried out to the seismic data of shallow-layer
Thin reliable seismic interpretation.Specifically, in step sl, when carrying out shallow seismic data processing to research area, mainly pass through
Following means improve signal-to-noise ratio：
Combining excision scanning superposition rationally to retain shallow-layer along direction in space application main road collection seismic data subhead graticule has
Information is imitated, signal-to-noise ratio technology is improved by binning homogenization, prestack and improves shallow-layer data signal-to-noise ratio, frequency technology is opened up using poststack and improves
The resolution ratio of data.
More specifically, subhead graticule combines excision scanning superimposing technique along direction in space application main road collection, rationally retain shallow
Layer effective information.The difficult point that shallow earthquake data moves school excision is to have cut greatly completely, but shallow-layer just disappears；Dynamic school is cut small
Stretching directly affects Overlay.Therefore, effective shallow-layer information in order to obtain scans uniform velocity analysis method in application speed
Under conditions of Accurate Analysis speed, using subhead graticule along direction in space big trace gather define cut off on the basis of, carry out excision sweep
Retouch superposition.It must accomplish in processing：First to greatest extent with small excision；Second cuts off nmo stretching completely.This requires
Accomplish finely to cut off and fine velocity analysis successive ignition in processing.
More specifically, improving the signal-to-noise ratio of shallow-layer data by binning homogenization technology.Shallow-layer data degree of covering is low, believes
It makes an uproar than low.By the degree of covering of shallow-layer data being properly increased, to properly increase shallow-layer in the minimum homogenizing distance parameter of selection
The signal-to-noise ratio of data.Binning homogenization is specifically done in the range of 0-800 meters of geophone offsets, homogenizing distance is 24 meters.It is that face element is equal
Change the superposition of front and back, it can be seen that the signal-to-noise ratio of shallow-layer is increased after homogenizing.
More specifically, improving the signal-to-noise ratio that signal-to-noise ratio technology improves shallow-layer data by prestack.It is protected in both the above method
On the basis of staying shallow-layer information and improving shallow-layer signal-to-noise ratio, the linear disturbance, surface wave interference, arteries and veins being commonly present also are rejected in prestack
Punching interference and powerful random noise.Strong amplitude interference for these abnormal strong amplitudes and part, if do not rejected, except to folded
Outside adduction migration result adversely affects, also there is great side effect to the fidelity processing of amplitude.
For example, the signal-to-noise ratio of shallow-layer data can be improved using high-precision Radon converter technique multiple suppression.This method
Main thought be the primary wave smoothing with after the dynamic correction of wave velocity, multiple wave undercorrection and be a parabola passes through
Radon transformation is completely separable with multiple wave in domain primary wave and is removed, then carries out inverse transformation, reaches the mesh of Multiple attenuation
, this method can preferably multiple suppression.It is especially preferable to shortcut pressing result, and the effect with the anti-alias of high-fidelity
Fruit.
More specifically, opening up the resolution ratio that frequency technology improves data using poststack.By applying prestack inverse Q filtering, earth's surface one
Cause property deconvolution, THE MULTICHANNEL PREDICTIVE DECONVOLUTION, the frequency band of data are effectively widened；Using prestack multiple domain high-fidelity noise elimination technology,
Prestack high-fidelity multiple elimination technology has been suppressed various strong jamming waves, the signal-to-noise ratio of data present in data very well and has been obtained very
It is big to improve.
Specifically, in step s 2, structure interpretation is carried out to the three-dimensional data for completing the processing of shallow-layer data, it is right by layer position
Than obtaining construction and the fracture in research area.
By the seismic data signal-to-noise ratio with higher that step S1 is handled, data frequency band is effectively widened, in this base
The available fine reliable explanation results of structure interpretation are carried out on plinth.
More specifically, it is more than exploration drilling number of the sandstone-type uranium mineralization with respect containing mining area and intensive, drilling data can be made full use of
Carry out composite traces horizon calibration.It the characteristics of according to research area's seismic data, according to seismic reflection configuration, wave group feature, utilizes
The means such as section, slice, block drift, carry out the explanation of layer position；Identification, relevant and water are combined with temporal profile using normal profile
The means such as the extraction of flat slice carry out fault interpretation.To obtain fine seismic data interpretation result.
Specifically, in step s3, according to well completion data, the statistics porosity ranges in research area are obtained.Only in appropriateness
Aperture ranges in, enough oxygen-containing uranium containing waters could sufficiently be acted on reducing substances, so that uranium is precipitated, formed uranium deposit,
Therefore excessively high porosity and too low porosity are all unfavorable at mine.
More specifically, being to be counted to the porosity of research area's target zone, and sum up reservoir porosity and uranium mineralization
Linear relationship between grade, and then the Appreciation gist as porosity inversion result.For example, in uranium ore block containing mining area, complete well
Abundant information, applied analysis chemically examines means and carries out lacunarity analysis to rock core information, by statistical result it can be seen that uranium ore
The porosity of abundance zone reservoir is mostly in 20%-30%.To the mineral deposit with certain particularity, can be united according to reapective features
Meter.
But rock core information is relatively limited in actual exploration, makes full use of limited rock core information test porosity,
It is compared with log, establishes the correlativity of log and core data, be generalized to the prospect pit for lacking rock core information.
Adopting said method can use well-log information and calculate the reservoir porosity of all prospect pits in research area, thus more acurrate
The porosity ranges for counting uranium ore abundance zone reservoir.
Specifically, in step s 4, according to the construction in research area and fracture, inverting is carried out to the porosity in research area, and
It is compared with the statistics porosity ranges in research area, obtains the uranium ore enrichment region in research area.
Porosity inversion can accurately depict the porosity distribution characteristics of metallogenic belt reservoir.Uranium deposit enrichment region is general
All there is more prospect pit, drilling data is abundant, and the porosity inversion of reservoir can be carried out using Application of Logging-constrained Inversion.Well logging is about
Beam inverting belongs to earthquake post-stack inversion, obtains Acoustic Impedance Data using convolution model based on being a kind of poststack data by earthquake
Inversion method.
More specifically, the inversion method generates initial reflection coefficient series based on convolution model, according to well-log information
R1(t), post-stack seismic data obtains reflection seimogram x (t), extracted using the wavelet extraction of individual well multiple tracks or statistical wavelet etc.
Method obtains primary earthquake wavelet b1(t), by y1(t)=b1(t)R1(t) artificial synthesized record y is obtained1(t), to x (t) and y1
(t) cross-correlation is done,Wherein RxyFor related coefficient.By modifying y1(t) make RxyAs far as possible
Greatly, that is, primary earthquake wavelet b is modified1(t) and initial reflection coefficient series R1(t).Work as RxyWhen reaching satisfactory value, to density, speed
Degree and reflection coefficient establish initial model using interpolating methods such as anti-square distance, triangular mesh and Ke Lijin, generate initial
Wave Impedance Data Volume.Above procedure is a forward modeling process, it passes through modification log, initial reflection coefficient and initial son
Wave generates reasonable impedance initial value model.On the basis of initial model, to the wave impedance of all interpolations according to conjugate gradient
Method carries out limited times modification in certain variation range, reaches objective function e1(t)=x (t)-y1(t) minimal point generates most
Whole inverting section.Wherein e1(t) degree of agreement of representative model and earthquake record.Above procedure is a refutation process, it
The wave impedance for modifying initial model, its purpose is to final mask and earthquake record are coincide.
Application of Logging-constrained Inversion is organically to combine earthquake with well logging, breaches seismic resolution in traditional sense
Limitation theoretically can be obtained resolution ratio identical with well-log information, be the key technology of fine description lithology.Inverting is required from ground
3 shake wavelet, well-log information and initial model aspects do elaboration.Wavelet is the bridge of building well logging and Earthquakes,
Good wavelet should waveform stabilization, energy is concentrated mainly on the main lobe of wavelet, and side-lobe energy is small and decays rapidly.Well logging money
Material, especially sound wave and density log, are the basic foundations for establishing the basic data and geologic interpretation of initial model, it should be noted that disappear
Unless the influence of geologic(al) factor.Initial model is the fundamental way for reducing final result multi-solution, needs to believe with known geology
Breath constantly comparison, establishes the surge impedance model as close possible to actual formation situation.
Finally, the statistics porosity ranges that area will be studied, the porosity of the 3-D seismics area of coverage obtained with seismic inversion
Predicted value is compared, and the overlapping region of the two is uranium ore enrichment region.
Embodiment 1
A kind of method for predicting the shop Qian Jia uranium deposit is present embodiments provided, following steps are specifically included：
Predict that the uranium deposit has drilled hundreds of mouthfuls of prospect pits in prior prospect, with the continuous improvement of degree of prospecting, flared end at
Power is gradually reduced, and exploration cost rises obvious.
Step 1：Signal-to-noise ratio technology is improved using binning homogenization, prestack and opens up frequency technical method using poststack, to the uranium
Mineral deposit 400km2Seismic data carry out shallow seismic data processing, improve data resolution reduce signal-to-noise ratio.So that exploration mesh
Layer Upper Cretaceous Yao Jia group reflectance signature it is clear and legible.
Step 2：Hundreds of mouthfuls of prospect pits of the uranium deposit finishing drilling are carried out with the production of composite traces, and layer position is demarcated.
The method that reflectance signature and wave resistance characteristic use section and slice according to seismic data combine is to exploration target zone Yao Jia group
Fine explanation is carried out, construction feature is implemented.
Step 3：Analysis is sampled to the rock core of the uranium deposit, the rock core for obtaining available core porosity routinely divides
Analysis report, each rock core conventional analysis fruit is compared with the porosity logging curve (POR) of the well, establishes the uranium deposit hole
The correlativity of porosity log and core porosity data, by the hole of the correlativity of foundation and the drilling well of every mouthful of the mineral deposit
Degree log combines, and counts the reservoir porosity range of all drilled well Yao Jia groups in the mineral deposit, statistical result shows to reach
Core porosity range to industrial standard is 20%-30% (as shown in Figure 3).
Step 4：Since the mineral deposit drilled well number is more, the porosity of Yao Jia group is carried out using the method for Application of Logging-constrained Inversion
Prediction can show the porosity prediction of Yao Jia group by slice (Fig. 4) as a result, the Regional Representative of 20-30 is in being sliced
Porosity is in the region of 20%-30%, and according to the statistical result in step 3, the region 20-30 is exactly the uranium ore enrichment predicted
Area.
Exploration deployment has been carried out according to the prediction result, and exploration success ratio is close to 80%.
Above embodiments explanation directs the shop Qian Jia mineral deposit using the method for seismic technology prediction uranium ore enrichment region
Flared end exploration, the northern deployment prospect pit mouths up to a hundred in encrypted area, success rate close to 80%, relatively before 40% success rate improve
40 percentage points, exploration cost has been saved, has saved exploration fund.
Claims (10)
1. a kind of prediction technique of uranium ore enrichment region, which is characterized in that the prediction technique includes the following steps：
Step 1：Shallow seismic data processing is carried out to research area；
Step 2：Structure interpretation is carried out to the three-dimensional data for completing shallow seismic data processing to be studied by faults comparison
The construction in area and fracture；
Step 3：According to well completion data, the statistics porosity ranges in research area are obtained；
Step 4：According to it is described research area construction and fracture, to research area porosity carry out inverting, and with the research area
Statistics porosity ranges compare, obtain research area uranium ore enrichment region.
2. prediction technique according to claim 1, which is characterized in that in said step 1, the shallow earthquake data
Processing includes：Combining excision scanning superposition rationally to retain shallow-layer along direction in space application main road collection seismic data subhead graticule has
Information is imitated, signal-to-noise ratio technology is improved by binning homogenization, prestack and improves shallow-layer data signal-to-noise ratio, frequency technology is opened up using poststack and improves
The resolution ratio of data.
3. prediction technique according to claim 2, which is characterized in that it is high-precision that the prestack, which improves signal-to-noise ratio technology,
Radon converter technique multiple suppression.
4. prediction technique according to claim 1, which is characterized in that in the step 2, press when structure interpretation
It is carried out according to following steps：
The explanation of layer position is carried out using section, slice, block drift means according to seismic reflection configuration, wave group feature；Utilize routine
Section combines identification with temporal profile, relevant and dropping cut slice extraction means carry out fault interpretation, obtains fine earthquake
Data interpretation result.
5. prediction technique according to claim 1, which is characterized in that in the step 3, obtain the statistics in research area
When porosity ranges, follow the steps below：
The porosity of metallogenic belt section containing mine rock core is obtained using lab analysis means；
Using well completion data, the curve of log and core porosity is established；
The statistics porosity ranges in research area are obtained according to the curve.
6. prediction technique according to claim 1, which is characterized in that in the step 4, to the porosity in research area
Carried out when inverting using Application of Logging-constrained Inversion.
7. a kind of prediction meanss of uranium ore enrichment region, which is characterized in that the prediction meanss include：
Seism processing module, for carrying out shallow seismic data processing to research area；
Structure interpretation module is right by layer position for carrying out structure interpretation to the three-dimensional data for completing shallow seismic data processing
Than obtaining construction and the fracture in research area；
Porosity module is obtained, for obtaining the statistics porosity ranges in research area according to well completion data；
Prediction module, for according to it is described research area construction and fracture, to research area porosity carry out inverting, and with it is described
The statistics porosity ranges in research area compare, and obtain the uranium ore enrichment region in research area.
8. prediction meanss according to claim 7, which is characterized in that described shallow in the seism processing module
Layer seism processing include：Combine excision scanning superposition reasonable along direction in space application main road collection seismic data subhead graticule
Retain shallow-layer effective information, signal-to-noise ratio technology is improved by binning homogenization, prestack and improves shallow-layer data signal-to-noise ratio, is opened up using poststack
The resolution ratio of frequency technology raising data.
9. prediction meanss according to claim 8, which is characterized in that it is high-precision that the prestack, which improves signal-to-noise ratio technology,
Radon converter technique multiple suppression.
10. prediction meanss according to claim 7, which is characterized in that in the structure interpretation module, carry out construction solution
It is followed the steps below when releasing：
The explanation of layer position is carried out using section, slice, block drift means according to seismic reflection configuration, wave group feature；Utilize routine
Section combines identification with temporal profile, relevant and dropping cut slice extraction means carry out fault interpretation, obtains fine earthquake
Data interpretation result.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810932273.8A CN108897041B (en) | 2018-08-16 | 2018-08-16 | Prediction method and device for uranium ore enrichment area |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201810932273.8A CN108897041B (en) | 2018-08-16 | 2018-08-16 | Prediction method and device for uranium ore enrichment area |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN108897041A true CN108897041A (en) | 2018-11-27 |
| CN108897041B CN108897041B (en) | 2020-06-09 |
Family
ID=64354412
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201810932273.8A Active CN108897041B (en) | 2018-08-16 | 2018-08-16 | Prediction method and device for uranium ore enrichment area |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN108897041B (en) |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN109580453A (en) * | 2018-12-24 | 2019-04-05 | 核工业北京地质研究院 | The method for calculating sandstone-type uranium mineralization with respect sand body porosity based on optimum logging interpretation |
| CN111239815A (en) * | 2020-01-20 | 2020-06-05 | 吉林大学 | Sandstone-type uranium reservoir mineralization and deposition element extraction method based on three-dimensional seismic attributes |
| CN111257926A (en) * | 2018-12-03 | 2020-06-09 | 核工业二0八大队 | Method for predicting ancient valley uranium reservoir by using old seismic data |
| CN111852467A (en) * | 2020-07-28 | 2020-10-30 | 核工业北京地质研究院 | Method and system for delineating extension range of sandstone uranium ore body |
| CN112764122A (en) * | 2019-10-21 | 2021-05-07 | 中国石油化工股份有限公司 | Delineation method of ancient river channel type sandstone uranium reservoir |
| CN113514886A (en) * | 2021-07-22 | 2021-10-19 | 核工业北京地质研究院 | Geological-seismic three-dimensional prediction method for beneficial part of sandstone-type uranium deposit mineralization |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090028000A1 (en) * | 2007-07-26 | 2009-01-29 | O'brien Thomas B | Method and process for the systematic exploration of uranium in the athabasca basin |
| CN103954995A (en) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | Sand body reorganization method in sandstone-type uranium deposit exploration |
| CN107144873A (en) * | 2017-04-12 | 2017-09-08 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect 3D seismic data observation procedure |
| CN107576982A (en) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect Comprehensive Seismic Prediction method |
| CN107664772A (en) * | 2017-07-28 | 2018-02-06 | 核工业北京地质研究院 | A kind of Forecasting Methodology for sandstone-type uranium mineralization with respect favorab1e sandbadies and Favorable Zones |
-
2018
- 2018-08-16 CN CN201810932273.8A patent/CN108897041B/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20090028000A1 (en) * | 2007-07-26 | 2009-01-29 | O'brien Thomas B | Method and process for the systematic exploration of uranium in the athabasca basin |
| CN103954995A (en) * | 2014-04-22 | 2014-07-30 | 核工业北京地质研究院 | Sand body reorganization method in sandstone-type uranium deposit exploration |
| CN107144873A (en) * | 2017-04-12 | 2017-09-08 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect 3D seismic data observation procedure |
| CN107664772A (en) * | 2017-07-28 | 2018-02-06 | 核工业北京地质研究院 | A kind of Forecasting Methodology for sandstone-type uranium mineralization with respect favorab1e sandbadies and Favorable Zones |
| CN107576982A (en) * | 2017-09-01 | 2018-01-12 | 核工业北京地质研究院 | A kind of sandstone-type uranium mineralization with respect Comprehensive Seismic Prediction method |
Non-Patent Citations (7)
| Title |
|---|
| 卢佳岚，等: "测井约束地震反演技术在泌阳凹陷赵凹地区地震资料解释中的应用", 《河南石油》 * |
| 吴曲波，等: "砂岩型铀矿地震勘探技术应用现状与发展", 《物探与化探》 * |
| 吴曲波，等: "砂岩型铀矿地震综合预测技术研究-以齐哈日格图凹陷铀矿为例", 《地球物理学进展》 * |
| 孙国红，等: "喇萨杏油田萨零组油层有效孔隙度与原始含油饱和度精细解释", 《大庆石油地质与开发》 * |
| 李立尧，等: "煤系承压含水层静储量计算", 《煤炭技术》 * |
| 薛志恒: "浅层地震勘探在砂岩铀矿勘查中的应用", 《河南理工大学学报(自然科学版)》 * |
| 谈顺佳，等: "有利含铀砂岩体的地球物理响应特征研究-以松辽盆地钱家店铀矿床为例", 《地球物理学进展》 * |
Cited By (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN111257926A (en) * | 2018-12-03 | 2020-06-09 | 核工业二0八大队 | Method for predicting ancient valley uranium reservoir by using old seismic data |
| CN111257926B (en) * | 2018-12-03 | 2022-07-26 | 核工业二0八大队 | Method for predicting ancient valley uranium reservoir by using old seismic data |
| CN109580453A (en) * | 2018-12-24 | 2019-04-05 | 核工业北京地质研究院 | The method for calculating sandstone-type uranium mineralization with respect sand body porosity based on optimum logging interpretation |
| CN112764122A (en) * | 2019-10-21 | 2021-05-07 | 中国石油化工股份有限公司 | Delineation method of ancient river channel type sandstone uranium reservoir |
| CN112764122B (en) * | 2019-10-21 | 2023-12-12 | 中国石油化工股份有限公司 | Method for confining uranium reservoir of ancient river channel sandstone |
| CN111239815A (en) * | 2020-01-20 | 2020-06-05 | 吉林大学 | Sandstone-type uranium reservoir mineralization and deposition element extraction method based on three-dimensional seismic attributes |
| CN111239815B (en) * | 2020-01-20 | 2021-07-13 | 吉林大学 | Sandstone-type uranium reservoir mineralization and deposition element extraction method based on three-dimensional seismic attributes |
| CN111852467A (en) * | 2020-07-28 | 2020-10-30 | 核工业北京地质研究院 | Method and system for delineating extension range of sandstone uranium ore body |
| CN111852467B (en) * | 2020-07-28 | 2023-04-07 | 核工业北京地质研究院 | Method and system for delineating extension range of sandstone uranium ore body |
| CN113514886A (en) * | 2021-07-22 | 2021-10-19 | 核工业北京地质研究院 | Geological-seismic three-dimensional prediction method for beneficial part of sandstone-type uranium deposit mineralization |
| CN113514886B (en) * | 2021-07-22 | 2021-12-10 | 核工业北京地质研究院 | Geological-seismic three-dimensional prediction method for beneficial part of sandstone-type uranium deposit mineralization |
Also Published As
| Publication number | Publication date |
|---|---|
| CN108897041B (en) | 2020-06-09 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN108897041A (en) | A kind of prediction technique and device of uranium ore enrichment region | |
| Bauer et al. | Deep structure of the Namibia continental margin as derived from integrated geophysical studies | |
| CN102466815A (en) | Triassic clastic rock petroleum reservoir identification method | |
| Maries et al. | Cross-profile seismic data acquisition, imaging, and modeling of iron-oxide deposits: A case study from Blötberget, south-central Sweden | |
| Beckel et al. | The cross-dip correction as a tool to improve imaging of crooked-line seismic data: a case study from the post-glacial Burträsk fault, Sweden | |
| CN107269269A (en) | The quick look method of unconformable strata layer position | |
| Farrell et al. | Refraction statics | |
| Rafiq et al. | Reservoir characterization using microseismic facies analysis integrated with surface seismic attributes | |
| Ding et al. | Reverse time migration (RTM) imaging of iron-oxide deposits in the Ludvika mining area, Sweden | |
| Sokolov et al. | Empirical study of sediment-filled basin response: The case of Taipei City | |
| Chandler et al. | Utility of the horizontal-to-vertical spectral ratio passive seismic method for estimating thickness of Quaternary sediments in Minnesota and adjacent parts of Wisconsin | |
| Kim et al. | Imaging of reflection seismic energy for mapping shallow fracture zones in crystalline rocks | |
| Gao et al. | Emplacement of intrusions of the Tarim Flood Basalt Province and their impacts on oil and gas reservoirs: A 3D seismic reflection study in Yingmaili fields, Tarim Basin, northwest China | |
| Wang et al. | Surface waves in the western Taiwan coastal plain from an aftershock of the 1999 Chi-Chi, Taiwan, earthquake | |
| Chacón-Hernández et al. | Analysis of shear wave splitting parameters in los Humeros geothermal field, Puebla, Mexico | |
| Stephenson et al. | High-resolution seismic reflection surveys and modeling across an area of high damage from the 1994 Northridge earthquake, Sherman Oaks, California | |
| Poppeliers et al. | The seismic response of a steep slope: High-resolution observations with a dense, three-component seismic array | |
| Behera et al. | Imaging Mesozoic sediments in Deccan Volcanic province of India: inferences from seismic and gravity studies | |
| Martinez et al. | Utility of ground-penetrating radar in near-surface, high-resolution imaging of Lansing-Kansas City (Pennsylvanian) limestone reservoir analogs | |
| Mata García et al. | Seismic characterisation using artificial intelligence algorithms for rock property distribution in static modeling | |
| Ebel et al. | A geophysical and field survey for the source region of the 1638 New Hampshire earthquake | |
| Rashed et al. | High‐resolution shallow seismic and ground penetrating radar investigations revealing the evolution of the Uemachi Fault system, Osaka, Japan | |
| Zhao et al. | Interference between seismic body wave and secondary surface wave resulting in the peak collapse ratios of buildings | |
| RU2093860C1 (en) | Method of seismic test of change of capacitive properties and position of productivity contour of oil and gas pool process of its development | |
| Rucker | Constraining boundaries and extent of the Charleston uplift, Northeast New Madrid Seismic Zone, using shallow seismic reflection methods |
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 |