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

Abstract

The invention discloses a method and system for delineating the extension range of a sandstone uranium ore body. The invention belongs to the field of sandstone uranium ore exploration. The method comprises the following steps: performing well logging and mud logging in an industrial hole to obtain well logging data, geological stratification and lithology histogram; taking cores from the ore-bearing target layer of the industrial hole to measure the compressional wave velocity, S-wave velocity and density, calculate P-wave impedance and S-wave impedance, and make petrophysical quantity boards for distinguishing ore bodies and surrounding rocks; collect 3D seismic data around industrial holes, complete data processing and interpretation, and use seismic pre-stack reflection to complete data processing and interpretation. Obtain the P-wave impedance distribution and S-wave impedance distribution of the target layer, delineate the distribution range of ore-bearing rocks on the P-wave impedance distribution map and S-wave impedance distribution map according to the rock physical quantity plate, and obtain the extension range of the sandstone uranium ore body. This method provides a technical method for the delineation of the extension range of sandstone uranium deposits, which is environmentally friendly and speeds up the exploration cycle of sandstone uranium deposits.

Description

A method and system for delineating the extension range of a sandstone uranium ore body

technical field

The invention relates to the technical field of sandstone uranium ore exploration, in particular to a method and system for delineating the extension range of a sandstone uranium ore body.

Background technique

In the exploration of sandstone uranium ore, after discovering an industrial hole with mining value, in order to determine the extension range of the sandstone uranium ore body, it is necessary to lay out a large number of drill holes to investigate the extension of the ore body, delineate the range of the ore body, and invest in While costing a lot of money, it will also cause damage to the environment. How to delineate the extension range of the sandstone uranium ore body without laying a large number of drill holes has become an urgent technical problem to be solved.

Contents of the invention

The object of the present invention is to provide a method and system for delineating the extension range of sandstone uranium ore bodies, so as to realize the delineation of the extension range of sandstone uranium ore bodies without laying a large number of drill holes.

To achieve the above object, the present invention provides the following scheme:

A method for delineating the extension range of a sandstone uranium ore body, said delineating method comprising the following steps:

Perform well logging and mud logging operations on industrial holes to obtain well logging data and mud logging data; the logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data includes geological layers and rock formations of industrial holes sex histogram;

Sampling is carried out in the ore-bearing target layer of the industrial hole to obtain multiple ore-bearing samples and multiple surrounding rock samples;

Based on the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample, a petrophysical quantity board for distinguishing ore bodies and surrounding rocks is made;

Acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole through a seismograph;

Using seismic data processing software, performing superposition processing and migration processing on the three-dimensional seismic data to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and result data;

Interpreting the three-dimensional seismic processing data by using the well logging data and the mud logging data to obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data including layers and fracture positions of ore-bearing target layers;

Using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume of the 3D seismic data, and obtain the compressional wave impedance distribution and shear wave impedance of the ore-bearing target layer distributed;

According to the distribution of longitudinal wave impedance and shear wave impedance of the ore-bearing target layer, the extension range of the sandstone uranium ore body is delineated by using the petrophysical quantity plate.

Optionally, performing logging and mud logging operations on industrial holes to obtain logging data and mud logging data specifically includes:

Use logging tools to measure sound waves and densities at different depths in industrial boreholes;

Drill and core industrial holes to obtain cores;

Determine the geological stratification of industrial holes according to the lithology of each section of the core;

A lithology histogram is built according to the lithology of each geological layer of the core.

Optionally, according to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample, making a petrophysical quantity board for distinguishing ore bodies and surrounding rocks specifically includes:

The longitudinal wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample were measured by ultrasonic pulse transmission method;

The density of each ore-bearing sample and each surrounding rock sample is measured separately by the volumetric method;

Multiply the compressional wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by the density of each ore-bearing sample and each surrounding rock sample to obtain the compressional wave of each ore-bearing sample and each surrounding rock sample impedance and shear wave impedance;

The coordinate system is established by taking the shear wave impedance as the abscissa axis of the petrophysical quantity plate, and taking the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate;

projecting the intersection points of the longitudinal wave impedance and the shear wave impedance of each ore-bearing sample and each surrounding rock sample onto the coordinate system to obtain an intersection map;

