CN115291288B - While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking - Google Patents

While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking Download PDF

Info

Publication number
CN115291288B
CN115291288B CN202211194628.0A CN202211194628A CN115291288B CN 115291288 B CN115291288 B CN 115291288B CN 202211194628 A CN202211194628 A CN 202211194628A CN 115291288 B CN115291288 B CN 115291288B
Authority
CN
China
Prior art keywords
drilling
porosity
pulse
neutron
neutrons
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
CN202211194628.0A
Other languages
Chinese (zh)
Other versions
CN115291288A (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.)
China University of Petroleum East China
Sinopec Jingwei Co Ltd
Original Assignee
China University of Petroleum East China
Sinopec Jingwei Co Ltd
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 China University of Petroleum East China, Sinopec Jingwei Co Ltd filed Critical China University of Petroleum East China
Priority to CN202211194628.0A priority Critical patent/CN115291288B/en
Publication of CN115291288A publication Critical patent/CN115291288A/en
Application granted granted Critical
Publication of CN115291288B publication Critical patent/CN115291288B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V5/00Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity
    • G01V5/04Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging
    • G01V5/08Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays
    • G01V5/10Prospecting or detecting by the use of nuclear radiation, e.g. of natural or induced radioactivity specially adapted for well-logging using primary nuclear radiation sources or X-rays using neutron sources
    • 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

Abstract

The invention discloses an intelligent processing method for porosity of pulse neutrons while drilling based on iron neutron marking, and relates to the technical field of geophysical logging in a mine field. The method includes the steps of measuring by a pulse-while-drilling neutron porosity logging instrument, measuring a standard graduated well by the pulse-while-drilling neutron porosity logging instrument, simulating by a Monte Carlo numerical simulation method, establishing a graduation relation between a simulation value and an actually-measured value of the pulse-while-drilling neutron porosity logging instrument, simulating by the Monte Carlo numerical simulation method to obtain pulse-while-drilling neutron porosity samples in different logging-while-drilling environments, constructing a sample database and a multilayer sensing neural network, training the multilayer sensing neural network by the sample database to calculate the pulse-while-drilling neutron porosity and verify the accuracy of a calculation result, inputting measured data of the pulse-while-drilling neutron porosity logging instrument after graduation processing into the multilayer sensing neural network, calculating the pulse-while-drilling neutron porosity, and improving the calculation accuracy of the pulse-while-drilling neutron porosity.

Description

