CN111852467B - Method and system for delineating extension range of sandstone uranium ore body - Google Patents

Method and system for delineating extension range of sandstone uranium ore body Download PDF

Info

Publication number
CN111852467B
CN111852467B CN202010736639.1A CN202010736639A CN111852467B CN 111852467 B CN111852467 B CN 111852467B CN 202010736639 A CN202010736639 A CN 202010736639A CN 111852467 B CN111852467 B CN 111852467B
Authority
CN
China
Prior art keywords
data
ore
seismic
logging
wave impedance
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
Application number
CN202010736639.1A
Other languages
Chinese (zh)
Other versions
CN111852467A (en
Inventor
李子伟
曹成寅
黄昱丞
吴曲波
郭江川
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Beijing Research Institute of Uranium Geology
Original Assignee
Beijing Research Institute of Uranium Geology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Research Institute of Uranium Geology filed Critical Beijing Research Institute of Uranium Geology
Priority to CN202010736639.1A priority Critical patent/CN111852467B/en
Publication of CN111852467A publication Critical patent/CN111852467A/en
Application granted granted Critical
Publication of CN111852467B publication Critical patent/CN111852467B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • E21B49/02Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells by mechanically taking samples of the soil
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • E21B47/12Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling
    • E21B47/14Means for transmitting measuring-signals or control signals from the well to the surface, or from the surface to the well, e.g. for logging while drilling using acoustic waves
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/02Agriculture; Fishing; Mining
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30Assessment of water resources

Abstract

The invention discloses a method and a system for delineating the extension range of sandstone uranium ore bodies. The invention belongs to the sandstone uranium mine exploration field. The method comprises the following steps: logging in an industrial hole to obtain logging data, geological stratification and lithology histograms; coring a target layer containing ores in an industrial hole, measuring the longitudinal wave velocity, the transverse wave velocity and the density of a rock sample, calculating the longitudinal wave impedance and the transverse wave impedance, and manufacturing a rock physical quantity plate for distinguishing an ore body from surrounding rocks; acquiring three-dimensional seismic data around an industrial hole, finishing data processing and interpretation, acquiring longitudinal wave impedance distribution and transverse wave impedance distribution of a target layer through seismic prestack inversion on the basis, and circling the distribution range of ore-containing rocks on the longitudinal wave impedance distribution diagram and the transverse wave impedance distribution diagram according to a rock physical quantity plate to acquire the extension range of sandstone uranium ore bodies. The method provides a technical method for delineating the extension range of the sandstone uranium ore body, is green and environment-friendly, and accelerates the exploration period of the sandstone uranium ore.

Description

Method and system for delineating extension range of sandstone uranium ore body
Technical Field
The invention relates to the technical field of sandstone uranium ore exploration, in particular to a method and a system for delineating the extension range of a sandstone uranium ore body.
Background
In sandstone uranium ore exploration, after an industrial hole with mining value is found, in order to determine the extension range of a sandstone uranium ore body, a large number of drill holes are required to be arranged to investigate the extension of the ore body, the range of the ore body is defined, and the environment can be damaged while a large amount of expenditure is invested. How to realize the delineation of the extension range of sandstone uranium ore bodies under the condition of not arranging a large number of drill holes becomes a technical problem to be solved urgently.
Disclosure of Invention
The invention aims to provide a method and a system for delineating the extension range of a sandstone uranium ore body, so as to realize the delineation of the extension range of the sandstone uranium ore body without arranging a large number of drill holes.
In order to achieve the purpose, the invention provides the following scheme:
a delineation method for the extension range of sandstone uranium ore bodies comprises the following steps:
carrying out logging and logging operations on the industrial hole to obtain logging data and logging data; the logging data comprises sound waves and densities at different depths of the industrial hole, and the logging data comprises geological stratification and lithology histograms of the industrial hole;
sampling an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples;
according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample, manufacturing a rock physical quantity plate for distinguishing an ore body from surrounding rocks;
acquiring three-dimensional seismic data of an area to be surveyed around the industrial hole through a seismograph;
performing stack processing and migration processing on the three-dimensional seismic data by using seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data comprises pure wave data and achievement data;
interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, wherein the seismic interpretation data comprise a layer position and a fracture position of a mineral-containing target layer;
performing seismic prestack inversion on a seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mineral-containing target layer;
and according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-bearing target layer, the extension range of the sandstone uranium ore body is defined by the rock physical quantity plate.
Optionally, the logging and logging operations are performed on the industrial hole, and logging data are obtained, which specifically include:
measuring sound waves and density of different depths of the industrial hole by using a logging instrument;
drilling and coring are carried out on the industrial hole to obtain a rock core;
determining geological stratification of the industrial hole according to lithology of each section of the rock core;
a lithology histogram is created from the lithology of each geological layer of the core.
Optionally, the manufacturing of the rock physical quantity plate for distinguishing the ore body from the surrounding rock according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample specifically includes:
respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method;
measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method;
respectively multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample with the density of each ore-containing sample and each surrounding rock sample to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and longitudinal wave impedance as an ordinate axis of the rock physical quantity plate;
respectively projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample to the coordinate system to obtain an intersection graph;
and (4) enclosing the area containing the ore sample on the cross plot to obtain a rock physical quantity plate for distinguishing the ore body from the surrounding rock.
Optionally, the three-dimensional seismic processing data is interpreted by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, where the seismic interpretation data includes a horizon and a fracture position of a mineral-bearing target layer, and the method specifically includes:
adopting three-dimensional seismic data interpretation software to perform synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data, and completing calibration of the horizon of the mineral-containing target layer of the seismic data;
and tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer, and determining the fracture position of the ore-containing target layer according to the homodromous characteristics of the three-dimensional seismic processing data.
Optionally, the performing seismic prestack inversion on the seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer specifically includes:
processing the three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume;
and performing seismic prestack inversion on the seismic prestack gather data body by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer.
A delineation system for sandstone uranium ore body extension, the delineation system comprising:
the logging and logging module is used for performing logging and logging operations on the industrial hole to acquire logging data and logging data; the logging data comprises sound waves and densities at different depths of the industrial hole, and the logging data comprises geological stratification and lithology histograms of the industrial hole;
the sample acquisition module is used for sampling in an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples;
the rock physical gauge plate establishing module is used for manufacturing a rock physical gauge plate for distinguishing ore bodies and surrounding rocks according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
the three-dimensional seismic data acquisition module is used for acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole through a seismograph;
the three-dimensional seismic data processing module is used for carrying out stack processing and migration processing on the three-dimensional seismic data by utilizing seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data comprises pure wave data and achievement data;
the three-dimensional seismic data interpretation module is used for interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, and the seismic interpretation data comprise a horizon and a fracture position of a mineral-containing target layer;
the three-dimensional seismic data inversion module is used for performing seismic prestack inversion on a seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mineral-containing target layer;
and the delineating module is used for delineating the extension range of the sandstone uranium ore body by utilizing the rock physical quantity plate according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-bearing target layer.
Optionally, the logging and logging module specifically includes:
the logging submodule is used for measuring sound waves and density of different depths of the industrial hole by using a logging instrument;
the drilling coring sub-module is used for performing drilling coring on the industrial hole to obtain a core;
the geological stratification determining submodule is used for determining the geological stratification of the industrial hole according to the lithology of each section of the rock core;
and the lithology histogram establishing sub-module is used for establishing the lithology histogram according to the lithology of each geological layer of the rock core.
Optionally, the rock physical quantity plate building module specifically includes:
the wave velocity measuring submodule is used for respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method;
the density measurement submodule is used for respectively measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method;
the impedance calculation submodule is used for multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by the density of each ore-containing sample and each surrounding rock sample respectively to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
the coordinate system establishing submodule is used for establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and longitudinal wave impedance as an ordinate axis of the rock physical quantity plate;
the intersection point projection submodule is used for projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample onto the coordinate system respectively to obtain an intersection graph;
and the mineral-containing rock impedance distribution range delineating sub-module is used for delineating the region of the mineral-containing sample on the cross plot to obtain a rock physical quantity plate for distinguishing the ore body from the surrounding rock.
Optionally, the three-dimensional seismic data interpretation module specifically includes:
the horizon determining submodule is used for performing synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data by adopting three-dimensional seismic data interpretation software to finish the calibration of horizons of mineral-containing target horizons of the seismic data;
and the fracture position determining submodule is used for tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer and determining the fracture position of the ore-containing target layer according to the same-direction axial characteristics of the three-dimensional seismic processing data.
Optionally, the three-dimensional seismic data inversion module specifically includes:
the data processing submodule is used for processing the three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume;
and the impedance inversion submodule is used for performing seismic prestack inversion on the seismic prestack gather data volume by utilizing the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a method and a system for delineating the extension range of sandstone uranium ore bodies. The method comprises the following steps: logging in an industrial hole to obtain logging data, geological stratification and lithology histograms; coring a target layer containing ores in an industrial hole, measuring the longitudinal wave velocity, the transverse wave velocity and the density of a rock sample, calculating the longitudinal wave impedance and the transverse wave impedance, and manufacturing a rock physical quantity plate for distinguishing an ore body from surrounding rocks; acquiring three-dimensional seismic data around an industrial hole, finishing data processing and interpretation, acquiring longitudinal wave impedance distribution and transverse wave impedance distribution of a target layer through seismic prestack inversion on the basis, and circling the distribution range of ore-containing rocks on the longitudinal wave impedance distribution diagram and the transverse wave impedance distribution diagram according to a rock physical quantity plate to acquire the extension range of sandstone uranium ore bodies. The method provides a technical method for delineating the extension range of the sandstone uranium ore body, is green and environment-friendly, and accelerates the exploration period of the sandstone uranium ore.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings needed to be used in the embodiments will be briefly described below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings without inventive exercise.
Fig. 1 is a flow chart of a method for delineating the extension range of a sandstone uranium ore body provided by the invention.
Detailed Description
The invention aims to provide a method and a system for delineating the extension range of a sandstone uranium ore body, so as to realize the delineation of the extension range of the sandstone uranium ore body without arranging a large number of drill holes.
In order to make the aforementioned objects, features and advantages of the present invention comprehensible, embodiments accompanied with figures are described in further detail below.
In order to solve the technical problems, the invention provides that after an industrial hole is found, ground exploration work, namely three-dimensional seismic exploration, is deployed at the periphery to replace drilling and quickly define the extension range of sandstone uranium ore bodies, so that the purpose is realized, the damage to the earth surface and the underground is reduced, and the requirement of green exploration is met.
As shown in fig. 1, the invention provides a method for delineating the extension range of sandstone uranium ore body, which comprises the following steps:
101, performing logging and logging operations on an industrial hole to obtain logging data and logging data; the well logging data includes sound waves and densities at different depths of the industrial hole, and the well logging data includes geological stratification and lithology histograms of the industrial hole.
Step 101 is to perform logging and logging operations on the industrial hole, acquiring logging data and logging data, specifically comprising: measuring sound waves and density at different depths of the industrial hole by using a logging instrument; drilling and coring the industrial hole to obtain a core; determining geological stratification of the industrial hole according to lithology of each section of the rock core; a lithology histogram is created from the lithology of each geological layer of the core.
Specifically, the logging data in step 101 is collected in the well through the logging tool, and the logging data is collected by sliding a parameter probe of the logging tool downhole, and measuring data every 0.05 m.
The logging data in step 101 is that geologists take cores according to the drilling, determine geological stratification according to regional stratum characteristics by observing and describing the lithology of each section of core, and draw a lithology histogram by using software.
And 102, sampling in an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples.
The sampling in the step 102 is well core sampling, the sampling depth range covers a target layer containing ore, the sampling depth interval is uniform, the ore body is subjected to encrypted sampling, and the number of the ore-containing samples and the number of the surrounding rock samples are not less than 10 respectively.
And 103, manufacturing a rock physical quantity plate for distinguishing ore bodies and surrounding rocks according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample.
103, according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample, manufacturing a rock physical quantity plate for distinguishing an ore body from surrounding rocks, which specifically comprises the following steps: respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method; measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method; respectively multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample with the density of each ore-containing sample and each surrounding rock sample to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample; establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and taking longitudinal wave impedance as an ordinate axis of the rock physical quantity plate; respectively projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample to the coordinate system to obtain an intersection graph; and (4) enclosing the area containing the ore sample on the cross plot to obtain a rock physical quantity plate for distinguishing the ore body from the surrounding rock.
Specifically, the wave velocity measurement in step 103 is performed by an elastic parameter measurement method, an ultrasonic pulse transmission method is used to measure the longitudinal wave velocity and the transverse wave velocity of the sample (the ore-containing sample and the surrounding rock sample), and an integral measurement method is used to measure the density of the sample (the ore-containing sample and the surrounding rock sample).
The calculation method of the longitudinal wave impedance and the transverse wave impedance in step 103 is as follows: longitudinal wave impedance = longitudinal wave velocity x density; shear wave impedance = shear wave velocity × density.
The manufacturing of the rock physical quantity plate in the step 103 is realized by intersecting longitudinal wave impedance and transverse wave impedance, the abscissa axis of the rock physical quantity plate is transverse wave impedance, the ordinate axis is longitudinal wave impedance, the impedance calculation results of the ore-containing sample and the surrounding rock sample are projected onto an intersection graph, and the areas of the ore-containing sample and the surrounding rock sample are respectively circled, so that the manufacturing of the rock physical quantity plate is completed.
And 104, acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole by using a seismograph.
Step 104, acquiring three-dimensional seismic data, namely acquiring field actual measurement data of the earthquake by using a seismometer, wherein a bin of an observation system is less than or equal to 10m multiplied by 10m, the covering times are not less than 48 times, the maximum offset distance is determined according to the exploration depth, and the calculation mode is as follows: the maximum offset is more than or equal to the exploration depth multiplied by 1.5.
And step 104, acquiring the three-dimensional seismic data, wherein the main frequency of the detector is not higher than 10Hz, and the seismic source can be a vibroseis or a dynamite.
105, stacking and shifting the three-dimensional seismic data by using seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and production data.
And 106, interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, wherein the seismic interpretation data comprise the position and the fracture position of the ore-containing target layer.
106, interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, wherein the seismic interpretation data comprise a layer position and a fracture position of a mineral-containing target layer, and the method specifically comprises the following steps: adopting three-dimensional seismic data interpretation software to perform synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data, and completing calibration of the horizon of the mineral-containing target layer of the seismic data; and tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer, and determining the fracture position of the ore-containing target layer according to the homodromous characteristics of the three-dimensional seismic processing data.
And 106, interpreting the seismic data, specifically, by using seismic data interpretation software, synthesizing the logging data into a seismic record, completing calibration of a target layer position of the seismic data, further performing layer position tracking on the seismic data, identifying a fracture position according to the same-direction axis characteristics of the seismic data, and finally interpreting the layer position and the fracture of the target layer on a seismic data body.
And 107, performing seismic prestack inversion on the seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-containing target layer.
Step 107, performing seismic prestack inversion on the seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer, specifically comprising: processing the three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume and obtain a seismic prestack gather data volume; and performing seismic prestack inversion on the seismic prestack gather data volume by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer.
The seismic prestack inversion in step 107 is to perform seismic prestack inversion on the seismic prestack gather data volume by using seismic prestack inversion software to obtain the longitudinal wave impedance parameter distribution and the transverse wave impedance parameter distribution of the target layer. The method comprises the steps of firstly, carrying out synthetic seismic recording operation on logging data, a seismic prestack gather data volume and seismic interpretation data, completing well seismic calibration of the seismic data, obtaining seismic data after well seismic calibration, then extracting seismic wavelets from the seismic data after well seismic calibration, then combining the seismic data and the logging data, establishing an inversion initial model, carrying out seismic inversion, obtaining a longitudinal wave impedance parameter and a transverse wave impedance parameter of a target layer through inversion calculation, specifically, setting the longitudinal wave impedance parameter and the transverse wave impedance parameter of the initial target layer, establishing an inversion initial model by using the longitudinal wave impedance parameter and the transverse wave impedance parameter of the initial target layer, carrying out convolution operation on the seismic and the inversion initial model to obtain forward seismic data, comparing the forward seismic data with the seismic data obtained through measurement, outputting the longitudinal wave impedance parameter and the transverse wave impedance parameter when the difference between the longitudinal wave impedance parameter and the transverse wave impedance parameter is within an allowable threshold range, adjusting the longitudinal wave impedance parameter and the transverse wave impedance parameter until the forward seismic data is reestablished, and obtaining the difference value of the forward seismic data within the allowable threshold range.
And 108, utilizing the rock physical quantity plate to circle the extension range of the sandstone uranium ore body according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-containing target layer.
In step 108, on the results of the longitudinal wave impedance and the transverse wave impedance of the target layer, the distribution range of the sandstone uranium ore body can be defined by using the rock physical quantity plate.
The invention also provides a delineation system for the extension range of sandstone uranium ore bodies, which comprises the following components:
the logging and logging module is used for performing logging and logging operations on the industrial hole to acquire logging data and logging data; the well logging data includes sound waves and densities at different depths of the industrial hole, and the well logging data includes geological stratification and lithology histograms of the industrial hole.
The logging and logging module specifically comprises: the logging submodule is used for measuring sound waves and density of different depths of the industrial hole by using a logging instrument; the drilling coring sub-module is used for performing drilling coring on the industrial hole to obtain a core; the geological stratification determining submodule is used for determining the geological stratification of the industrial hole according to the lithology of each section of the rock core; and the lithology histogram establishing sub-module is used for establishing the lithology histogram according to the lithology of each geological layer of the rock core.
The sample acquisition module is used for sampling in an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples;
and the rock physical gauge plate establishing module is used for manufacturing the rock physical gauge plate for distinguishing the ore body and the surrounding rock according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample.
The rock physical quantity plate establishing module specifically comprises: the wave velocity measuring submodule is used for respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method; the density measurement submodule is used for respectively measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method; the impedance calculation submodule is used for multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by the density of each ore-containing sample and each surrounding rock sample respectively to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample; the coordinate system establishing submodule is used for establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and longitudinal wave impedance as an ordinate axis of the rock physical quantity plate; the intersection point projection submodule is used for projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample to the coordinate system respectively to obtain an intersection graph; and the mineral-containing rock impedance distribution range delineating sub-module is used for delineating the region of the mineral-containing sample on the cross plot to obtain a rock physical quantity plate for distinguishing the ore body from the surrounding rock.
The three-dimensional seismic data acquisition module is used for acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole through a seismograph;
the three-dimensional seismic data processing module is used for carrying out stack processing and migration processing on the three-dimensional seismic data by utilizing seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data comprises pure wave data and achievement data;
and the three-dimensional seismic data interpretation module is used for interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, wherein the seismic interpretation data comprise the horizon and the fracture position of the ore-containing target layer.
The three-dimensional seismic data interpretation module specifically comprises: the horizon determining submodule is used for performing synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data by adopting three-dimensional seismic data interpretation software to finish the calibration of horizons of mineral-containing target horizons of the seismic data; and the fracture position determining submodule is used for tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer and determining the fracture position of the ore-containing target layer according to the same-direction axial characteristics of the three-dimensional seismic processing data.
And the three-dimensional seismic data inversion module is used for performing seismic prestack inversion on a seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mineral-containing target layer.
The three-dimensional seismic data inversion module specifically comprises: the data processing submodule is used for processing three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume and obtain a seismic prestack gather data volume; and the impedance inversion submodule is used for performing seismic prestack inversion on the seismic prestack gather data volume by utilizing the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mining target layer.
And the delineating module is used for delineating the extension range of the sandstone uranium ore body by utilizing the rock physical quantity plate according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-bearing target layer.
According to the specific embodiment provided by the invention, the invention discloses the following technical effects:
the invention discloses a method and a system for delineating the extension range of sandstone uranium ore bodies. The method comprises the following steps: logging in an industrial hole to obtain logging data, geological stratification and lithology histograms; coring a target layer containing ores in an industrial hole, measuring the longitudinal wave velocity, the transverse wave velocity and the density of a rock sample, calculating the longitudinal wave impedance and the transverse wave impedance, and manufacturing a rock physical quantity plate for distinguishing an ore body from surrounding rocks; acquiring three-dimensional seismic data around an industrial hole, finishing data processing and interpretation, acquiring longitudinal wave impedance distribution and transverse wave impedance distribution of a target layer through seismic prestack inversion on the basis, and circling the distribution range of ore-containing rocks on the longitudinal wave impedance distribution diagram and the transverse wave impedance distribution diagram according to a rock physical quantity plate to acquire the extension range of sandstone uranium ore bodies. The method provides a technical method for delineating the extension range of the sandstone uranium ore body, is green and environment-friendly, and accelerates the exploration period of the sandstone uranium ore.
The technical means of logging, rock sample measurement, three-dimensional seismic data acquisition, processing, interpretation, inversion and the like are all general technical means in the field, and various changes can be made on the premise of not departing from the purpose of the invention. As long as three-dimensional seismic data, logging data and rock sample elastic data are comprehensively utilized, the extension range of the sandstone uranium ore body can be defined by adopting the prior art commonly used in the field according to the scheme of the invention.
The embodiments in the present description are described in a progressive manner, each embodiment focuses on differences from other embodiments, and the same and similar parts among the embodiments are referred to each other.
The principle and the implementation manner of the present invention are explained by applying specific examples, the above description of the embodiments is only used to help understanding the method of the present invention and the core idea thereof, the described embodiments are only a part of the embodiments of the present invention, not all embodiments, and all other embodiments obtained by a person of ordinary skill in the art based on the embodiments of the present invention without creative efforts belong to the protection scope of the present invention.

Claims (8)

1. The method for delineating the extension range of the sandstone uranium ore body is characterized by comprising the following steps of:
carrying out logging and logging operations on the industrial hole to obtain logging data and logging data; the logging data comprises sound waves and densities at different depths of the industrial hole, and the logging data comprises geological stratification and lithology histograms of the industrial hole;
sampling in an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples;
according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample, a rock physical quantity plate for distinguishing an ore body from surrounding rocks is manufactured, and the method specifically comprises the following steps:
respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method;
measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method;
respectively multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample with the density of each ore-containing sample and each surrounding rock sample to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and longitudinal wave impedance as an ordinate axis of the rock physical quantity plate;
respectively projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample to the coordinate system to obtain an intersection graph;
a region containing an ore sample and a region of a surrounding rock sample are circled on the cross plot, and a rock physical quantity plate for distinguishing an ore body from the surrounding rock is obtained;
acquiring three-dimensional seismic data of a region to be surveyed around the industrial hole by using a seismograph;
performing stack processing and migration processing on the three-dimensional seismic data by using seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data comprises pure wave data and achievement data;
interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data; the seismic interpretation data comprises horizons and fracture positions of the ore-bearing target layer;
performing seismic prestack inversion on a seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mineral-containing target layer;
and according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-bearing target layer, the extension range of the sandstone uranium ore body is defined by the rock physical quantity plate.
2. The method for delineating the extension range of sandstone uranium ore bodies according to claim 1, wherein the operations of logging and logging are performed on an industrial hole to obtain logging data and logging data, and specifically comprises:
measuring sound waves and density at different depths of the industrial hole by using a logging instrument;
drilling and coring are carried out on the industrial hole to obtain a rock core;
determining geological stratification of the industrial hole according to lithology of each section of the rock core;
a lithology histogram is created from the lithology of each geological layer of the core.
3. The method for delineating the extension range of sandstone uranium ore bodies according to claim 1, wherein the three-dimensional seismic processing data is interpreted by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, wherein the seismic interpretation data comprise horizons and fracture positions of ore-bearing target layers, and specifically comprises:
adopting three-dimensional seismic data interpretation software to perform synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data, and completing calibration of the horizon of the mineral-containing target layer of the seismic data;
and tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer, and determining the fracture position of the ore-containing target layer according to the homodromous characteristics of the three-dimensional seismic processing data.
4. The method for delineating the extension range of sandstone uranium ore bodies according to claim 1, wherein the seismic prestack inversion is performed on a seismic prestack gather data body of three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain a longitudinal wave impedance distribution and a transverse wave impedance distribution of an ore-bearing target layer, and specifically comprises:
processing the three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume;
and performing seismic prestack inversion on the seismic prestack gather data volume by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer.
5. A delineation system for sandstone uranium ore body extension, the delineation system comprising:
the logging and logging module is used for performing logging and logging operations on the industrial hole to acquire logging data and logging data; the logging data comprises sound waves and densities at different depths of the industrial hole, and the logging data comprises geological stratification and lithology histograms of the industrial hole;
the sample acquisition module is used for sampling in an ore-containing target layer of the industrial hole to obtain a plurality of ore-containing samples and a plurality of surrounding rock samples;
the rock physical gauge plate establishing module is used for manufacturing a rock physical gauge plate for distinguishing ore bodies and surrounding rocks according to the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
the rock physical quantity plate establishing module specifically comprises:
the wave velocity measuring submodule is used for respectively measuring the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by adopting an ultrasonic pulse transmission method;
the density measurement sub-module is used for respectively measuring the density of each ore-containing sample and each surrounding rock sample by adopting a volume method;
the impedance calculation submodule is used for multiplying the longitudinal wave velocity and the transverse wave velocity of each ore-containing sample and each surrounding rock sample by the density of each ore-containing sample and each surrounding rock sample respectively to obtain the longitudinal wave impedance and the transverse wave impedance of each ore-containing sample and each surrounding rock sample;
the coordinate system establishing submodule is used for establishing a coordinate system by taking transverse wave impedance as an abscissa axis of the rock physical quantity plate and longitudinal wave impedance as an ordinate axis of the rock physical quantity plate;
the intersection point projection submodule is used for projecting intersection points of longitudinal wave impedance and transverse wave impedance of each ore-containing sample and each surrounding rock sample onto the coordinate system respectively to obtain an intersection graph;
the mineral-containing rock impedance distribution range delineating sub-module is used for delineating a region of a mineral-containing sample on the cross plot to obtain a rock physical quantity plate for distinguishing an ore body and surrounding rocks;
the three-dimensional seismic data acquisition module is used for acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole through a seismograph;
the three-dimensional seismic data processing module is used for carrying out stack processing and migration processing on the three-dimensional seismic data by utilizing seismic data processing software to obtain three-dimensional seismic processing data; the seismic processing data comprises pure wave data and result data;
the three-dimensional seismic data interpretation module is used for interpreting the three-dimensional seismic processing data by using the logging data and the logging data to obtain seismic interpretation data of the three-dimensional seismic data, and the seismic interpretation data comprise a horizon and a fracture position of a mineral-containing target layer;
the three-dimensional seismic data inversion module is used for performing seismic prestack inversion on a seismic prestack gather data volume of the three-dimensional seismic data by using the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of a mineral-containing target layer;
and the delineating module is used for delineating the extension range of the sandstone uranium ore body by utilizing the rock physical quantity plate according to the longitudinal wave impedance distribution and the transverse wave impedance distribution of the ore-bearing target layer.
6. The delineation system of sandstone uranium ore body extension according to claim 5, wherein the logging and logging module specifically comprises:
the logging submodule is used for measuring sound waves and density of different depths of the industrial hole by using a logging instrument;
the drilling coring sub-module is used for performing drilling coring on the industrial hole to obtain a core;
the geological stratification determining submodule is used for determining the geological stratification of the industrial hole according to the lithology of each section of the rock core;
and the lithology histogram establishing sub-module is used for establishing the lithology histogram according to the lithology of each geological layer of the rock core.
7. The delineation system of sandstone uranium ore body extension according to claim 5, wherein the three-dimensional seismic data interpretation module specifically comprises:
the horizon determining submodule is used for performing synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data by adopting three-dimensional seismic data interpretation software to finish the calibration of horizons of mineral-containing target horizons of the seismic data;
and the fracture position determining submodule is used for tracking the three-dimensional seismic processing data of the layer of the ore-containing target layer and determining the fracture position of the ore-containing target layer according to the same-direction axial characteristics of the three-dimensional seismic processing data.
8. The delineation system of sandstone uranium ore body extension according to claim 5, wherein the three-dimensional seismic data inversion module specifically comprises:
the data processing submodule is used for processing the three-dimensional seismic data by using three-dimensional seismic data processing software to obtain a seismic prestack gather data volume;
and the impedance inversion submodule is used for performing seismic prestack inversion on the seismic prestack gather data volume by utilizing the logging data, the logging data and the seismic interpretation data to obtain longitudinal wave impedance distribution and transverse wave impedance distribution of the ore-bearing target layer.
CN202010736639.1A 2020-07-28 2020-07-28 Method and system for delineating extension range of sandstone uranium ore body Active CN111852467B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010736639.1A CN111852467B (en) 2020-07-28 2020-07-28 Method and system for delineating extension range of sandstone uranium ore body

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010736639.1A CN111852467B (en) 2020-07-28 2020-07-28 Method and system for delineating extension range of sandstone uranium ore body

Publications (2)

Publication Number Publication Date
CN111852467A CN111852467A (en) 2020-10-30
CN111852467B true CN111852467B (en) 2023-04-07

Family

ID=72948852

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010736639.1A Active CN111852467B (en) 2020-07-28 2020-07-28 Method and system for delineating extension range of sandstone uranium ore body

Country Status (1)

Country Link
CN (1) CN111852467B (en)

Families Citing this family (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117250658B (en) * 2023-11-17 2024-02-09 核工业北京地质研究院 Method for creating seismic dataset of investigation region
CN117557401B (en) * 2024-01-12 2024-04-02 东华理工大学南昌校区 Geological big data-based uranium ore prospecting target area intelligent demarcating method

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977618A (en) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 Method for evaluating shale gas reservoir and finding dessert area
CN106054248A (en) * 2016-07-15 2016-10-26 河海大学 Earthquake rock physical inversion method based on large area tight reservoir
CN106368691A (en) * 2015-07-24 2017-02-01 中国石油化工股份有限公司 Method for predicting three-dimensional abnormal pore pressure based on rock physical seismic information
CN108897041A (en) * 2018-08-16 2018-11-27 中国石油天然气股份有限公司 A kind of prediction technique and device of uranium ore enrichment region

Family Cites Families (10)

* Cited by examiner, † Cited by third party
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
CN103364822B (en) * 2012-04-05 2016-08-24 核工业北京地质研究院 A kind of layered approach being applicable to sedimentary basin shallow layer stratum
CN103954995B (en) * 2014-04-22 2016-12-07 核工业北京地质研究院 A kind of sand body recognition methods in exploration of sandstone type uranium deposits
CN105093306B (en) * 2014-05-15 2018-03-09 中国石油化工股份有限公司 Reservoir is explained and thickness acquiring method automatically in a kind of geophysical exploration
CN105044770A (en) * 2015-07-06 2015-11-11 成都理工大学 Compact glutenite gas reservoir quantificational prediction method
CN105403909A (en) * 2015-08-13 2016-03-16 核工业北京地质研究院 Method for detecting paleochannel type uranium mine
CN108802812B (en) * 2017-04-28 2020-02-14 中国石油天然气股份有限公司 Well-seismic fusion stratum lithology inversion method
CN107576982B (en) * 2017-09-01 2019-05-17 核工业北京地质研究院 A kind of sandstone-type uranium mineralization with respect Comprehensive Seismic Prediction method
CN109738947B (en) * 2018-12-12 2020-11-20 核工业北京地质研究院 Physicochemical exploration combination method for ore exploration scenic spot of confined sandstone type uranium deposit
CN111158052B (en) * 2020-01-07 2021-07-13 吉林大学 Method for delineating distribution range beside uranium-bearing sand body well in inversion well based on three-dimensional reservoir

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104977618A (en) * 2014-04-09 2015-10-14 中国石油集团东方地球物理勘探有限责任公司 Method for evaluating shale gas reservoir and finding dessert area
CN106368691A (en) * 2015-07-24 2017-02-01 中国石油化工股份有限公司 Method for predicting three-dimensional abnormal pore pressure based on rock physical seismic information
CN106054248A (en) * 2016-07-15 2016-10-26 河海大学 Earthquake rock physical inversion method based on large area tight reservoir
CN108897041A (en) * 2018-08-16 2018-11-27 中国石油天然气股份有限公司 A kind of prediction technique and device of uranium ore enrichment region

Also Published As

Publication number Publication date
CN111852467A (en) 2020-10-30

Similar Documents

Publication Publication Date Title
Li et al. An overview of ahead geological prospecting in tunneling
CN104977618B (en) A kind of method evaluated shale gas reservoir and find dessert area
CN104853822A (en) Method for evaluating shale gas reservoir and searching sweet spot region
CN105629325B (en) Foreland basin alluvial fan meticulous depiction and Forecasting Methodology
CN105510993A (en) Foreland basin deep buried and compressed type complex gypsum-salt rock identification and distribution prediction method
US20190094397A1 (en) Surface detection and location of microseismic events and earthquakes without the use of a velocity model
CN111852467B (en) Method and system for delineating extension range of sandstone uranium ore body
CN105629308A (en) Phase control heterogeneous mechanical parameter crustal stress prediction method
CN105572745A (en) Seismic prospecting method of three-component slot waves under coalmine well
CN107132573A (en) A kind of method that application wavelet decomposition reconfiguration technique recognizes the lower lithological pool of strong impedance shielding
EP3959544A1 (en) Joint sensor orientation and velocity model calibration
Agoha et al. Weathered layer determination and its effects on engineering structures: case study of parts of the Niger delta
Anbazhagan et al. Evaluation of dynamic properties and ground profiles using MASW: correlation between Vs and N60
Lacazette et al. Geomechanical and flow simulation of hydraulic fractures using high-resolution passive seismic images
CN111983721A (en) Sandstone uranium ore mud-sand-mud geological structure identification method and system
CN110454155B (en) Method for determining age of quaternary stratum by using susceptibility logging method
Bergman et al. High-resolution reflection seismic imaging of the upper crust at Laxemar, southeastern Sweden
Yılmaz et al. Shallow seismic and geotechnical site surveys at the Turkish national grid for strong-motion seismograph stations
CN112946753B (en) Near-surface stratum structure analysis system and method
Shaibu et al. Uses of seismic refraction tomography in investigation of groundwater potential in Kaduna millennium city, Kaduna Nigeria
Lacazette et al. Passive Seismic Methods for Unconventional Resource Development
US11474272B2 (en) Methods and systems for identifying and plugging subterranean conduits
CN112305594B (en) Oil gas distribution determining method and system for heterogeneous reservoir
Frei Hybrid seismic surveying for detailed characterization of the shallow and intermediate depths subsurface
CN114594529A (en) Method for predicting buried hill reservoir based on azimuth seismic attribute change rate

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