The ore-bearing sample area is circled on the intersection map to obtain a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks.

Optionally, using the well logging data and the mud logging data to interpret the three-dimensional seismic processing data to obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data includes the horizon of the ore-bearing target layer and fracture locations, including:

Using three-dimensional seismic data interpretation software, performing a synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data, and completing the calibration of the horizon of the ore-bearing target layer of the seismic data;

The three-dimensional seismic processing data of the layer of the ore-bearing target layer is tracked, and the fracture position of the ore-bearing target layer is determined according to the synaxial characteristics of the three-dimensional seismic processing data.

Optionally, using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume of the 3D seismic data to obtain the ore-bearing target layer P-wave impedance distribution and shear-wave impedance distribution, specifically including:

Use the 3D seismic data processing software to process the 3D seismic data to obtain the seismic pre-stack gather data volume;

Using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume to obtain the distribution of compressional wave impedance and shear wave impedance of the ore-bearing target layer.

A delineation system for the extended range of a sandstone uranium ore body, the delineation system comprising:

The well logging and mud logging module is used to perform well logging and mud logging operations on industrial holes, and obtain well logging data and mud logging data; the well logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data Geological stratification and lithology histograms including industrial boreholes;

The sample acquisition module is used to take samples in the ore-bearing target layer of the industrial hole, and obtain multiple ore-bearing samples and multiple surrounding rock samples;

The petrophysical quantity plate building module is used to make a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks according to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample;

The three-dimensional seismic data acquisition module is used to collect the three-dimensional seismic data of the area to be investigated around the industrial hole through the seismograph;

The three-dimensional seismic data processing module is used to use seismic data processing software to perform superposition processing and migration processing on the three-dimensional seismic data to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and result data;

A three-dimensional seismic data interpretation module, configured to interpret the three-dimensional seismic processing data by using the well logging data and the mud logging data, and obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data including the ore-bearing target layer horizons and fracture locations;

The three-dimensional seismic data inversion module is used to perform seismic pre-stack inversion on the seismic pre-stack gather data volume of the three-dimensional seismic data by using the logging data, the mud logging data and the seismic interpretation data to obtain ore-bearing P-wave impedance distribution and shear-wave impedance distribution of the target layer;

The delineation module is used to delineate the extension range of the sandstone uranium ore body by using the petrophysical quantity board according to the distribution of the longitudinal wave impedance and the shear wave impedance of the ore-bearing target layer.

Optionally, the well logging and mud logging modules specifically include:

The logging sub-module is used to measure the acoustic wave and density of industrial holes at different depths with the logging instrument;

The drilling and coring sub-module is used for drilling and coring industrial holes to obtain cores;

The sub-module for determining geological stratification is used to determine the geological stratification of industrial holes according to the lithology of each section of the core;

The lithology histogram establishment sub-module is used for establishing the lithology histogram according to the lithology of each geological layer of the core.

Optionally, the petrophysical quantity board building module specifically includes:

The wave velocity measurement sub-module is used to measure the longitudinal wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by ultrasonic pulse transmission method;

The density measurement sub-module is used to measure the density of each ore-bearing sample and each surrounding rock sample respectively by the volumetric method;

The impedance calculation sub-module is used to multiply the compressional wave velocity and the shear wave velocity of each ore-bearing sample and each surrounding rock sample by the density of each ore-bearing sample and each surrounding rock sample to obtain each ore-bearing sample and each The compressional wave impedance and shear wave impedance of each surrounding rock sample;

The coordinate system establishment sub-module is used to establish a coordinate system with the shear wave impedance as the abscissa axis of the petrophysical quantity plate and the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate;

The intersection point projection sub-module is used to project the intersection point of the longitudinal wave impedance and the shear wave impedance of each ore-bearing sample and each surrounding rock sample onto the coordinate system to obtain an intersection map;

The ore-bearing rock impedance distribution range delineation sub-module is used to delimit the ore-bearing sample area on the intersection map, and obtain a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks.

Optionally, the 3D seismic data interpretation module specifically includes:

The horizon determination sub-module is used to use the three-dimensional seismic data interpretation software to perform synthetic seismic recording operation on the three-dimensional seismic processing data and the well logging data, and complete the calibration of the horizon of the ore-bearing target layer of the seismic data;

The sub-module for determining the fracture position is used to track the 3D seismic processing data of the layer of the ore-bearing target layer, and determine the fracture position of the ore-bearing target layer according to the concentric axis characteristics of the 3D seismic processing data.

Optionally, the 3D seismic data inversion module specifically includes:

The data processing sub-module is used to process the three-dimensional seismic data by using the three-dimensional seismic data processing software to obtain the seismic pre-stack gather data body;

The impedance inversion sub-module is used to perform seismic prestack inversion on the seismic prestack gather data volume by using the well logging data, the mud logging data and the seismic interpretation data, and obtain the ore-bearing target layer P-wave impedance distribution and shear-wave impedance distribution.

According to the specific embodiments provided by the invention, the invention discloses the following technical effects:

The invention discloses a method and system for delineating the extension range of a sandstone uranium ore body. The method comprises the following steps: performing well logging and mud logging in an industrial hole to obtain well logging data, geological stratification and lithology histogram; taking cores from the ore-bearing target layer of the industrial hole to measure the compressional wave velocity, S-wave velocity and density, calculate P-wave impedance and S-wave impedance, and make petrophysical quantity boards for distinguishing ore bodies and surrounding rocks; collect 3D seismic data around industrial holes, complete data processing and interpretation, and use seismic pre-stack reflection to complete data processing and interpretation. Obtain the P-wave impedance distribution and S-wave impedance distribution of the target layer, delineate the distribution range of ore-bearing rocks on the P-wave impedance distribution map and S-wave impedance distribution map according to the rock physical quantity plate, and obtain the extension range of the sandstone uranium ore body. This method provides a technical method for the delineation of the extension range of sandstone uranium deposits, which is environmentally friendly and speeds up the exploration cycle of sandstone uranium deposits.

Description of drawings

In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the accompanying drawings required in the embodiments. Obviously, the accompanying drawings in the following description are only some of the present invention. Embodiments, for those of ordinary skill in the art, other drawings can also be obtained according to these drawings without paying creative labor.

Fig. 1 is a flow chart of a method for delineating the extension range of a sandstone uranium ore body provided by the present invention.

Detailed ways

The object of the present invention is to provide a method and system for delineating the extension range of sandstone uranium ore bodies, so as to realize the delineation of the extension range of sandstone uranium ore bodies without laying a large number of drill holes.

In order to make the above objects, features and advantages of the present invention more comprehensible, the invention will be further described in detail below in conjunction with the accompanying drawings and specific embodiments.

In order to solve the above technical problems, the present invention proposes that after an industrial hole is discovered, ground survey work, i.e. 3D seismic survey, is deployed in the surrounding area to replace drilling and quickly delineate the extension range of the sandstone uranium ore body. and underground destruction, meeting the requirements of green exploration.

As shown in Figure 1, the present invention provides a kind of delineation method of sandstone uranium ore body extension range, and described delineation method comprises the following steps:

Step 101, perform logging and mud logging operations on industrial holes, and obtain logging data and mud logging data; the logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data includes geological analysis of industrial holes. Layer and lithology histograms.

In step 101, performing well logging and mud logging operations on industrial holes to obtain well logging data and mud logging data, specifically including: using a well logging instrument to measure sound waves and densities at different depths of industrial holes; drilling and coring industrial holes, Obtain the core; determine the geological stratification of the industrial hole according to the lithology of each section of the core; establish a lithology column diagram according to the lithology of each geological layer of the core.

Specifically, the logging data in step 101 is collected in the well by the logging tool, and the logging data acquisition process is to use the parameter probe of the logging tool to slide downhole, and measure a data every 0.05m.

The mud logging data in the step 101 is obtained by geologists based on drilling and coring, describing the lithology of each section of core through observation, determining geological stratification according to regional stratum characteristics, and using software to draw a lithology columnar diagram.

Step 102: Sampling is carried out in the ore-bearing target layer of the industrial hole to obtain a plurality of ore-bearing samples and a plurality of surrounding rock samples.