While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking
Technical Field
The invention relates to the technical field of geophysical logging in a mine field, in particular to an intelligent processing method for porosity of pulse neutrons while drilling based on iron neutron marking.
Background
The accurate measurement of the formation porosity is taken as the key of high-precision formation evaluation, because the neutron porosity logging instrument adopting the traditional chemical source is limited by a series of problems of health, safety, environment and the like in development and application, the neutron porosity logging based on the pulse neutron source becomes the development direction of the porosity logging technology, the emission of neutrons can be controlled by adopting the D-T pulse neutron source, and the radiation risk of the traditional chemical source in the aspects of transportation, construction, storage and the like is avoided.
Compared with the traditional chemical source, the high-energy fast neutrons energy released by the pulse neutron source is 14MeV, and the neutron energy-space-time field distribution is different from that of the traditional chemical source in the process of acting with the formation medium, so that the porosity measurement sensitivity and the influence factors of the high-energy fast neutrons energy released by the pulse neutron source are greatly different from that of the conventional compensated neutrons, and the research in the present stage finds that the high-energy fast neutrons emitted by the pulse neutron source have the phenomenon that the neutron moderation process tends to be saturated in the high-porosity formation, so that the sensitivity of detecting the formation porosity by utilizing thermal neutrons is obviously lower than that of the chemical source, and how to eliminate the phenomenon of neutron moderation saturation in the high-porosity formation is realized, thereby improving the pulse neutron porosity measurement sensitivity to be the key for accurately obtaining the formation porosity.
The logging-while-drilling technology is used as an important means applied to exploration and development of unconventional reservoirs, deep oil gas and deep water oil gas, the drill collar is always in a vibration state in the measurement-while-drilling process, so that a gap exists between the pulsed neutron porosity instrument and a well wall, and meanwhile mud, formation water mineralization and the like adopted in the measurement-while-drilling process can also influence the measurement result of the pulsed neutron porosity while drilling. The traditional neutron porosity processing method corrects the influence of the environment on the neutron porosity by establishing a plate, but in the actual measurement process, the instrument gap and the formation condition are usually difficult to directly obtain. Therefore, it is desirable to provide a method for correcting the influence of environmental factors by using the measurement information of the apparatus itself, so as to improve the measurement accuracy of the porosity of the pulsed neutrons while drilling.
Disclosure of Invention
Aiming at the problem that the measurement precision of a pulse neutron porosity while drilling result is low due to the influence of a measurement-while-drilling environment, the invention discloses an intelligent processing method for the pulse neutron porosity while drilling based on iron neutron marking.
In order to realize the purpose, the invention adopts the following technical scheme:
an intelligent processing method for porosity of pulse neutrons while drilling based on iron neutron marking adopts a pulse neutron porosity logging while drilling instrument to measure, and specifically comprises the following steps:
step 1), placing a pulse neutron porosity logging-while-drilling instrument in a standard calibration well for measurement, obtaining a thermal neutron count and a neutron time spectrum of a near neutron detector, a thermal neutron count and a neutron time spectrum of a far neutron detector and a gamma energy spectrum of a gamma detector, and calculating to obtain measured values of the pulse neutron porosity logging-while-drilling instrument, wherein the measured values comprise a measured value of a macroscopic capture cross section of a borehole medium, a measured value of a macroscopic capture cross section of a formation medium, a measured value of a thermal neutron count ratio, a measured value of a non-bomb gamma count, a measured value of a capture gamma count, a measured value of an iron element characteristic gamma count and a measured value of a chlorine element characteristic gamma count;
step 2), according to the instrument structure of the pulse neutron porosity logging-while-drilling instrument, in combination with the structural parameters of a standard graduated well, establishing a pulse neutron porosity value calculation model while drilling by using a Monte Carlo numerical simulation method, and simulating to obtain a simulated value of the pulse neutron porosity logging-while-drilling instrument, wherein the simulated value comprises a simulated value of a macroscopic capture cross section of a borehole medium, a simulated value of a macroscopic capture cross section of a formation medium, a simulated value of a thermal neutron count ratio, a simulated value of a non-elastic gamma count, a simulated value of a capture gamma count, a simulated value of an iron element characteristic gamma count and a simulated value of a chlorine element characteristic gamma count;
step 3), establishing a scale relation between the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument based on the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument;
step 4), changing the size of a drill collar of the pulse neutron porosity logging while drilling instrument according to the instrument structure of the pulse neutron porosity logging while drilling instrument, sequentially changing the lithology, the formation water mineralization, the drill collar gap, the mud type and the formation porosity of a formation medium under different drill collar sizes, reestablishing a pulse neutron porosity value calculation model while drilling, simulating and obtaining the thermal neutron count and the neutron time spectrum of a near neutron detector, the thermal neutron count and the neutron time spectrum of a far neutron detector and the gamma energy spectrum of a gamma detector under different logging while drilling environments, and forming a plurality of groups of pulse neutron porosity while drilling samples according to the drill collar sizes, the macro capture section of a borehole medium, the macro capture section of the formation medium, the counting ratio of iron-marked thermal neutrons, the non-elastic gamma count, the capture gamma count and the gamma count of chlorine element characteristics under different logging while drilling environments;
step 5), constructing a sample database based on the multiple groups of while-drilling pulse neutron porosity samples obtained in the step 4), and dividing the sample database into a training set and a verification set;
step 6), constructing a multilayer perception neural network, and calculating the porosity of the while-drilling pulse neutrons by using the training set to train the multilayer perception neural network;
step 7), verifying the accuracy of the porosity of the pulse-while-drilling neutrons by using the multi-layer perceptive neural network after training by using a verification set;
and 8) based on the scale relation between the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument established in the step 3), performing scale processing on a borehole medium macroscopic capture section, a formation medium macroscopic capture section, a thermal neutron count ratio, a non-ballistic gamma count, a capture gamma count, an iron element characteristic gamma count and a chlorine element characteristic gamma count which are actually measured by the pulse-while-drilling neutron porosity logging instrument, marking the thermal neutron count ratio by using the iron element characteristic gamma count, inputting the borehole medium macroscopic capture section, the formation medium macroscopic capture section, the non-ballistic gamma count, the capture gamma count, the chlorine element characteristic gamma count, the iron-marked thermal neutron count ratio and the drill collar size of the pulse-while-drilling neutron porosity logging instrument after scale processing into a multilayer sensing neural network, and calculating by using the multilayer sensing neural network to obtain the pulse-while-drilling neutron porosity.
Preferably, the pulse-while-drilling neutron porosity logging instrument is arranged in a drill collar, a slurry flow guide channel is arranged in the drill collar, the outer wall of a shell of the pulse-while-drilling neutron porosity logging instrument is attached to the inner wall of the drill collar, and a far neutron detector, a gamma detector, a near neutron detector, a tungsten-nickel-iron shield and a pulse neutron source are arranged in the shell at intervals from top to bottom.
Preferably, the near neutron detector and the far neutron detector are both in a cylinder structure, the diameter of the near neutron detector is set to be 4cm, the length of the near neutron detector is set to be 5cm, the diameter of the far neutron detector is set to be 6cm, the length of the far neutron detector is set to be 10cm, the diameter of the gamma detector is set to be 4cm, the length of the gamma detector is set to be 4cm, the source distance of the near neutron detector is set to be 35cm, the source distance of the far neutron detector is set to be 70cm, and the source distance of the gamma detector is set to be 50cm.
Preferably, the step 3) includes a calibration relation between the borehole medium macroscopic capture section simulation value and the measured value, a calibration relation between the formation medium macroscopic capture section simulation value and the measured value, a calibration relation between the thermal neutron count ratio simulation value and the measured value, a calibration relation between the non-bomb-gamma-count simulation value and the measured value, a calibration relation between the capture gamma-count simulation value and the measured value, a calibration relation between the iron element characteristic gamma-count simulation value and the measured value, and a calibration relation between the chlorine element characteristic gamma-count simulation value and the measured value.
Preferably, in the step 4), the drill collar sizes of the logging-while-drilling pulsed neutron porosity logging instrument are sequentially set to be 12.065cm, 17.145cm and 20.32cm, the lithology of the formation medium is sequentially set to be sandstone, limestone, dolomite, sand-ash mixed rock and sand-dolomite mixed rock, wherein the sand-ash mixed rock is formed by mixing sandstone and limestone, the sand-dolomite mixed rock is formed by mixing sandstone and dolomite, the water mineralization of the formation is sequentially set to be 0kppm, 5kppm, 10kppm, 20kppm, 50kppm, 100kppm and 200kppm, the drill collar gaps are sequentially set to be 0cm, 0.5cm, 1cm, 1.5cm and 2cm, the mud types include oil-based mud and water-based mud, and the formation porosity is sequentially set to be 0, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45% and 50%.
Preferably, in the step 4), the thermal neutron count ratio is marked by using the characteristic gamma count of the iron element, so as to obtain an iron-marked thermal neutron count ratio, as shown in formula (1):
Figure 100002_DEST_PATH_IMAGE002
(1)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE004
the counting ratio of the thermal neutrons is marked for the iron,
Figure 100002_DEST_PATH_IMAGE006
is the counting ratio of the thermal neutrons,
Figure 100002_DEST_PATH_IMAGE008
and (4) performing gamma counting for the characteristic iron element after normalization treatment.
Preferably, the step 6) specifically comprises the following steps:
6.1, constructing a multilayer perception neural network, wherein the multilayer perception neural network is provided with an input layer, a first hidden layer, a second hidden layer and an output layer, the input layer, the first hidden layer, the second hidden layer and the output layer are sequentially connected, 7 nodes are arranged in the input layer and are used for obtaining input parameters of the multilayer perception neural network, the input parameters comprise the size of a drill collar in a pulse neutron porosity while drilling sample, a macroscopic capture section of a borehole medium, a macroscopic capture section of a formation medium, a counting ratio of iron-marked thermal neutrons, a non-elastic gamma count, a capture gamma count and a chlorine element characteristic gamma count, 10 nodes are arranged in the first hidden layer, 15 nodes are arranged in the second hidden layer, the first hidden layer and the second hidden layer are connected through a ReLU activation function, the second hidden layer and the output layer are connected through a tanh activation function, and 1 node is arranged in the output layer and is used for outputting the calculated pulse neutron porosity while drilling;
step 6.2, setting a porosity precision value of the pulse neutrons while drilling;
step 6.3, randomly selecting a pulse neutron porosity while drilling sample in the training set, inputting the pulse neutron porosity while drilling sample into the multilayer perception neural network, enabling nodes in the input layer of the multilayer perception neural network to obtain input parameters in the pulse neutron porosity while drilling sample, activating the input parameters through a ReLU function, and transmitting the activated input parameters to the first hidden layer, wherein the output result of the first hidden layer is shown as a formula (2):
Figure 100002_DEST_PATH_IMAGE010
(2)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE012
as a result of the output of the first hidden layer,
Figure 100002_DEST_PATH_IMAGE014
is an input parameter of the multi-layer perceptive neural network,
Figure 100002_DEST_PATH_IMAGE016
is the size of a drill collar,
Figure 100002_DEST_PATH_IMAGE018
the ratio of the counts of the thermal neutrons to the iron marker,
Figure 100002_DEST_PATH_IMAGE020
the counting is carried out for the non-bullet gamma counting,
Figure 100002_DEST_PATH_IMAGE022
in order to capture the gamma counts,
Figure 100002_DEST_PATH_IMAGE024
is the gamma counting of the characteristic of the chlorine element,
Figure 100002_DEST_PATH_IMAGE026
for the macroscopic capture cross-section of the wellbore medium,
Figure 100002_DEST_PATH_IMAGE028
is a macroscopic trapping cross section of the formation medium,
Figure 100002_DEST_PATH_IMAGE030
is the weight of the first hidden layer,
Figure 100002_DEST_PATH_IMAGE032
is a first hidingBiasing of the layers;
the output result of the first hidden layer is activated by the ReLU function and then transmitted to the second hidden layer, and the output result of the second hidden layer is shown in formula (3):
Figure 100002_DEST_PATH_IMAGE034
(3)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE036
is the output result of the second hidden layer,
Figure 100002_DEST_PATH_IMAGE038
is the weight of the second hidden layer,
Figure 100002_DEST_PATH_IMAGE040
a bias for the second hidden layer;
the output result of the second hidden layer is activated by the tanh function and then transmitted to the output layer, and the porosity of the while-drilling pulse neutrons output by the output layer is as follows:
Figure 100002_DEST_PATH_IMAGE042
(4)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE044
the porosity of the pulse neutron output by the output layer while drilling,
Figure 100002_DEST_PATH_IMAGE046
as a weight of the output layer,
Figure 100002_DEST_PATH_IMAGE048
is the bias of the output layer;
step 6.4, calculating the formation porosity of the model by combining the porosity value of the while-drilling pulse neutrons corresponding to the while-drilling pulse neutron porosity sample according to the porosity of the while-drilling pulse neutrons output by the output layer, calculating the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network, and comparing the porosity error value with the preset precision value of the porosity of the while-drilling pulse neutrons; if the while-drilling pulse neutron porosity error value of the multilayer perception neural network is larger than the set while-drilling pulse neutron porosity precision value, updating the weight and bias of each layer in the multilayer perception neural network by using a logcosh loss function and adopting a back propagation and random gradient descent method, returning to the step 6.3, and continuing to train the multilayer perception neural network; and if the while-drilling pulse neutron porosity error value of the multilayer perception neural network is not less than the set while-drilling pulse neutron porosity precision value, stopping training the multilayer perception neural network, and entering the step 7).
Preferably, in step 6.4, the update formula of the weights and the biases of each layer in the multilayer perceptive neural network is as follows:
Figure 100002_DEST_PATH_IMAGE050
(5)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE052
to minimize the value of the variable at which the logcosh loss function takes a minimum,
Figure 100002_DEST_PATH_IMAGE054
as logcosh loss function.
Preferably, in the step 7), all the pulse-while-drilling neutron porosities in the verification set are input into the trained multilayer perceptive neural network, a calculation value of the pulse-while-drilling neutron porosity is calculated by using the multilayer perceptive neural network for each pulse-while-drilling neutron porosity sample, the formation porosity of the model is calculated by combining the pulse-while-drilling neutron porosity values corresponding to the pulse-while-drilling neutron porosity samples, an absolute error value of the pulse-while-drilling neutron porosity is calculated by using the multilayer perceptive neural network by using an absolute error function, and if the absolute error value is less than 1%, the calculation result of the pulse-while-drilling neutron porosity of the multilayer perceptive neural network is proved to be accurate;
the calculation formula of the absolute error value of the porosity of the while-drilling pulse neutrons is as follows:
Figure 100002_DEST_PATH_IMAGE056
(6)
in the formula (I), the compound is shown in the specification,
Figure 100002_DEST_PATH_IMAGE058
is the absolute error value of the porosity of the while-drilling pulse neutrons,
Figure 100002_DEST_PATH_IMAGE060
to verify the number of concentrated while-drilling pulsed neutron porosity samples,
Figure 100002_DEST_PATH_IMAGE062
in order to verify the serial number of the porosity sample of the concentrated while-drilling pulse neutron,
Figure 100002_DEST_PATH_IMAGE064
for multi-layer perceptive neural network according to
Figure 100002_DEST_PATH_IMAGE062A
A pulse-while-drilling neutron porosity value calculated for each pulse-while-drilling neutron porosity sample,
Figure 100002_DEST_PATH_IMAGE066
is as follows
Figure 100002_DEST_PATH_IMAGE062AA
And calculating the formation porosity of the model according to the pulse-while-drilling neutron porosity value corresponding to the pulse-while-drilling neutron porosity sample.
The invention has the beneficial effects that:
according to the intelligent processing method for the porosity of the pulse neutrons while drilling based on the iron neutron marker, provided by the invention, the iron element characteristic gamma counting is used as a marking factor to process the counting ratio of the thermal neutrons, so that the influence of the energy of a pulse neutron source on the distribution of a neutron field is eliminated, and the measurement sensitivity of a pulse neutron porosity logging-while-drilling instrument is improved.
Meanwhile, the method trains a multilayer perception neural network by using a pulse-while-drilling neutron porosity sample obtained by Monte Carlo simulation, inputs the drill collar size of the pulse-while-drilling neutron porosity instrument, the measured macroscopic capture cross section of a borehole medium, the macroscopic capture cross section of a stratum medium, the iron-marked thermal neutron count ratio, the non-elastic gamma count, the capture gamma count and the chlorine element characteristic gamma count into the multilayer perception neural network as input parameters, calculates and obtains the real pulse-while-drilling neutron porosity by using the multilayer perception neural network, automatically eliminates errors caused by measurement-while-drilling environments such as the borehole environment, the drill collar and the like, improves the calculation precision of the pulse-while-drilling neutron porosity, ensures the accuracy of the calculation result of the pulse-while-drilling neutron porosity, and lays a foundation for high-precision stratum evaluation.
Drawings
FIG. 1 is a schematic structural diagram of a pulse-neutron porosity while drilling logging instrument according to the present invention. In the figure, 1 is a stratum, 2 is a drill collar, 3 is a slurry diversion channel, 4 is a far neutron detector, 5 is a gamma detector, 6 is a near neutron detector, 7 is a tungsten-nickel-iron shield, and 8 is a pulse neutron source.
FIG. 2 is a graph of capture gamma energy for different formation porosity limestone formations.
FIG. 3 is a graph of the characteristic gamma counts of elemental iron as a function of formation porosity.
FIG. 4 is a graph of the response of thermal neutron count ratio and iron-labeled thermal neutron count ratio with formation porosity.
FIG. 5 is a plot of formation porosity sensitivity versus thermal neutron count ratio and iron-labeled thermal neutron count ratio.
Fig. 6 is a curve of logcosh loss function value with training times in the process of multi-layer perceptive neural network training.
FIG. 7 is a distribution frequency diagram for verifying absolute error values of concentrated while-drilling pulsed neutron porosity samples.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
After entering the stratum, fast neutrons emitted by the pulse neutron source react with the nuclei of stratum substances through inelastic scattering, elastic scattering, radiation capture and the like, and the energy of the neutrons is reduced to become epithermal neutrons and thermal neutrons. The distribution of thermal neutron flux in the formation depends on the distance traveled by the neutrons
Figure DEST_PATH_IMAGE068
And the deceleration length of the formation medium
Figure DEST_PATH_IMAGE070
The formation density also has negligible effect on the thermal neutron distribution, so the distribution of thermal neutrons in the formation is affected by both the hydrogen index and the formation density.
According to the theory of double-group diffusion, the ratio of the counts of thermal neutrons recorded by two neutron detectors in an infinite homogeneous medium
Figure DEST_PATH_IMAGE072
The approximate description is:
Figure DEST_PATH_IMAGE074
(7)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE076
is the thermal neutron count rate of the near neutron detector,
Figure DEST_PATH_IMAGE078
is the thermal neutron count rate of the far neutron detector,
Figure DEST_PATH_IMAGE080
is the source distance of the near-neutron detector,
Figure DEST_PATH_IMAGE082
is the source distance of the far neutron detector,
Figure DEST_PATH_IMAGE070A
is the deceleration length of the formation medium.
Because the energy of the pulse neutron source is higher than that of an americium-beryllium neutron source, high-energy neutrons emitted by the pulse neutron source mainly have an inelastic scattering effect with a rock framework after entering a stratum and are greatly influenced by the stratum density, and the porosity sensitivity measured by the pulse neutron porosity logging method in a high-porosity stratum is lower than that of an Am-Be neutron source.
The method is characterized in that the characteristic gamma counting of iron element is only related to the flux distribution of thermal neutrons, the characteristic gamma counting of silicon element and the characteristic gamma counting of calcium element are mainly related to the lithology of the stratum, and the characteristic gamma counting of chlorine element and the macroscopic capture cross section are mainly related to the mineralization of stratum water.
The invention provides an intelligent processing method of pulse-while-drilling neutron porosity based on iron neutron marking, which is characterized in that a pulse-while-drilling neutron porosity logging instrument is adopted to measure in a stratum, as shown in figure 1, the pulse-while-drilling neutron porosity logging instrument is arranged in a drill collar 2, a mud diversion channel 3 is arranged in the drill collar 2, the outer wall of a shell of the pulse-while-drilling neutron porosity logging instrument is attached to the inner wall of the drill collar 2, and a far neutron detector 4, a gamma detector 5, a near neutron detector 6, a tungsten-nickel-iron shield 7 and a pulse neutron source 8 are arranged in the shell at intervals from top to bottom.
In this embodiment, the pulse neutron source 8 is set as a D-T pulse neutron generator, the D-T pulse neutron generator emits fast neutrons into the formation 1 in a pulse form, the near neutron detector and the far neutron detector are both in a cylinder structure, the diameter of the near neutron detector is set to be 4cm, the length of the near neutron detector is set to be 5cm, the diameter of the far neutron detector is set to be 6cm, the length of the far neutron detector is set to be 10cm, the diameter of the gamma detector is set to be 4cm, the length of the gamma detector is set to be 4cm, the source distance of the near neutron detector is set to be 35cm, the source distance of the far neutron detector is set to be 70cm, and the source distance of the gamma detector is set to be 50cm.
The invention provides an intelligent processing method for porosity of pulse neutrons while drilling based on iron neutron marking, which specifically comprises the following steps:
step 1, placing a pulse neutron porosity logging-while-drilling instrument in a standard calibration well for measurement, setting the iron element energy window range of a gamma detector to be 6.4-7.8MeV, the chlorine element energy window range to be 1.02-1.23MeV, setting the non-bomb gamma recording time period of the gamma detector to be 0-20 mu s, setting the capture gamma recording time period to be 60-160 mu s, obtaining the thermal neutron count and neutron time spectrum of a near neutron detector, the thermal neutron count and neutron time spectrum of a far neutron detector and the gamma energy spectrum of the gamma detector, and calculating to obtain the measured value of the pulse neutron porosity logging-while-drilling instrument, the method comprises an actual measurement value of a macroscopic capture section of a borehole medium, an actual measurement value of a macroscopic capture section of a formation medium, an actual measurement value of a thermal neutron count ratio, an actual measurement value of a non-bomb gamma counting, an actual measurement value of a capture gamma counting, an actual measurement value of an iron element characteristic gamma counting and an actual measurement value of a chlorine element characteristic gamma counting, wherein the macroscopic capture section of the borehole medium and the macroscopic capture section of the formation medium are determined by utilizing neutron time spectrums of a near neutron detector and a far neutron detector, the thermal neutron count ratio is determined by utilizing the thermal neutron counting of the near neutron detector and the far neutron detector, and the non-bomb gamma counting, the capture gamma counting, the iron element characteristic gamma counting and the chlorine element characteristic gamma counting are determined by utilizing a gamma energy spectrum.
And 2, establishing a pulse neutron porosity while drilling numerical calculation model by using a Monte Carlo numerical simulation method according to the instrument structure of the pulse neutron porosity while drilling instrument and combining the structural parameters of a standard graduated well, and simulating to obtain a simulated value of the pulse neutron porosity while drilling instrument, wherein the simulated value comprises a simulated value of a macroscopic capture cross section of a borehole medium, a simulated value of a macroscopic capture cross section of a formation medium, a simulated value of a thermal neutron count ratio, a simulated value of a non-elastic gamma count, a simulated value of a capture gamma count, a simulated value of an iron element characteristic gamma count and a simulated value of a chlorine element characteristic gamma count. FIG. 2 is a capture gamma energy spectrum of a limestone formation under different formation porosity conditions, wherein the capture gamma energy spectrum of the limestone formation with the formation porosity of 5%, the capture gamma energy spectrum of the limestone formation with the formation porosity of 20% and the capture gamma energy spectrum of the limestone formation with the porosity of 40% are shown, and the capture gamma energy spectra under different formation porosity conditions are compared, so that the capture gamma counts of the capture gamma energy spectra under the energy range of 6.5 to 7.8MeV are reduced along with the increase of the formation porosity, and the difference between the capture gamma counts in other energy sections is smaller. Fig. 3 shows the change rule of the characteristic gamma count of the iron element with the porosity of the formation, which can be obtained from fig. 3, and even in a high-porosity formation, the characteristic gamma count of the iron element still changes with the porosity of the formation, and by analyzing fig. 2 and fig. 3, the characteristic gamma count of the iron element can reflect the distribution of the thermal neutron field in the formation. Therefore, the influence of the neutron source type on the neutron field distribution can be eliminated by marking the thermal neutron counting ratio by utilizing the characteristic gamma counting of the iron element, and the sensitivity of the instrument on the porosity measurement of the pulse neutrons is improved. Marking the thermal neutron count ratio by utilizing iron element characteristic gamma counting to obtain the iron-marked thermal neutron count ratio as shown in a formula (1):
Figure DEST_PATH_IMAGE002A
(1)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE004A
the counting ratio of the thermal neutrons is marked for the iron,
Figure DEST_PATH_IMAGE006A
is the counting ratio of the thermal neutrons,
Figure DEST_PATH_IMAGE008A
and (4) performing gamma counting for the characteristic iron element after normalization treatment.
FIG. 4 is a graph showing the response law of the thermal neutron count ratio and the iron-labeled thermal neutron count ratio with the formation porosity, which can be obtained from FIG. 4, wherein the thermal neutron count ratios all increase with the increase of the formation porosity, the thermal neutron count ratio remains unchanged after the formation porosity exceeds 20%, and the iron-labeled thermal neutron count ratio has an obvious change law in the whole formation porosity distribution range, a change curve of the formation porosity relative sensitivity of the thermal neutron count ratio with the formation porosity and a change curve of the formation porosity relative sensitivity of the iron-labeled thermal neutron count ratio with the formation porosity are obtained through calculation, as shown in FIG. 5, the formation porosity relative sensitivity of the iron-labeled thermal neutron count ratio is higher than that of the thermal neutron count ratio, and when the formation porosity is 30%, the formation porosity relative sensitivity of the iron-labeled thermal neutron count ratio is about 2.7 times that of the thermal neutron count ratio.
And 3, establishing a scale relation between the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument based on the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument, wherein the scale relation between the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument comprises a scale relation between a macroscopic capturing cross section analog value and the measured value of the borehole medium, a scale relation between a macroscopic capturing cross section analog value and the measured value of the formation medium, a scale relation between a thermal neutron count ratio analog value and the measured value, a scale relation between a non-bomb gamma counting analog value and the measured value, a scale relation between a capturing gamma counting analog value and the measured value, a scale relation between an iron element characteristic gamma counting analog value and the measured value and a scale relation between a chlorine element characteristic gamma counting analog value and the measured value.
And 4, in addition to the source energy difference, the distribution difference of the neutron field in the stratum can be caused by environmental factors such as lithology, stratum water mineralization and the like, and the problem that iron marking factors are difficult to obtain exists in the actual measurement process, so that the measurement sensitivity and the measurement accuracy of the pulse neutron porosity logging-while-drilling instrument under different measurement-while-drilling environments are improved.
According to the instrument structure of the pulse-while-drilling neutron porosity logging instrument, the size of a drill collar of the pulse-while-drilling neutron porosity logging instrument is changed, the lithology, the formation water mineralization, the drill collar gap, the mud type and the formation porosity of a formation medium are sequentially changed under the condition of different drill collar sizes, a pulse-while-drilling neutron porosity value calculation model is reestablished, the thermal neutron count and the neutron time spectrum of a near neutron detector, the thermal neutron count and the neutron time spectrum of a far neutron detector and the gamma energy spectrum of a gamma detector under different logging-while-drilling environments are obtained in a simulated mode, and a plurality of groups of pulse-while-drilling neutron porosity samples are formed according to the drill collar sizes, the macro capture cross section of a borehole medium, the macro capture cross section of the formation medium, the counting ratio of iron-marked thermal neutrons, the non-elastic gamma count, the capture gamma count and the gamma counting of chlorine element characteristics under different logging-while-drilling environments.
In the embodiment, the drill collar sizes of the pulse neutron porosity logging-while-drilling instrument are sequentially set to be 12.065cm, 17.145cm and 20.32cm, after the drill collar sizes are fixed, the lithology of a formation medium, the formation water mineralization, the drill collar gap, the mud type and the formation porosity in the model are calculated by sequentially modifying the value of the pulse neutron porosity while-drilling value, wherein the lithology of the formation medium is sequentially set to be sandstone, limestone, dolomite and sand-lime mixed rock (formed by mixing sandstone and limestone, the proportion of the limestone is increased from 0 to 100%, and is increased by 10% each time), the formation water mineralization is sequentially set to be 0kppm, 5kppm, 10kppm, 20kppm, 50kppm, 100kppm and 200kppm, the drill collar gap is sequentially set to be 0cm, 0.5cm, 1cm, 1.5cm and 2cm, the water-based mud type is sequentially set to be 0kppm, 20kppm, 50kppm, 100kppm and 200kppm, the mud porosity is sequentially set to be 0cm, 0.5cm, 25%, 40%, 35%, and 35%. Different measurement-while-drilling environments are simulated by modifying the calculation model of the porosity value of the pulse-while-drilling neutrons, so that the sizes of drill collars, the macroscopic capture cross section of a borehole medium, the macroscopic capture cross section of a formation medium, the counting ratio of iron-marked thermal neutrons, the non-bomb gamma counting, the capture gamma counting and the chlorine element characteristic gamma counting in the different measurement-while-drilling environments are obtained, and 10000 groups of measurement-while-drilling pulse neutron porosity samples are obtained in the embodiment.
And 5, constructing a sample database based on the multiple groups of pulse-while-drilling neutron porosity samples obtained in the step 4, and dividing the sample database into a training set and a verification set, wherein the number of pulse-while-drilling neutron porosity samples in the training set accounts for 85% of the total number of pulse-while-drilling neutron porosity samples in the sample database, and the number of pulse-while-drilling neutron porosity samples in the verification set accounts for 15% of the total number of pulse-while-drilling neutron porosity samples in the sample database.
Step 6, constructing a multilayer perception neural network, and calculating the porosity of the while-drilling pulse neutrons by using the training set to train the multilayer perception neural network, wherein the method specifically comprises the following steps:
6.1, constructing a multilayer perception neural network, wherein the multilayer perception neural network is provided with an input layer, a first hidden layer, a second hidden layer and an output layer, the input layer, the first hidden layer, the second hidden layer and the output layer are sequentially connected, 7 nodes are arranged in the input layer and used for obtaining input parameters of the multilayer perception neural network, the input parameters comprise drill collar size, borehole medium macroscopic capture cross section, formation medium macroscopic capture cross section, iron marker thermal neutron count ratio, non-bomb gamma counting, capture gamma counting and chlorine element characteristic gamma counting in a while-drilling pulse neutron porosity sample, 10 nodes are arranged in the first hidden layer, 15 nodes are arranged in the second hidden layer, the first hidden layer and the second hidden layer are connected through a ReLU activation function, the second hidden layer and the output layer are connected through a tanh activation function, and 1 node is arranged in the output layer and used for outputting calculated while-drilling pulse neutron porosity.
And 6.2, setting a porosity precision value of the pulse neutrons while drilling.
Step 6.3, randomly selecting a pulse neutron porosity while drilling sample in the training set, inputting the pulse neutron porosity while drilling sample into the multilayer perception neural network, enabling nodes in the input layer of the multilayer perception neural network to obtain input parameters in the pulse neutron porosity while drilling sample, activating the input parameters through a ReLU function, and transmitting the activated input parameters to the first hidden layer, wherein the output result of the first hidden layer is shown as a formula (2):
Figure DEST_PATH_IMAGE010A
(2)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE012A
as a result of the output of the first hidden layer,
Figure DEST_PATH_IMAGE014A
is an input parameter of the multi-layer perceptive neural network,
Figure DEST_PATH_IMAGE016A
is the size of a drill collar,
Figure DEST_PATH_IMAGE018A
the ratio of the counts of the thermal neutrons to the iron marker,
Figure DEST_PATH_IMAGE020A
the counting is carried out for the non-bullet gamma counting,
Figure DEST_PATH_IMAGE022A
in order to capture the gamma counts,
Figure DEST_PATH_IMAGE024A
is the gamma counting of the characteristic of the chlorine element,
Figure DEST_PATH_IMAGE026A
is a macroscopic capture cross-section for the wellbore medium,
Figure DEST_PATH_IMAGE028A
is a macroscopic trapping cross section of the formation medium,
Figure DEST_PATH_IMAGE030A
is the weight of the first hidden layer,
Figure DEST_PATH_IMAGE032A
is the bias of the first hidden layer.
The output result of the first hidden layer is activated by the ReLU function and then transmitted to the second hidden layer, and the output result of the second hidden layer is shown in formula (3):
Figure 525891DEST_PATH_IMAGE034
(3)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE036A
as a result of the output of the second hidden layer,
Figure DEST_PATH_IMAGE038A
is the weight of the second hidden layer and,
Figure DEST_PATH_IMAGE040A
is the bias of the second hidden layer.
The output result of the second hidden layer is activated by the tanh function and then transmitted to the output layer, and the porosity of the while-drilling pulse neutrons output by the output layer is as follows:
Figure DEST_PATH_IMAGE042A
(4)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE044A
the porosity of the pulse neutron output by the output layer while drilling,
Figure DEST_PATH_IMAGE046A
as a weight of the output layer,
Figure DEST_PATH_IMAGE048A
is the bias of the output layer.
Step 6.4, calculating the formation porosity of the model by combining the porosity value of the while-drilling pulse neutrons corresponding to the while-drilling pulse neutron porosity sample according to the porosity of the while-drilling pulse neutrons output by the output layer, calculating the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network, and comparing the porosity error value with the preset precision value of the porosity of the while-drilling pulse neutrons; if the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network is larger than the set porosity precision value of the while-drilling pulse neutrons, updating the weight and bias of each layer in the multilayer perception neural network by using a formula (5) based on a logcosh loss function through a back propagation and random gradient descent method, returning to the step 6.3, and continuing to train the multilayer perception neural network; and if the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network is not less than the set porosity precision value of the while-drilling pulse neutrons, stopping training the multilayer perception neural network, and entering the step 7.
In this embodiment, as shown in fig. 6, when the training frequency of the multilayer perceptual neural network reaches 400 times, the loss function value is stable, the error value of the pulse-while-drilling neutron porosity of the multilayer perceptual neural network is not less than the set precision value of the pulse-while-drilling neutron porosity, and the calculation precision of the pulse-while-drilling neutron porosity of the multilayer perceptual neural network is stable.
And 7, verifying the accuracy of the porosity of the pulse-while-drilling neutrons by using the multi-layer perceptive neural network after training by using a verification set, inputting all the porosity of the pulse-while-drilling neutrons in the verification set into the multi-layer perceptive neural network after training, respectively calculating the porosity of the pulse-while-drilling neutrons by using the multi-layer perceptive neural network for each sample of the porosity of the pulse-while-drilling neutrons to obtain a calculated value of the porosity of the pulse-while-drilling neutrons, calculating the absolute error value of the porosity of the pulse-while-drilling neutrons by using a formula (6) in combination with the porosity of the pulse-while-drilling neutrons corresponding to the sample of the porosity of the pulse-while-drilling neutrons, and proving that the calculation result of the porosity of the pulse-while-drilling neutrons in the multi-layer perceptive neural network is accurate if the absolute error value is less than 1%.
In the embodiment, a distribution frequency diagram of absolute error values of concentrated while-drilling pulse neutron porosity samples is verified, as shown in fig. 7, the distribution frequency diagram can be obtained from fig. 7, the true value that the porosity of the while-drilling pulse neutrons is close to the porosity of the while-drilling pulse neutrons is obtained through calculation of the multilayer perception neural network trained by the method, the absolute error between the porosity of the while-drilling pulse neutrons calculated by the multilayer perception neural network after training and the true value of the formation porosity is basically smaller than 1%, and the accuracy of the multilayer perception neural network used for calculating the porosity of the while-drilling pulse neutrons after training by the method is verified.
And 8, based on the scale relation between the simulated value and the measured value of the pulse-while-drilling neutron porosity logging instrument established in the step 3, performing scale processing on a borehole medium macroscopic capture section, a formation medium macroscopic capture section, a thermal neutron count ratio, a non-elastic gamma count, a capture gamma count, an iron element characteristic gamma count and a chlorine element characteristic gamma count which are actually measured by the pulse-while-drilling neutron porosity logging instrument, marking the thermal neutron count ratio by using the iron element characteristic gamma count, inputting the borehole medium macroscopic capture section, the formation medium macroscopic capture section, the non-elastic gamma count, the capture gamma count, the chlorine element characteristic gamma count, the iron mark thermal neutron count ratio and the drill collar size of the pulse-while-drilling neutron porosity logging instrument after scale processing into a multilayer perception neural network, and calculating by using the multilayer perception neural network to obtain the pulse-while-drilling neutron porosity.
It is to be understood that the above description is not intended to limit the present invention, and the present invention is not limited to the above examples, and those skilled in the art may make modifications, alterations, additions or substitutions within the spirit and scope of the present invention.