The sampling in step 102 is drilling core sampling, the sampling depth range covers the ore-bearing target layer, the sampling depth interval is uniform, and the ore body is densely sampled, and the ore-bearing samples and surrounding rock samples are not less than 10 pieces respectively.

Step 103, according to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample, make a petrophysical quantity board for distinguishing ore bodies and surrounding rocks.

In step 103, according to the longitudinal wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample, making a petrophysical quantity board for distinguishing the ore body and surrounding rock, specifically includes: using the ultrasonic pulse transmission method to measure each containing The compressional-wave velocity and shear-wave velocity of each ore-bearing sample and each surrounding-rock sample were measured; the density of each ore-bearing sample and each surrounding-rock sample was measured by the volumetric method; Velocity and shear wave velocity are multiplied by the density of each ore-bearing sample and each surrounding rock sample to obtain the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample; take the shear wave impedance as the abscissa of the petrophysical quantity plate Axis, taking the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate, a coordinate system is established; the intersection points of the longitudinal wave impedance and the shear wave impedance of each ore-bearing sample and each surrounding rock sample are respectively projected onto the coordinate system to obtain the intersection Figure; Circle the area containing the ore sample on the said intersection map to obtain a petrophysical quantity plate for distinguishing the ore body and the surrounding rock.

Specifically, the wave velocity measurement in step 103 adopts an elastic parameter measurement method. The ultrasonic pulse transmission method is used to measure the longitudinal wave velocity and the shear wave velocity of the sample (mineral sample and surrounding rock sample), and the volume product method is used to measure the sample (mineral sample and surrounding rock sample). surrounding rock sample) density.

The calculation methods of the longitudinal wave impedance and the shear wave impedance in step 103 are: longitudinal wave impedance = longitudinal wave velocity × density; shear wave impedance = shear wave velocity × density.

The making of the petrophysical quantity plate in the step 103 is realized by intersecting the longitudinal wave impedance and the shear wave impedance. The abscissa axis of the petrophysical quantity plate is the shear wave impedance, and the ordinate axis is the longitudinal wave impedance. The impedance of the ore-bearing sample and the surrounding rock sample is calculated The results are projected on the intersection map, and the areas of ore-bearing samples and surrounding rock samples are circled respectively, that is, the production of petrophysical quantity plates is completed.

In step 104, the three-dimensional seismic data of the area to be surveyed around the industrial hole is collected by a seismograph.

The 3D seismic data acquisition in step 104 is to collect seismic field measured data through seismographs, the bins of the observation system are ≤10m×10m, the number of coverage times is not less than 48, and the maximum offset is determined according to the exploration depth. The calculation method is: maximum offset Shift distance ≥ exploration depth × 1.5.

For the three-dimensional seismic data acquisition in step 104, the main frequency of the geophone is not higher than 10 Hz, and the vibrator or the explosive source can be selected as the seismic source.

Step 105, using seismic data processing software to perform stacking and migration processing on the 3D seismic data to obtain 3D seismic processing data; the seismic processing data includes pure wave data and result data.

Step 106, using the well logging data and the mud logging data to interpret the 3D seismic processing data to obtain seismic interpretation data of the 3D seismic data, the seismic interpretation data including the horizon and fracture position of the ore-bearing target layer .

In step 106, using the well logging data and the mud logging data to interpret the three-dimensional seismic processing data to obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data includes layers and fractures of ore-bearing target layers location, specifically including: using 3D seismic data interpretation software to perform synthetic seismic recording operations on the 3D seismic processing data and the well logging data, and to complete the calibration of the horizon of the ore-bearing target layer of the seismic data; According to the three-dimensional seismic processing data of the stratum, the fracture position of the ore-bearing target layer is determined according to the feature of the event axis of the three-dimensional seismic processing data.

Seismic data interpretation in step 106, the specific method is to use seismic data interpretation software to perform synthetic seismic records on well logging data to complete the calibration of the target horizon of seismic data, and then perform horizon tracking on seismic data, and at the same time according to The epicenter feature of the seismic data identifies the fault location, and finally interprets the horizon and fault of the target layer on the seismic data volume.