Claims (8)

1. An intelligent processing method for porosity of pulse neutrons while drilling based on iron neutron marking is characterized in that a pulse neutron porosity logging while drilling instrument is adopted for measurement, and comprises the following steps:
1) Placing a pulse neutron porosity logging-while-drilling instrument in a standard calibration well for measurement to obtain an actual measurement value of the pulse neutron porosity logging-while-drilling instrument;
2) Establishing a while-drilling pulse neutron porosity numerical calculation model by using a Monte Carlo numerical simulation method to obtain a simulation value of a while-drilling pulse neutron porosity logging instrument;
3) Establishing a scale relation between the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument based on the analog value and the measured value of the pulse-while-drilling neutron porosity logging instrument;
4) Changing the size of a drill collar of a pulse-while-drilling neutron porosity logging instrument, sequentially changing the lithology, the water mineralization of the stratum, the gap of the drill collar, the mud type and the porosity of the stratum under different drill collar sizes, reestablishing a pulse-while-drilling neutron porosity value calculation model, and simulating to obtain the size of the drill collar, the macroscopic capture section of the borehole medium, the macroscopic capture section of the stratum medium, the counting ratio of iron-marked thermal neutrons, the non-elastic gamma counting, the capture gamma counting and the chlorine element characteristic gamma counting under different logging-while-drilling environments to form a plurality of groups of pulse-while-drilling neutron porosity samples;
5) Constructing a sample database based on a plurality of groups of while-drilling pulse neutron porosity samples, and dividing the sample database into a training set and a verification set;
6) Constructing a multilayer perception neural network, and calculating the porosity of the pulse neutrons while drilling by using a training set to train the multilayer perception neural network;
7) Verifying the accuracy of the porosity of the neutrons in the while-drilling pulse by using the multi-layer perception neural network after training by using a verification set;
8) Based on the scale relation between the analog value and the measured value of the pulse neutron porosity logging-while-drilling instrument, after a borehole medium macroscopic capture section, a formation medium macroscopic capture section, a non-elastic gamma count, a capture gamma count, a chlorine element characteristic gamma count and an iron-marked thermal neutron count ratio which are actually measured by the pulse neutron porosity logging-while-drilling instrument are obtained through scale processing, the obtained values and the size of the drill collar are input into a multilayer perception neural network together, and the pulse neutron porosity-while-drilling is obtained through calculation;
in the step 4), the iron-labeled thermal neutron counting ratio is obtained by labeling the thermal neutron counting ratio by iron element characteristic gamma counting, as shown in formula (1):
Figure DEST_PATH_IMAGE002
(1)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE004
the counting ratio of the thermal neutrons is marked for the iron,
Figure DEST_PATH_IMAGE006
is the counting ratio of the thermal neutrons,
Figure DEST_PATH_IMAGE008
and (4) carrying out gamma counting on the iron element characteristics after normalization treatment.
2. The intelligent processing method for porosity of pulse-while-drilling neutrons based on iron neutron marking as claimed in claim 1, wherein the pulse-while-drilling neutron porosity logging instrument is arranged in a drill collar, a mud diversion channel is arranged in the drill collar, the outer wall of a shell of the pulse-while-drilling neutron porosity logging instrument is closely attached to the inner wall of the drill collar, and a far neutron detector, a gamma detector, a near neutron detector, a ferrotungsten shield and a pulse neutron source are arranged in the shell at intervals from top to bottom.
3. The intelligent processing method for porosity of pulse-while-drilling neutrons based on iron neutron marking as claimed in claim 2, wherein the near neutron detector and the far neutron detector are both in a cylinder structure, the diameter of the near neutron detector is set to be 4cm, the length of the near neutron detector is set to be 5cm, the diameter of the far neutron detector is set to be 6cm, the length of the far neutron detector is set to be 10cm, the diameter of the gamma detector is set to be 4cm, the length of the gamma detector is set to be 4cm, the source distance of the near neutron detector is set to be 35cm, the source distance of the far neutron detector is set to be 70cm, and the source distance of the gamma detector is set to be 50cm.
4. The method as claimed in claim 1, wherein in the step 3), the calibration relationship between the simulated value and the measured value of the LWD tool includes a calibration relationship between the simulated value and the measured value of the macroscopic capture cross-section of the borehole medium, a calibration relationship between the simulated value and the measured value of the macroscopic capture cross-section of the formation medium, a calibration relationship between the simulated value and the measured value of the thermal neutron count ratio, a calibration relationship between the simulated value and the measured value of the non-bomb gamma count, a calibration relationship between the simulated value and the measured value of the capture gamma count, a calibration relationship between the simulated value and the measured value of the iron element gamma count, and a calibration relationship between the simulated value and the measured value of the chlorine element gamma count.
5. The intelligent iron neutron marker-based pulsed neutron porosity while drilling processing method as claimed in claim 1, wherein in step 4), the drill collar sizes of the pulsed neutron porosity while drilling logging tool are set to be 12.065cm, 17.145cm and 20.32cm in sequence, and the lithology of the formation medium is set to be sandstone, limestone, dolomite, sand-lime mixed rock and sand-dolomite mixed rock in sequence, wherein the sand-lime mixed rock is formed by mixing sandstone and limestone, the sand-lime mixed rock is formed by mixing sandstone and dolomite, the water mineralization of the formation is set to be 0kppm, 5kppm, 10kppm, 20kppm, 50kppm, 100kppm and 200kppm in sequence, the drill collar gaps are set to be 0cm, 0.5cm, 1cm, 1.5cm and 2cm in sequence, the mud types comprise oil-based mud and water-based mud, and the formation porosity is set to be 0, 5%, 10%, 15%, 20%, 25%, 30%, 35%, 40%, 45% and 50% in sequence.
6. The intelligent processing method for porosity of pulsed neutrons while drilling based on iron neutron marking as claimed in claim 1, wherein the step 6) comprises:
6.1, constructing a multilayer perception neural network, wherein the multilayer perception neural network is provided with an input layer, a first hidden layer, a second hidden layer and an output layer, the input layer, the first hidden layer, the second hidden layer and the output layer are sequentially connected, 7 nodes are arranged in the input layer and are used for obtaining input parameters of the multilayer perception neural network, the input parameters comprise the size of a drill collar in a pulse neutron porosity while drilling sample, a macroscopic capture section of a borehole medium, a macroscopic capture section of a formation medium, a counting ratio of iron-marked thermal neutrons, a non-elastic gamma count, a capture gamma count and a chlorine element characteristic gamma count, 10 nodes are arranged in the first hidden layer, 15 nodes are arranged in the second hidden layer, the first hidden layer and the second hidden layer are connected through a ReLU activation function, the second hidden layer and the output layer are connected through a tanh activation function, and 1 node is arranged in the output layer and is used for outputting the calculated pulse neutron porosity while drilling;
step 6.2, setting a porosity precision value of the pulse neutrons while drilling;
step 6.3, randomly selecting a pulse neutron porosity while drilling sample in the training set, inputting the pulse neutron porosity while drilling sample into the multilayer perception neural network, enabling nodes in the input layer of the multilayer perception neural network to obtain input parameters in the pulse neutron porosity while drilling sample, activating the input parameters through a ReLU function, and transmitting the activated input parameters to the first hidden layer, wherein the output result of the first hidden layer is shown as a formula (2):
Figure DEST_PATH_IMAGE010
(2)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE012
as a result of the output of the first hidden layer,
Figure DEST_PATH_IMAGE014
is an input parameter of the multi-layer perceptive neural network,
Figure DEST_PATH_IMAGE016
is the size of a drill collar,
Figure DEST_PATH_IMAGE018
the ratio of the counts of the thermal neutrons to the iron marker,
Figure DEST_PATH_IMAGE020
is a non-elastic gamma meterThe number of the first and second groups is counted,
Figure DEST_PATH_IMAGE022
in order to capture the gamma counts,
Figure DEST_PATH_IMAGE024
is the gamma counting of the characteristic of the chlorine element,
Figure DEST_PATH_IMAGE026
for the macroscopic capture cross-section of the wellbore medium,
Figure DEST_PATH_IMAGE028
is a macroscopic trapping cross section of the formation medium,
Figure DEST_PATH_IMAGE030
is the weight of the first hidden layer,
Figure DEST_PATH_IMAGE032
a bias for the first hidden layer;
the output result of the first hidden layer is activated by the ReLU function and then transmitted to the second hidden layer, and the output result of the second hidden layer is shown in formula (3):
Figure DEST_PATH_IMAGE034
(3)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE036
as a result of the output of the second hidden layer,
Figure DEST_PATH_IMAGE038
is the weight of the second hidden layer,
Figure DEST_PATH_IMAGE040
a bias for the second hidden layer;
and the output result of the second hidden layer is activated by the tanh function and then transmitted to the output layer, and the porosity of the pulse-while-drilling neutrons output by the output layer is as follows:
Figure DEST_PATH_IMAGE042
(4)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE044
the porosity of the pulse neutron output by the output layer while drilling,
Figure DEST_PATH_IMAGE046
as a weight of the output layer,
Figure DEST_PATH_IMAGE048
is the bias of the output layer;
step 6.4, calculating the formation porosity of the model by combining the porosity value of the while-drilling pulse neutrons corresponding to the while-drilling pulse neutron porosity sample according to the porosity of the while-drilling pulse neutrons output by the output layer, calculating the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network, and comparing the porosity error value with the preset precision value of the porosity of the while-drilling pulse neutrons; if the porosity error value of the while-drilling pulse neutrons of the multilayer perception neural network is larger than the set porosity precision value of the while-drilling pulse neutrons, updating the weight and bias of each layer in the multilayer perception neural network by using a logcosh loss function and adopting a back propagation and random gradient descent method, returning to the step 6.3, and continuing to train the multilayer perception neural network; and if the while-drilling pulse neutron porosity error value of the multilayer perception neural network is not less than the set while-drilling pulse neutron porosity precision value, stopping training the multilayer perception neural network, and entering the step 7).
7. The intelligent processing method for porosity of pulsed neutrons while drilling based on iron neutron marking as claimed in claim 6, wherein in step 6.4, the weight and bias updating formula of each layer in the multilayer perceptive neural network is as follows:
Figure DEST_PATH_IMAGE050
(5)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE052
to minimize the value of the variable at which the logcosh loss function takes a minimum,
Figure DEST_PATH_IMAGE054
as logcosh loss function.
8. The method for intelligently processing porosity of pulse-while-drilling neutrons based on iron neutron labeling according to claim 1, wherein in step 7), all porosity of pulse-while-drilling neutrons in the validation set are input into the trained multilayer perceptive neural network, a calculation value of the porosity of the pulse-while-drilling neutrons is obtained by using the multilayer perceptive neural network for each sample of the porosity of the pulse-while-drilling neutrons, and an absolute error value of the porosity of the pulse-while-drilling neutrons is calculated by using the multilayer perceptive neural network by combining the formation porosity of the calculation model of the porosity value of the pulse-while-drilling neutrons corresponding to the sample of the porosity of the pulse-while-drilling neutrons, and if the absolute error value is less than 1%, the calculation result of the porosity of the pulse-while-drilling neutrons in the multilayer perceptive neural network is proved to be accurate;
the calculation formula of the absolute error value of the porosity of the while-drilling pulse neutrons is as follows:
Figure DEST_PATH_IMAGE056
(6)
in the formula (I), the compound is shown in the specification,
Figure DEST_PATH_IMAGE058
is the absolute error value of the porosity of the neutron in the while-drilling pulse,
Figure DEST_PATH_IMAGE060
to verify the number of concentrated while-drilling pulsed neutron porosity samples,
Figure DEST_PATH_IMAGE062
in order to verify the serial number of the porosity sample of the concentrated while-drilling pulse neutron,
Figure DEST_PATH_IMAGE064
for multi-layer perceptive neural network according to
Figure DEST_PATH_IMAGE062A
A pulse-while-drilling neutron porosity value calculated for each pulse-while-drilling neutron porosity sample,
Figure DEST_PATH_IMAGE066
is a first
Figure DEST_PATH_IMAGE062AA
And calculating the formation porosity of the model according to the pulse-while-drilling neutron porosity value corresponding to the pulse-while-drilling neutron porosity sample.
CN202211194628.0A 2022-09-29 2022-09-29 While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking Active CN115291288B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211194628.0A CN115291288B (en) 2022-09-29 2022-09-29 While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211194628.0A CN115291288B (en) 2022-09-29 2022-09-29 While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking

Publications (2)

Publication Number Publication Date
CN115291288A CN115291288A (en) 2022-11-04
CN115291288B true CN115291288B (en) 2022-12-30

Family

ID=83834251

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211194628.0A Active CN115291288B (en) 2022-09-29 2022-09-29 While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking

Country Status (1)

Country Link
CN (1) CN115291288B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116500694B (en) * 2023-06-28 2023-09-01 中海油田服务股份有限公司 Post-sleeve physical quantity inversion method, post-sleeve physical quantity inversion device, computing equipment and storage medium

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4737636A (en) * 1986-11-17 1988-04-12 Halliburton Company Simultaneous neutron lifetime and oxygen activation and fluid sampling methods and apparatus to locate fluid flow in a cased well borehole
CN115012920A (en) * 2022-06-08 2022-09-06 中国石油大学(华东) Controllable neutron source multi-spectrum logging instrument and method based on double CLYC double-particle detectors

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4625111A (en) * 1985-02-04 1986-11-25 Halliburton Company Method and apparatus for pulsed neutron capture measurement of borehole conditions and formation hydrogen index
US5345077A (en) * 1991-07-24 1994-09-06 Mobil Oil Corporation Method and apparatus for producing a porosity log of a subsurface formation corrected for detector standoff
US9568638B2 (en) * 2010-06-30 2017-02-14 Schlumberger Technology Corporation Identification of neutron capture from a pulsed neutron logging tool
US9052405B2 (en) * 2013-05-15 2015-06-09 Schlumberger Technology Corporation Borehole fluid effect correction for pulsed neutron porosity measurements
CN103696765B (en) * 2013-11-06 2016-08-17 中国石油大学(华东) Double LaBr based on controllable neutron source3detector elements spectrometer tool and logging method
CN104747179A (en) * 2013-12-31 2015-07-01 中国石油化工集团公司 Stratum density measuring while drilling instrument based on deuterium-tritium accelerator neutron source
CN107462929B (en) * 2017-07-25 2019-04-26 中国石油大学(华东) Cupro-nickel mineral products level measuring arrangement and method in a kind of well
CN110439545B (en) * 2019-08-02 2023-04-25 中国石油天然气集团有限公司 Environment correction method for logging instrument with neutron porosity of controllable source while drilling
CN112904436B (en) * 2021-04-14 2022-07-05 中国石油大学(华东) Porosity measurement method combining element yield and thermal neutron count ratio
CN115061219B (en) * 2022-08-17 2022-11-08 北京派特杰奥科技有限公司 Fracture type reservoir prediction identification method and system based on petroleum and natural gas detection

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4737636A (en) * 1986-11-17 1988-04-12 Halliburton Company Simultaneous neutron lifetime and oxygen activation and fluid sampling methods and apparatus to locate fluid flow in a cased well borehole
CN115012920A (en) * 2022-06-08 2022-09-06 中国石油大学(华东) Controllable neutron source multi-spectrum logging instrument and method based on double CLYC double-particle detectors