Step 107: Using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume of the 3D seismic data, and obtain the P-wave impedance distribution of the ore-bearing target layer and shear wave impedance distribution.

In step 107, using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume of the 3D seismic data, and obtain the compressional wave impedance of the ore-bearing target layer distribution and shear wave impedance distribution, specifically including: using 3D seismic data processing software to process 3D seismic data to obtain seismic prestack gather data volume, and obtaining seismic prestack gather data volume; using the logging data, the logging Well data and the seismic interpretation data, performing seismic pre-stack inversion on the seismic pre-stack gather data volume, to obtain the distribution of compressional wave impedance and shear wave impedance of the ore-bearing target layer.

The seismic pre-stack inversion in step 107 refers to using seismic pre-stack inversion software to carry out seismic pre-stack inversion on the seismic pre-stack gather data volume to obtain the distribution of longitudinal wave impedance parameters and shear wave impedance parameters of the target layer distributed. The specific method is, first, carry out synthetic seismic recording operation on well logging data, seismic pre-stack gather data volume and seismic interpretation data, complete the well-seismic calibration of seismic data, obtain the seismic data after well-seismic calibration, and then obtain the seismic data from well-seismic Seismic wavelets are extracted from the calibrated seismic data, combined with seismic data and well logging data, the inversion initial model is established, seismic inversion is carried out, and the compressional wave impedance parameters and shear wave impedance parameters of the target layer are obtained through inversion calculation. , set the initial P-wave impedance parameters and S-wave impedance parameters of the target layer, use the initial P-wave impedance parameters and S-wave impedance parameters of the target layer to establish an initial inversion model, and perform convolution operations on the seismic wavelet and the initial inversion model to obtain positive The simulated seismic data is compared with the measured seismic data. When the difference between the two is within the allowable threshold range, the P-wave impedance parameters and S-wave impedance parameters are output. When the difference is not within the allowable threshold range, adjust P-wave impedance parameters and shear-wave impedance parameters, the inversion model is re-established until the difference between the forward seismic data and the measured seismic data is within the allowable threshold.

Step 108, according to the distribution of compressional wave impedance and shear wave impedance of the ore-bearing target layer, use the petrophysical quantity board to delineate the extension range of the sandstone uranium ore body.

In step 108, based on the results of compressional wave impedance and shear wave impedance of the target layer, the distribution range of the sandstone uranium ore body can be delineated by using the petrophysical quantity board.

The present invention also provides a delineation system for the extension range of sandstone uranium ore body, the delineation system includes:

The well logging and mud logging module is used to perform well logging and mud logging operations on industrial holes, and obtain well logging data and mud logging data; the well logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data Includes geological stratification and lithology histograms for industrial bores.

The well logging and mud logging module specifically includes: a well logging sub-module for measuring sound waves and densities at different depths of industrial holes with a well logging instrument; a drilling and coring sub-module for drilling and coring industrial holes, Obtain the rock core; determine the geological stratification sub-module, which is used to determine the geological stratification of the industrial hole according to the lithology of each section of the core; establish a sub-module according to the lithology of each geological layer of the core Create a lithology histogram.

The sample acquisition module is used to take samples in the ore-bearing target layer of the industrial hole, and obtain multiple ore-bearing samples and multiple surrounding rock samples;

The petrophysical quantity plate building module is used to make a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks according to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample.

The petrophysical quantity plate establishment module specifically includes: a wave velocity measurement submodule, which is used to measure the longitudinal wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by ultrasonic pulse transmission method; the density measurement submodule is used for The density of each ore-bearing sample and each surrounding rock sample is measured by the volumetric method; the impedance calculation sub-module is used to compare the compressional wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample with each Multiply the density of the ore sample and each surrounding rock sample to obtain the longitudinal wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample; the coordinate system establishment sub-module is used to use the shear wave impedance as the abscissa of the petrophysical quantity plate Axis, taking the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate to establish a coordinate system; the intersection point projection sub-module is used to project the intersection point of the longitudinal wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample to the On the coordinate system, a cross map is obtained; the ore-bearing rock impedance distribution range delineation sub-module is used to circle the area of ore-bearing samples on the cross map, and obtain a petrophysical quantity board for distinguishing ore bodies and surrounding rocks.