Also Published As

Publication number Publication date
CN115291288A (en) 2022-11-04

Similar Documents

Publication Publication Date Title
US8476584B2 (en) Methods for sourceless density downhole measurement using pulsed neutron generator
US6032102A (en) Method and apparatus for measuring well characteristics and formation properties
US9268056B2 (en) Neutron porosity based on one or more gamma ray detectors and a pulsed neutron source
CN105874354A (en) System and method for making downhole measurements
US10725201B2 (en) Compensated neutron correction for contributions outside the petrophysical model
CN108643890B (en) Method for determining formation porosity
US6150655A (en) Inferential measurement of photoelectric absorption cross-section of geologic formations from neutron-induced, gamma-ray spectroscopy
Wahl et al. The thermal neutron decay time log
CN103375166A (en) Method for determining content of a plurality of mineral constituents in stratum
CN102084271B (en) Absolute elemental concentrations from nuclear spectroscopy
MX2013006661A (en) Gas detection and quantification method using a pulsed neutron logging tool.
CN112904436B (en) Porosity measurement method combining element yield and thermal neutron count ratio
CN115291288B (en) While-drilling pulse neutron porosity intelligent processing method based on iron neutron marking
CN106250619B (en) Method and device for determining mineral content of stratum
Wang et al. Neutron transport correction and density calculation in the neutron-gamma density logging
Liu et al. A method to improve the sensitivity of neutron porosity measurement based on DT source
Eyvazzadeh et al. Modern carbon/oxygen logging methodologies: comparing hydrocarbon saturation determination techniques
Fu et al. A new method of gas reservoir evaluation based on neutron cross section logging
CN110469324A (en) A kind of calculating density of earth formations method based on pulsed neutron log
CN115267930A (en) High-sensitivity neutron porosity measurement method based on D-T pulse neutron source
CN100492055C (en) A chlorine spectrometry logging method
CN115012920A (en) Controllable neutron source multi-spectrum logging instrument and method based on double CLYC double-particle detectors
US5272629A (en) Method for determining the slowing down length and the porosity of a formation surrounding a borehole
RU152169U1 (en) Borehole Device with Neutral Measuring Probes
CN111638559B (en) Fast neutron scattering cross section characterization method based on pulse neutron logging

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