The three-dimensional seismic data acquisition module is used to collect the three-dimensional seismic data of the area to be investigated around the industrial hole through the seismograph;

The three-dimensional seismic data processing module is used to use seismic data processing software to perform superposition processing and migration processing on the three-dimensional seismic data to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and result data;

A three-dimensional seismic data interpretation module, configured to interpret the three-dimensional seismic processing data by using the well logging data and the mud logging data, and obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data including the ore-bearing target layer layers and fracture locations.

The three-dimensional seismic data interpretation module specifically includes: a layer determination sub-module, which is used to use three-dimensional seismic data interpretation software to perform synthetic seismic recording operations on the three-dimensional seismic processing data and the logging data, and complete the seismic data interpretation. Calibration of the horizon of the ore-bearing target layer; the fracture position determination sub-module is used to track the 3D seismic processing data of the layer of the ore-bearing target layer, and determine the location of the ore-bearing target layer according to the synaxial characteristics of the 3D seismic processing data break location.

The three-dimensional seismic data inversion module is used to perform seismic pre-stack inversion on the seismic pre-stack gather data volume of the three-dimensional seismic data by using the logging data, the mud logging data and the seismic interpretation data to obtain ore-bearing P-wave impedance distribution and shear-wave impedance distribution of the target layer.

The three-dimensional seismic data inversion module specifically includes: a data processing sub-module, which is used to process three-dimensional seismic data using three-dimensional seismic data processing software to obtain a seismic pre-stack gather data volume, and obtain a seismic pre-stack gather data volume; The impedance inversion sub-module is used to perform seismic prestack inversion on the seismic prestack gather data volume by using the well logging data, the mud logging data and the seismic interpretation data, and obtain the ore-bearing target layer P-wave impedance distribution and shear-wave impedance distribution.

The delineation module is used to delineate the extension range of the sandstone uranium ore body by using the petrophysical quantity board according to the distribution of the longitudinal wave impedance and the shear wave impedance of the ore-bearing target layer.

According to the specific embodiments provided by the invention, the invention discloses the following technical effects:

The invention discloses a method and system for delineating the extension range of a sandstone uranium ore body. The method comprises the following steps: performing well logging and mud logging in an industrial hole to obtain well logging data, geological stratification and lithology histogram; taking cores from the ore-bearing target layer of the industrial hole to measure the compressional wave velocity, S-wave velocity and density, calculate P-wave impedance and S-wave impedance, and make petrophysical quantity boards for distinguishing ore bodies and surrounding rocks; collect 3D seismic data around industrial holes, complete data processing and interpretation, and use seismic pre-stack reflection to complete data processing and interpretation. Obtain the P-wave impedance distribution and S-wave impedance distribution of the target layer, delineate the distribution range of ore-bearing rocks on the P-wave impedance distribution map and S-wave impedance distribution map according to the rock physical quantity plate, and obtain the extension range of the sandstone uranium ore body. This method provides a technical method for the delineation of the extension range of sandstone uranium deposits, which is environmentally friendly and speeds up the exploration cycle of sandstone uranium deposits.

The technical means such as well logging, mud logging, rock sample measurement, three-dimensional seismic data acquisition, processing, interpretation, and inversion mentioned in the present invention are all general technical means in the field, and can be obtained without departing from the purpose of the present invention. Make various changes. As long as the three-dimensional seismic data, well logging data, mud logging data and rock sample elastic data are comprehensively used, according to the scheme of the patent of the present invention, the extension range of the sandstone uranium ore body can be realized by adopting the common existing technology in this field. Circled.

Each embodiment in this specification is described in a progressive manner, each embodiment focuses on the difference from other embodiments, and the same and similar parts of each embodiment can be referred to each other.

In this paper, specific examples are used to illustrate the principle and implementation of the invention. The description of the above embodiments is only used to help understand the method of the present invention and its core idea. The described embodiments are only part of the embodiments of the present invention. , not all of the embodiments, based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work, all belong to the protection scope of the present invention.

Claims (8)

1. a delineation method of sandstone uranium ore body extension range, it is characterized in that, described delineation method comprises the steps: Perform well logging and mud logging operations on industrial holes to obtain well logging data and mud logging data; the logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data includes geological layers and rock formations of industrial holes sex histogram; Sampling is carried out in the ore-bearing target layer of the industrial hole to obtain multiple ore-bearing samples and multiple surrounding rock samples; According to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample, a petrophysical quantity board for distinguishing ore bodies and surrounding rocks is made, including: The longitudinal wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample were measured by ultrasonic pulse transmission method; The density of each ore-bearing sample and each surrounding rock sample is measured separately by the volumetric method; Multiply the compressional wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by the density of each ore-bearing sample and each surrounding rock sample to obtain the compressional wave of each ore-bearing sample and each surrounding rock sample impedance and shear wave impedance; The coordinate system is established by taking the shear wave impedance as the abscissa axis of the petrophysical quantity plate, and taking the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate; projecting the intersection points of the longitudinal wave impedance and the shear wave impedance of each ore-bearing sample and each surrounding rock sample onto the coordinate system to obtain an intersection map; Circle the area of the ore-bearing sample and the area of the surrounding rock sample on the intersection map, and obtain the petrophysical quantity plate for distinguishing the ore body and the surrounding rock; Acquiring three-dimensional seismic data of the area to be surveyed around the industrial hole through a seismograph; Using seismic data processing software, performing superposition processing and migration processing on the three-dimensional seismic data to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and result 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 includes layers and fracture positions of ore-bearing target layers; Using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume of the 3D seismic data, and obtain the compressional wave impedance distribution and shear wave impedance of the ore-bearing target layer distributed; According to the distribution of longitudinal wave impedance and shear wave impedance of the ore-bearing target layer, the extension range of the sandstone uranium ore body is delineated by using the petrophysical quantity plate. 2. The delineation method of the sandstone uranium ore body extension range according to claim 1, characterized in that, said logging and mud logging operations are carried out to industrial holes to obtain logging data and mud logging data, specifically comprising: Use logging tools to measure sound waves and densities at different depths in industrial boreholes; Drilling and coring industrial holes to obtain cores; Determine the geological stratification of industrial holes according to the lithology of each section of the core; A lithology histogram is built according to the lithology of each geological layer of the core. 3. The method for delineating the extension range of a sandstone uranium ore body according to claim 1, characterized in that, using the logging data and the logging data to interpret the three-dimensional seismic processing data to obtain a three-dimensional Seismic interpretation data of seismic data, the seismic interpretation data includes horizons and fracture positions of ore-bearing target layers, specifically including: Using three-dimensional seismic data interpretation software, performing a synthetic seismic recording operation on the three-dimensional seismic processing data and the logging data, and completing the calibration of the horizon of the ore-bearing target layer of the seismic data; The three-dimensional seismic processing data of the layer of the ore-bearing target layer is tracked, and the fracture position of the ore-bearing target layer is determined according to the synaxial characteristics of the three-dimensional seismic processing data. 4. the delineation method of sandstone uranium ore body extension range according to claim 1, is characterized in that, described utilizing described well logging data, described mud logging data and described seismic interpretation data, to the three-dimensional seismic data Seismic pre-stack inversion is performed on the seismic pre-stack gather data volume to obtain the distribution of P-wave impedance and S-wave impedance of the ore-bearing target layer, including: Use the 3D seismic data processing software to process the 3D seismic data to obtain the seismic pre-stack gather data volume; Using the well logging data, the mud logging data and the seismic interpretation data, perform seismic prestack inversion on the seismic prestack gather data volume to obtain the distribution of compressional wave impedance and shear wave impedance of the ore-bearing target layer. 5. A system for delineating the extension range of a sandstone uranium ore body, characterized in that the system for delineating includes: The well logging and mud logging module is used to perform well logging and mud logging operations on industrial holes, and obtain well logging data and mud logging data; the well logging data includes sound waves and densities at different depths of industrial holes, and the mud logging data Geological stratification and lithology histograms including industrial boreholes; The sample acquisition module is used to take samples in the ore-bearing target layer of the industrial hole, and obtain multiple ore-bearing samples and multiple surrounding rock samples; The petrophysical quantity plate building module is used to make a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks according to the compressional wave impedance and shear wave impedance of each ore-bearing sample and each surrounding rock sample; The petrophysical quantity plate building module specifically includes: The wave velocity measurement sub-module is used to measure the longitudinal wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by ultrasonic pulse transmission method; The density measurement sub-module is used to measure the density of each ore-bearing sample and each surrounding rock sample respectively by the volumetric method; The impedance calculation sub-module is used to multiply the compressional wave velocity and shear wave velocity of each ore-bearing sample and each surrounding rock sample by the density of each ore-bearing sample and each surrounding rock sample to obtain each ore-bearing sample and each The compressional wave impedance and shear wave impedance of each surrounding rock sample; The coordinate system establishment sub-module is used to establish a coordinate system with the shear wave impedance as the abscissa axis of the petrophysical quantity plate and the longitudinal wave impedance as the ordinate axis of the petrophysical quantity plate; The intersection point projection sub-module is used to project the intersection point of the longitudinal wave impedance and the shear wave impedance of each ore-bearing sample and each surrounding rock sample onto the coordinate system to obtain an intersection map; The ore-bearing rock impedance distribution range delimitation sub-module is used to circle the ore-bearing sample area on the intersection map, and obtain a petrophysical quantity plate for distinguishing ore bodies and surrounding rocks; The three-dimensional seismic data acquisition module is used to collect the three-dimensional seismic data of the area to be investigated around the industrial hole through the seismograph; The three-dimensional seismic data processing module is used to use seismic data processing software to perform superposition processing and migration processing on the three-dimensional seismic data to obtain three-dimensional seismic processing data; the seismic processing data includes pure wave data and result data; A three-dimensional seismic data interpretation module, configured to interpret the three-dimensional seismic processing data by using the well logging data and the mud logging data, and obtain seismic interpretation data of the three-dimensional seismic data, the seismic interpretation data including the ore-bearing target layer horizons and fracture locations; The three-dimensional seismic data inversion module is used to perform seismic pre-stack inversion on the seismic pre-stack gather data body of the three-dimensional seismic data by using the well logging data, the mud logging data and the seismic interpretation data to obtain ore-bearing P-wave impedance distribution and shear-wave impedance distribution of the target layer; The delineation module is used to delineate the extension range of the sandstone uranium ore body by using the petrophysical quantity board according to the distribution of the longitudinal wave impedance and the shear wave impedance of the ore-bearing target layer. 6. The delineation system of the sandstone uranium ore body extension range according to claim 5, characterized in that, the well logging and mud logging modules specifically include: The logging sub-module is used to measure the acoustic wave and density of industrial holes at different depths with the logging instrument; The drilling and coring sub-module is used for drilling and coring industrial holes to obtain cores; The sub-module for determining geological stratification is used to determine the geological stratification of industrial holes according to the lithology of each section of the core; The lithology histogram establishment sub-module is used for establishing the lithology histogram according to the lithology of each geological layer of the core. 7. The delineation system of the sandstone uranium ore body extension range according to claim 5, wherein the three-dimensional seismic data interpretation module specifically includes: The horizon determination sub-module is used to use the three-dimensional seismic data interpretation software to perform synthetic seismic recording operation on the three-dimensional seismic processing data and the well logging data, and complete the calibration of the horizon of the ore-bearing target layer of the seismic data; The sub-module for determining the fracture position is used to track the 3D seismic processing data of the horizon of the ore-bearing target layer, and determine the fracture position of the ore-bearing target layer according to the concentric axis characteristics of the 3D seismic processing data. 8. The delineation system of the sandstone uranium ore body extension range according to claim 5, wherein the three-dimensional seismic data inversion module specifically includes: The data processing sub-module is used to process the three-dimensional seismic data by using the three-dimensional seismic data processing software to obtain the seismic pre-stack gather data body; The impedance inversion sub-module is used to perform seismic prestack inversion on the seismic prestack gather data volume by using the well logging data, the mud logging data and the seismic interpretation data, and obtain the ore-bearing target layer P-wave impedance distribution and shear-wave impedance distribution.

Images