Saturation multi-probe instrument nuclear calculation and processing automation method and system, terminal and medium

Abstract

The invention discloses a saturation multi-probe instrument nuclear calculation and processing automation method, a saturation multi-probe instrument nuclear calculation and processing automation system, a saturation multi-probe instrument nuclear calculation and processing automation terminal and a saturation multi-probe instrument nuclear calculation and processing automation medium, and belongs to the technical field of nuclear logging in oil-gas exploration. Data generated by the model is extracted, converted and integrated into a simulation spectrum database, and saturation parameters of the simulation data are automatically generated into the simulation spectrum database; meanwhile, saturation actual measurement logging data are preprocessed, saturation parameter results are automatically calculated, the water-containing saturation, the gas-containing saturation and the oil-containing saturation under each depth are solved, and the whole system basically achieves the automatic processing flow; according to the method, the dependence on manual processing is reduced, the consumption of human resources is reduced, the influence of human subjective judgment on results is reduced, the consistency of the results in the oil-gas exploration field is improved, important saturation parameters for oil-gas exploration are generated in batches, and the intelligent level in the oil-gas field is improved.

Description

Technical Field

[0001] The present invention belongs to the technical field of nuclear logging in oil and gas exploration, and in particular relates to a method, system, terminal and medium for saturation multi-probe instrument nuclear calculation and processing automation. Background Art

[0002] Pulsed neutron logging is an important nuclear logging technique. A pulsed neutron source emits neutrons into the formation surrounding the wellbore. Detectors receive the radiation signals generated by the neutrons interacting with the formation material to determine formation saturation information. This technique was initially used to measure parameters such as the macroscopic thermal neutron absorption cross section (CCS) or the neutron-induced gamma-ray spectrum to distinguish between brine and hydrocarbons, thereby calculating water saturation in the formation pores. Typical pulsed neutron logging methods include neutron lifetime logging (also known as thermal neutron decay time Sigma (Σ) logging) and carbon-oxygen ratio logging (C / O), which measures the ratio of gamma ray production of carbon to oxygen. The Σ method is the oldest and most traditional saturation calculation method. Given a known mineralization in the formation pores, the Σ value is calculated based on the decay rate of thermal neutron counts or captured gamma ray counts over time. Formations containing brine have higher Σ values, while those containing oil and gas have lower Σ values. Combining the Σ curves obtained from logging with calculated parameters such as porosity, lithologic profiles, and shale content, the water saturation and oil and gas saturation in the formation can be calculated using material balance principles and weighted methods. Pulsed neutron logging can quantitatively characterize formation fluid saturation after well completion, playing an irreplaceable role in reservoir dynamic monitoring and remaining oil evaluation. For example, repeated pulsed neutron logging can monitor the changing trends of water saturation during waterflooding development in an oilfield, thereby determining oil recovery efficiency and remaining oil distribution. This logging information is widely used in reservoir monitoring, production and injection stimulation, enhanced oil recovery, and the discovery of missed remaining oil reservoirs, and is of great significance to oil and gas field development.

[0003] At present, when using pulsed neutron logging data to calculate formation saturation, the industry often uses some traditional methods to process data and calculate parameters. The main methods include: (1) Chart comparison method: Based on the standard interpretation chart provided by the logging service company, the logging instrument readings (such as Σ value or carbon-oxygen ratio, etc.) are compared with the curves or cross-plots on the chart, and the corresponding water saturation values ​​are manually found. Engineers need to find the closest solution on the chart based on the shape and reading of the logging curve and read the corresponding saturation results. (2) Manual regression method: Interpreters usually perform regression correction on the logging response and actual saturation based on limited formation core analysis or oil test and production data to obtain an empirical formula suitable for the target well. (3) Excel calculation: The pulsed neutron logging data is exported to a spreadsheet, and the engineer calculates it in Excel according to the saturation formula. This method requires manual setting of relevant formulas and parameters, and substituting the logging values ​​of each layer into the formula to calculate the water saturation. The above traditional methods have been widely used for quite a long time in the past, and interpreters use these methods to evaluate and report the saturation of single wells.

[0004] Although the above traditional methods can complete the calculation of saturation to a certain extent, they also have obvious shortcomings and limitations: (1) Dependence on manual labor and low efficiency: Whether it is chart comparison or manual regression, they are highly dependent on manual operation and expert experience. Interpreters need to manually compare charts or adjust parameters well by well and layer by layer, which is extremely time-consuming and inefficient when processing data from a large number of wells. (2) Accuracy and consistency are affected by human factors: Manual reading of charts may cause errors due to personal judgment, and different interpreters may have different interpretations of the same data, resulting in difficulty in ensuring the accuracy and consistency of the calculation results. The adjustment of parameters in the manual regression process is also subjective and easy to introduce bias. Overall, the manual process increases the possibility of accumulated operational errors. (3) It is difficult to adapt to batch and multi-well data processing: With the increase in the number of wells in oil field development and the increase in the amount of logging data, it is difficult to meet the needs by using Excel or charts for calculation on a well-by-well basis. When batch processing of multi-well data is required, traditional methods lack effective batch calculation functions and are difficult to quickly process massive data, which restricts the overall efficiency of reservoir evaluation. (4) Lack of intelligent and automated methods: Existing methods rely primarily on manual experience and lack the assistance of computer intelligent algorithms. They are unable to adaptively optimize saturation calculation models based on different reservoir and formation conditions. For example, when the formation water salinity or lithology changes, the current process requires manual re-adjustment of parameters and lacks automatic adjustment and learning capabilities. This non-intelligent approach is no longer able to cope with the complex and changing needs of reservoir logging interpretation. Summary of the Invention

[0005] In response to the problems existing in the above-mentioned traditional data processing and calculation methods, the present invention provides an automated method, system, terminal and medium for saturation multi-probe instrument core calculation and processing, which can automatically convert the data signal of pulse neutron logging directly into formation saturation results, reduce manual intervention, and improve processing efficiency and accuracy. It is of great significance to improving the efficiency and quality of oil reservoir logging interpretation, and can effectively make up for the shortcomings of current technical means and meet the needs of modern oil field development for fast, accurate and batch logging interpretation.

[0006] In order to achieve the above purpose, the technical methods adopted by the present invention are as follows:

[0007] The method for automating calculation and processing of saturation multi-probe instrument cores includes the following steps:

[0008] Step 1: Through the user interface, select the instrument parameters of the nuclear physics logging instrument to build a nuclear physics logging instrument model. Then, select the formation information, well parameters, and nuclear physics simulation parameters to automatically generate a nuclear logging simulation batch model. The model is uploaded to the high-performance computing platform for simulation calculation and generates a nuclear physics simulation output file.

[0009] Step 2: Determine whether the nuclear physics simulation output file includes energy spectrum and time spectrum data. If so, directly extract the energy spectrum and time spectrum data and integrate them into the batch processing template file; otherwise, convert the flux data in the nuclear physics simulation output file into a gamma spectrum and integrate it into the batch processing template file; and then generate a simulation spectrum database in the batch processing template file;

[0010] The energy spectrum includes inelastic gamma ray spectrum, capture gamma ray spectrum, thermal neutron and epithermal neutron counting spectrum;

[0011] Step 3: Quickly load the simulated spectrum database and select the saturation parameters to be calculated through the user interface, including carbon-oxygen ratio (C / O), Sigma, fast neutron scattering cross section (FNXS), multi-detector inelastic gamma count ratio (MultiINE), and multi-detector captured gamma count ratio (MultiCAP);

[0012] Step 4: Select the calculation algorithm corresponding to the selected saturation parameter through the user interface, automatically generate the saturation parameter results and the corresponding fitting curve plate results, and save them into the batch processing template file, which together with the simulation spectrum database serves as the saturation simulation database;

[0013] Step 5: Quickly load the saturation logging data, including the total spectrum of the burst period, the energy spectrum of the captured ray, the time spectrum, the formation porosity and lithology;

[0014] Step 6: Preprocess the saturation well logging data to obtain corrected saturation well logging data;

[0015] Step 7: Quickly load the saturation simulation database and the corrected saturation well data, select the saturation parameters to be calculated and their corresponding calculation algorithms through the user interface, generate saturation parameter results for the saturation well data at various well depths, and calculate the water saturation, gas saturation, or oil saturation for the saturation well data at various well depths based on the corresponding fitting curve chart results in the saturation simulation database.

[0016] Furthermore, the formation information described in step 1 includes lithology and pore fluid properties, the well parameters include instrument position, casing parameters, number of tubing layers, cement sheath parameters and wellbore fluid parameters, and the nuclear physics simulation parameters include the measurement mode, particle number and thread number of inelastic ray energy spectrum and capture ray energy spectrum.

[0017] Furthermore, the lithology includes conventional lithology formations (sandstone, limestone, dolomite) and complex lithology formations (based on the rock skeleton density and the mass ratio of each element in the rock skeleton input by the user), and the pore fluid properties include the medium type and density of water, oil and gas.

[0018] Furthermore, the nuclear well logging simulation batch processing model is a Monte Carlo nuclear physics simulation batch processing model.

[0019] Furthermore, in step 2, the flux data is converted into a gamma spectrum using the detector response function matrix, as follows:

[0020] Y=XA

[0021] Where, represents flux data; represents the gamma spectrum; It represents the detector response function matrix, which is calculated by simulating the detector performance of nuclear physics logging instruments using nuclear physics simulation software; n is the number of channels in the energy spectrum.

[0022] Furthermore, the calculation model for FNXS in step 4 is:

[0023] FNXS i =ω6×ln(N I3 )+ω7

[0024]

[0025] Where m represents the total number of lithologies or minerals contained in the formation; FNXS i is the comprehensive FNXS of the formation under the i-th lithology or mineral; N I3represents the total counts in all windows of the non-elastic ray spectrum of the ultra-long-range detector under the i-th lithology or mineral; ln(·) represents the natural logarithm; ω6 and ω7 represent the inversion coefficients of the comprehensive FNXS of the formation under the i-th lithology or mineral, which are obtained by calculating the relationship between the counts in the non-elastic ray spectrum window of the ultra-long-range detector under the i-th lithology or mineral and the FNXS based on the simulated spectrum database; V i It represents the volume proportion of the i-th lithology or mineral in the formation.

[0026] Furthermore, the specific process of preprocessing in step 6 is:

[0027] Step 6.1: perform depth averaging and filtering on the total spectrum of the burst period and the captured ray energy spectrum in the saturation logging data, and then perform filtering at each depth;

[0028] Step 6.2: Based on the total spectrum of the burst period and the captured ray spectrum obtained in step 6.1, obtain the inelastic ray spectrum;

[0029] Step 6.3: Correct the inelastic ray energy spectrum and the capture ray energy spectrum to obtain a corrected inelastic gamma ray energy spectrum and a corrected capture gamma ray energy spectrum.

[0030] Furthermore, the fitting curve plate results in step 4 and step 7 include water saturation fitting curve plate results, gas saturation fitting curve plate results, and oil saturation fitting curve plate results.

[0031] Furthermore, the specific process of calculating the gas saturation at each well depth from the saturation logging well data in step 7 is as follows:

[0032] The saturation parameters required for calculation include FNXS, MultiINE and MultiCAP. The FNXS saturation parameter result corresponding to the gas-saturated formation in the saturation simulation database is FNXS. g The FNXS saturation parameter corresponding to the water-saturated formation is FNXS w The MultiINE saturation parameter result for the gas-saturated formation is g The MultiINE saturation parameter result for the water-saturated formation is w The MultiCAP saturation parameter result for gas-saturated formations is MultiCAP g The MultiCAP saturation parameter result corresponding to the water-saturated formation is MultiCAP w , and the FNXS saturation parameter result corresponding to the saturation logging data is recorded as FNXS meas , the MultiINE saturation parameter result is MultiINE meas, the MultiCAP saturation parameter result is MultiCAP meas ;

[0033] Calculate the FNXS gas saturation S corresponding to the saturation logging data g1 、MultiINE gas saturation S g2 and MultiCAP gas saturation S g3 ;

[0034]

[0035] Then the gas saturation S is calculated g :

[0036] S g =ω1×S g1 +ω2×S g2 +ω3×S g3

[0037] ω1+ω2+ω3=1

[0038] Where ω1, ω2 and ω3 are weight coefficients.

[0039] Furthermore, FNXS g 、MultiINE g and MultiCAP g The greater the slope of the corresponding fitting curve chart result, the greater the corresponding weight coefficient value.

[0040] Furthermore, the fast loading refers to loading by adopting multi-threaded asynchronous loading, data format conversion, memory mapping, multi-process independent processing, enhanced data caching and the like.

[0041] The saturation multi-probe instrument nuclear calculation and processing automation system is used to implement the saturation multi-probe instrument nuclear calculation and processing automation method, specifically including a nuclear physics simulation module, a simulation spectrum database generation module, a saturation parameter and algorithm selection module, a saturation simulation database generation module, a preprocessing module and a saturation calculation module; wherein:

[0042] The nuclear physics simulation module is used to automatically generate a nuclear well logging simulation batch model, upload it to a high-performance computing platform for simulation to generate a nuclear physics simulation output file, and output it to a simulation spectrum database generation module;

[0043] The simulation spectrum database generation module is used to extract the energy spectrum and time spectrum data in the nuclear physics simulation output file and integrate them into the batch processing template file, or convert the flux data into the gamma energy spectrum and integrate it into the batch processing template file to generate a simulation spectrum database, which is output to the saturation parameter selection module in the simulation mode;

[0044] The saturation parameter selection module is used to select the saturation parameters to be calculated and their corresponding calculation algorithms in the simulation mode or the measurement mode, and is connected to the saturation simulation database generation module;

[0045] The saturation simulation database generation module is used to automatically execute the saturation parameter calculation algorithm in the simulation mode, save the saturation parameter results calculated based on the simulation spectrum database into a batch processing template file, and use it together with the simulation spectrum database as the saturation simulation database. The calculation results corresponding to the saturation well logging data are output to the saturation calculation module.

[0046] The preprocessing module is used to preprocess the loaded saturation well logging data to obtain the corrected saturation well logging data and output it to the saturation parameter selection module in the measurement mode;

[0047] The saturation calculation module is used to calculate the water saturation, gas saturation or oil saturation at each well depth based on the saturation logging well data.

[0048] A terminal includes a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the saturation multi-probe instrument core calculation and processing automation method.

[0049] A computer-readable storage medium stores a computer program which, when executed by a processor, implements the steps of the saturation multi-probe instrument core calculation and processing automation method.

[0050] Compared with the prior art, the present invention has the following beneficial effects:

[0051] The present invention proposes an automated method, system, terminal, and medium for calculating and processing saturation multi-probe instruments. A nuclear logging simulation batch processing model is generated by user input of formation information, wellbore information, and nuclear physics simulation parameters. The data generated by the model is extracted, processed, converted, and integrated into a simulation spectrum database. The simulation data is then automatically processed for saturation parameters to generate a simulation spectrum database. Simultaneously, the measured saturation logging data is preprocessed, and the saturation parameter results are automatically calculated to obtain water, gas, and oil saturations at each depth. The entire system essentially implements an automated processing flow.

[0052] The present invention is beneficial to reducing dependence on manual processing, reducing the consumption of human resources, reducing the impact of human subjective judgment on the results, improving the consistency of results in the oil and gas exploration field, batch-generating important saturation parameters for oil and gas exploration, and improving the level of intelligence in the oil and gas field. BRIEF DESCRIPTION OF THE DRAWINGS

[0053] In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the embodiments or the description of the prior art. Obviously, the drawings described below are only some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0054] Figure 1 This is a flow chart of the automated method for calculating and processing saturation multi-probe instrument cores proposed in Example 1;

[0055] Figure 2 Schematic diagram of the software architecture for executing the automated method for calculating and processing saturation multi-probe instrument cores in Example 1;

[0056] Figure 3 1 is a schematic diagram of the main interface of the software for executing the automated method for calculating and processing saturation multi-probe instrument cores in Example 1;

[0057] Figure 4 1 is a schematic diagram of the parameter setting interface of the nuclear well logging simulation batch processing model in Example 1;

[0058] Figure 5 This is a schematic diagram of the software interface for loading the simulation spectrum database, selecting saturation parameters, and calculating saturation parameter results in Example 1;

[0059] Figure 6 This is a schematic diagram of the preprocessing interface for saturation well logging data in Example 1;

[0060] Figure 7 Schematic diagram of the software interface for saturation calculation in Example 1. DETAILED DESCRIPTION

[0061] In order to enable those skilled in the art to better understand the solutions of the present invention, the technical solutions in the embodiments of the present invention will be clearly and completely described below in conjunction with the accompanying drawings of the embodiments of the present invention. Obviously, 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 ordinary technicians in this field without making any creative efforts shall fall within the scope of protection of the present invention.

[0062] The terms "first," "second," and so on, used in the specification and claims of the present invention and in the accompanying drawings, are used to distinguish between different objects, not to describe a specific order. Furthermore, the terms "including," "having," and any variations thereof, are intended to cover non-exclusive inclusions. For example, a process, method, apparatus, product, or device comprising a series of steps or units is not limited to the listed steps or units but may optionally include steps or units not listed, or may optionally include other steps or units inherent to the process, method, product, or device. The terms "automatic" and "automated" refer to the use of computer programs, scripts, and other methods that only require input, such as processing user input, nuclear physics simulation results, nuclear physics simulation databases, saturation measurement logging data, and so on, to output the information required by the user. The term "window" in energy spectrum refers to a selected energy range, and in time spectrum refers to a selected time range. The terms "merge" and "integrate" refer to storing data obtained from different processing steps in the same file. The terms "script" and "algorithm" are computer programs that implement the desired functionality. "Saturation measured well data" refers to the data collected by the saturation instrument during deep-sea and deep-earth instrument operations, including energy spectrum, time spectrum, and porosity φ and lithology data calculated from the two spectra in addition to saturation. It also includes other data such as element test logging generated during the operation.

[0063] References herein to "embodiments" mean that a particular feature, structure, or characteristic described in connection with the embodiments may be included in at least one embodiment of the present invention. The appearance of this phrase in various places in the specification does not necessarily refer to the same embodiment, nor does it constitute a separate or alternative embodiment that is mutually exclusive of other embodiments. It is understood, both explicitly and implicitly, by those skilled in the art that the embodiments described herein may be combined with other embodiments.

[0064] In this application, the word "exemplary" is used to mean "serving as an example, illustration, or illustration." Any embodiment described in this application as "exemplary" is not necessarily to be construed as preferred or advantageous over other embodiments. The following description is given to enable any person skilled in the art to implement and use the present application. In the following description, details are listed for the purpose of explanation. It should be understood that one of ordinary skill in the art can recognize that the present application can be implemented without using these specific details. In other instances, well-known structures and processes are not elaborated in detail to avoid obscuring the description of the present application with unnecessary details. Therefore, the present application is not intended to be limited to the embodiments shown, but is consistent with the widest scope consistent with the principles and features disclosed in this application.

[0065] It should be noted that since the method of the embodiment of the present application is executed in a computer device, the processing objects of each computer device exist in the form of data or information. For example, the saturation measured logging data is actually the corresponding data information measured by the saturation instrument in the deep sea and deep earth. It can be understood that the nuclear physics simulation database, saturation parameter results, standard energy spectrum, inelastic ray energy spectrum, capture ray energy spectrum, time spectrum, etc. in the subsequent embodiments are all corresponding data and information that exist for processing by the computer device, and the details will not be repeated here.

[0066] The present invention discloses an automated method, system, terminal, and medium for calculating and processing saturation data for multi-probe instruments. These methods, systems, and media help reduce reliance on manual processing, minimize human resource consumption, and mitigate the impact of subjective judgment on results. These methods improve the consistency of results in oil and gas exploration, enable batch generation of important saturation parameters for oil and gas exploration, and enhance the level of intelligence in the oil and gas field. These methods are described in detail below.

[0067] Example 1

[0068] This embodiment proposes an automated method for saturation multi-probe instrument core calculation and processing, which is implemented based on the system software installed on the computer hardware platform and the developed system environment. The system software installed on the computer hardware platform is the user interface for the automated calculation and processing proposed in this embodiment; the computer hardware platform can be a local server or a cloud server, which is not limited in this embodiment; the system environment is installed as a dependent environment required for the software interface to run, such as the code version required for Python / C++ / C# and the development environment relied on to implement scientific computing, drawing, data processing and other functions. This embodiment does not limit the implementation code of the interface used and the dependent development environment.

[0069] In this embodiment, the software architecture for executing the saturation multi-probe instrument core calculation and processing automation method is as follows: Figure 2 As shown, the main interface of the system software is as follows Figure 3 As shown in the figure, there are five main functional areas, among which the nuclear physics simulation model setting corresponds to step 1, the simulation spectrum database corresponds to step 2, the saturation parameter calculation corresponds to steps 3 and 4, the saturation measured logging data processing corresponds to step 6, and the measured formation saturation calculation corresponds to step 7.

[0070] The process of the saturation multi-probe instrument core calculation and processing automation method is as follows: Figure 1 As shown, the following steps are included:

[0071] Step 1: Through the user interface, select the instrument parameters of the nuclear physics logging instrument to build a nuclear physics logging instrument model. Then, select the formation information, well parameters, and Monte Carlo nuclear physics simulation parameters to automatically generate a Monte Carlo nuclear physics simulation batch model. The model is uploaded to the high-performance computing platform for simulation calculation and generates a nuclear physics simulation output file.

[0072] It should be noted that in actual situations, other parameters can be expanded according to actual needs. Figure 4 Typical parameters such as nuclear physics logging instrument model, formation lithology, pore fluid, well information, simulation information, etc. are given. For example, typical parameters of the nuclear physics logging instrument model include source information (source type, source direction, source intensity, source size), detector information (size, material), shield information (size, material) and the relative positions of each object. In the high-fidelity model, it also includes information such as the casing, circuit cavity, and boron sleeve. The instrument information is pre-stored in the batch file. The formation lithology includes conventional reservoir lithology (sandstone, limestone, dolomite) and complex lithology (the user can enter the rock skeleton density and the mass ratio of each element in the rock skeleton). The pore fluid includes the medium type and density of conventional water, oil, and gas. Monte Carlo nuclear physics simulation parameters (non-elastic / capture measurement mode, particle number, thread number), wellbore information (casing, number and thickness of tubing layers, cement sheath, wellbore fluid, instrument position). Among the typical parameters mentioned above, parameters such as multiple lithologies and multiple types of pore fluids can be simultaneously configured within the interface and automatically generated into a batch file to accommodate the complex and changing conditions in actual well logging environments. Due to the large number of simulation parameters, the simulation is relatively time-consuming and requires uploading to a high-performance computing platform for simulation. If time consumption is not a concern, uploading to any computer hardware platform equipped with Monte Carlo nuclear physics simulation is acceptable. By developing a first script, all user-entered information can be automatically generated into a batch template file.

[0073] Step 2: Determine whether the nuclear physics simulation output file includes energy spectrum and time spectrum data. If so, develop and execute a second script to locate the required data location in the nuclear physics output file, directly extract the energy spectrum and time spectrum data therein, and integrate them into the batch processing template file. The energy spectrum includes the inelastic gamma ray energy spectrum, the captured gamma ray energy spectrum, the thermal neutron and epithermal neutron counting energy spectrum, etc.; otherwise, add a conversion function between the extraction and integration functions, develop and execute a third script that implements the conversion function, convert the flux data in the nuclear physics simulation output file into a gamma energy spectrum, and integrate it into the batch processing template file; and then generate a simulation spectrum database in the batch processing template file;

[0074] As an example, the detector response function matrix is ​​used to convert the flux data into a gamma spectrum, and the calculation method is the same as the first calculation model:

[0075] Y=XA

[0076] Where, represents flux data; represents the gamma spectrum; n is the number of channels in the spectrum; represents the detector response function matrix, which is calculated by nuclear physics simulation software based on the detector performance of the nuclear physics well logging instrument. The specific acquisition process is as follows: in the nuclear physics simulation software, the energy spectrum data under each energy channel is simulated, specifically the flux data of the energy channel is set to 1, and the flux data of other energy channels is set to 0, and the detection performance of the nuclear physics well logging instrument is simulated to obtain it. Then, the energy spectrum data are arranged in the order of the energy channels to form the detector response function matrix A.

[0077] Step 3: Quickly load the simulated spectrum database and select the saturation parameters to be calculated through the user interface, including C / O, Sigma, FNXS, MultiINE, and MultiCAP;

[0078] The rapid loading of the simulation spectrum database adopted in this embodiment uses a data cache method (data objects are defined as abstract types) to speed up data loading;

[0079] It should be noted that C / O, Sigma, FNXS, MultiINE, and MultiCAP are typical parameters that need to be measured in the common measurement of water, oil, and gas saturation. In addition, parameters such as hydrogen index, silicon-calcium ratio, and carbon-hydrogen ratio can also indicate saturation information. The present invention can expand on uncommon parameters.

[0080] Step 4: Select the calculation algorithm corresponding to the selected saturation parameter through the user interface, automatically generate the saturation parameter results and the corresponding fitting curve plate results, and save them into the batch processing template file, which together with the simulation spectrum database serves as the saturation simulation database;

[0081] The software interface of step 3 and step 4 is as follows Figure 5 As shown, the software is divided into three areas, corresponding to the rapid loading of the simulation spectrum database in step 2, the selection of the saturation parameters to be calculated, and the selection of the calculation algorithm corresponding to the saturation parameters to be calculated in step 4, the calculation of the saturation parameters, the generation and storage of the database and plate results of the saturation parameter results.

[0082] Typical saturation parameter C / O, a calculation method (multi-detector weighting) such as the second calculation model:

[0083]

[0084] Where, CO weightThe saturation parameter result representing the saturation parameter C / O, N Cnear Represents the counts in the carbon window of the near-detector inelastic ray spectrum, N Onear Represents the counts in the oxygen window of the near-detector non-elastic ray energy spectrum, N Cfar represents the counts in the carbon window of the far detector's inelastic ray spectrum, N Ofar represents the counts in the oxygen window of the far detector's non-elastic energy spectrum, CO near Indicates the CO ratio calculated by the detector, CO far represents the CO ratio calculated by the remote detector, σ near Indicates the CO ratio near the detector under the weighted CO condition, σ far Represents the CO concentration of the far detector under weighted CO conditions. The carbon window is a window with an appropriate energy width selected near the inelastic characteristic peak of carbon, and the oxygen window is a window with an appropriate energy width selected near the inelastic characteristic peak of oxygen. ω4 and ω5 are the first and second weighting coefficients, respectively. The near and far detectors refer to the detectors closest and second closest to the source in a saturation multi-probe instrument.

[0085] It should be noted that the first weight coefficient and the second weight coefficient can be set by the user or obtained based on historical data, which is not limited in this embodiment; preferably, the ω4 is 1 and the ω5 is 1.

[0086] A typical saturation parameter Sigma is calculated using a sliding window method, such as the third calculation model. The most important thing about this model is that the counts in the next time window in the time spectrum are 1 / e of the counts in the previous time window:

[0087]

[0088] Where N1 represents the count in the previous time window, N2 represents the count in the next time window, and Δt represents the width of the time window.

[0089] Typical saturation parameters MultiINE, calculated as the fourth calculation model:

[0090]

[0091] Where N I1 It represents the total counts in all windows of the non-elastic ray spectrum of a detector in the multi-detector, N J1 Indicates the total counts in all windows of the non-elastic ray spectrum of another detector in a multi-detector.

[0092] The typical saturation parameter MultiCAP is calculated as in the fifth calculation model:

[0093]

[0094] Where N I2 Represents the total counts in all windows of the capture ray energy spectrum of a detector in the multi-probe, N J2 Represents the total counts in all windows of the capture ray energy spectrum of another detector in a multi-detector.

[0095] The typical saturation parameter FNXS is calculated as in the sixth calculation model:

[0096] FNXS i =ω6×ln(N I3 )+ω7

[0097]

[0098] Where m represents the total number of lithologies or minerals contained in the formation; FNXS i is the comprehensive FNXS of the formation under the i-th lithology or mineral; N I3 represents the total counts in all windows of the non-elastic ray spectrum of the ultra-long-range detector under the i-th lithology or mineral; ln(·) represents the natural logarithm; ω6 and ω7 represent the inversion coefficients of the comprehensive FNXS of the formation under the i-th lithology or mineral, which are obtained by calculating the relationship between the counts in the non-elastic ray spectrum window of the ultra-long-range detector under the i-th lithology or mineral and the FNXS based on the simulated spectrum database; V i It represents the volume proportion of the i-th lithology or mineral in the formation.

[0099] It should be noted that the above is a typical calculation method for saturation parameter results. The system still integrates multiple calculation methods for users to choose from. This embodiment does not make any specific limitations on the calculation method.

[0100] The fitting curve plate results include water saturation fitting curve plate results, gas saturation fitting curve plate results and oil saturation fitting curve plate results.

[0101] Step 5: Quickly load the saturation logging data, including the total spectrum of the burst period, the captured ray energy spectrum, the time spectrum and the formation porosity;

[0102] It should be noted that in the saturation well logging data, the data volume is usually large. In order to load the input information more quickly, fast loading is required. In this embodiment, fast loading uses enhanced data caching and data format conversion. The specific fast loading method is not limited in this embodiment.

[0103] Step 6: Preprocess the saturation well logging data to obtain corrected saturation well logging data;

[0104] It should be noted that the high temperature and high pressure environment in the deep sea and deep earth will affect the performance of the instrument. At the same time, due to the differences in the detectors and process errors of each actual instrument, it is impossible to guarantee that each instrument will measure the same data in the same environment. Therefore, it is necessary to pre-process the saturation measurement logging data. At the same time, under the measurement of actual instruments, only the total spectrum of the burst period can be obtained, and the inelastic ray energy spectrum cannot be directly obtained. Therefore, pre-processing of the actual well logging data is necessary.

[0105] The specific process of the pretreatment is:

[0106] Step 6.1: perform depth averaging and filtering on the total spectrum of the burst period and the captured ray energy spectrum in the saturation logging data, and then perform filtering at each depth;

[0107] It should be noted that the number of data points obtained by the actual logging instrument is related to the measurement velocity and sampling rate. Therefore, it is necessary to average and filter the saturation measured logging data at the depth level. In order to make the energy spectrum smoother, the energy spectrum needs to be filtered.

[0108] Step 6.2: Based on the total spectrum of the burst period and the energy spectrum of the captured ray obtained in step 6.1, develop and execute a fourth script to obtain the energy spectrum of the inelastic ray; the fourth script includes but is not limited to a fixed coefficient method and a double exponential fitting method;

[0109] Figure 6 The fixed coefficient method and the double exponential fitting method shown are typical applications of the fourth script, wherein a typical algorithm (fixed coefficient method) in the fourth script is as the seventh calculation model:

[0110]

[0111] Where Net represents the inelastic ray energy spectrum, Total represents the total spectrum during the burst period, Cap represents the capture ray energy spectrum, and ω8 represents the third weight coefficient. The total spectrum during the burst period refers to the energy spectrum received by the detector during the pulse time of the pulsed neutron tube.

[0112] It should be noted that the third weight coefficient may be set by the user or obtained based on historical data, and this embodiment does not limit this.

[0113] Step 6.3: Develop and execute a fifth script for correcting the inelastic ray spectrum and the capture ray spectrum to obtain a corrected inelastic gamma ray spectrum and a corrected capture gamma ray spectrum; the steps of the fifth script include peak finding, spectrum calibration, and spectrum drift correction;

[0114] It should be noted that the peak search refers to finding the characteristic peak of the element. Each element has a unique characteristic energy, and this characteristic energy is the characteristic peak. The energy spectrum scale refers to fitting the relationship between the actual channel address and energy. The spectrum drift correction is to correct the characteristic peak back to the channel address where it should be located. Figure 6 H, Si, Ca and Fe are the main elements for peak finding and calibration.

[0115] Step 7: Quickly load the saturation simulation database and the corrected saturation well data. Figure 7 The software is divided into four parts: the left side includes the input saturation simulation database results and the input saturation measured logging data;

[0116] Select the saturation parameters to be calculated and their corresponding calculation algorithms through the user interface. Figure 7 The energy spectrum mode and time spectrum mode in the saturation well data are used to generate the saturation parameter results at each well depth; combined with the corresponding fitting curve results in the saturation simulation database, the saturation parameter results of the saturation well data are scattered and projected onto the fitting curve results. Figure 7 The saturation parameters are displayed graphically in the graph. Furthermore, the saturation parameter results of the saturation well data, the fitting curve chart results, the formation porosity φ in the saturation well data, the formation lithology and the corresponding macroscopic capture cross section, and the element logging data are used to calculate the water saturation, gas saturation or oil saturation at each well depth of the saturation well data. Specifically:

[0117] Step 7.1: For C / O, determine the oil saturation S based on the formation porosity φ, saturation parameter results C / O, silicon-calcium ratio and other parameters in the saturation well data. o ;

[0118] Among them, a typical calculation method is the eighth calculation model:

[0119]

[0120] Where φ represents the formation porosity, ω9, ω 10 、ω 11 、ω 12 and ω 13 Expressed as the fourth weight coefficient, the fifth weight coefficient, the sixth weight coefficient, the seventh weight coefficient, and the eighth weight coefficient, they are determined by factors such as the saturation simulation database and the saturation measured logging data.

[0121] Step 7.2: For Sigma, calculate the macroscopic capture cross section Σ of the mudstone in the saturation measured logging data. sh , mudstone volume, oil and gas macroscopic capture cross section Σ h, macroscopic capture cross section of rock skeleton Σ ma , the macroscopic capture cross section of water in the formation pores Σ w , formation porosity φ and other parameters to jointly determine the water saturation S w ;

[0122] Among them, a typical calculation method is the ninth calculation model:

[0123]

[0124] Where, Σ represents the desired Sigma saturation parameter result;

[0125] Step 7.3: For MultiINE, MultiCAP, and FNXS, interpolate the pores according to the conditions of water and gas saturation under the simulated spectrum database to obtain the gas saturation S. g ;

[0126] Among them, a typical calculation method is the tenth calculation model:

[0127]

[0128] The saturation parameters required for calculation include FNXS, MultiINE and MultiCAP. The FNXS saturation parameter result corresponding to the gas-saturated formation in the saturation simulation database is FNXS. g The FNXS saturation parameter corresponding to the water-saturated formation is FNXS w The MultiINE saturation parameter result for the gas-saturated formation is g The MultiINE saturation parameter result for the water-saturated formation is w The MultiCAP saturation parameter result for gas-saturated formations is MultiCAP g The MultiCAP saturation parameter result corresponding to the water-saturated formation is MultiCAP w , and the FNXS saturation parameter result corresponding to the saturation logging data is recorded as FNXS meas , the MultiINE saturation parameter result is MultiINE meas , the MultiCAP saturation parameter result is MultiCAP meas ;

[0129] Calculate the FNXS gas saturation S corresponding to the saturation logging data g1 、MultiINE gas saturation S g2 and MultiCAP gas saturation S g3 ;

[0130] Then the gas saturation S is calculated g :

[0131] S g =ω1×S g1 +ω2×S g2 +ω3×S g3

[0132] ω1+ω2+ω3=1

[0133] In the formula, ω1, ω2 and ω3 are all weight coefficients, which can be set by the user or obtained based on historical data, and are not limited in this embodiment.

[0134] Example 2

[0135] This embodiment proposes a saturation multi-probe instrument nuclear calculation and processing automation system for implementing the saturation multi-probe instrument nuclear calculation and processing automation method described in Example 1. The system specifically includes a nuclear physics simulation module, a simulation spectrum database generation module, a saturation parameter and algorithm selection module, a saturation simulation database generation module, a preprocessing module, and a saturation calculation module.

[0136] The nuclear physics simulation module is used to automatically generate a nuclear well logging simulation batch model, upload it to a high-performance computing platform for simulation to generate a nuclear physics simulation output file, and output it to a simulation spectrum database generation module;

[0137] The simulation spectrum database generation module is used to extract the energy spectrum and time spectrum data in the nuclear physics simulation output file and integrate them into the batch processing template file, or convert the flux data into the gamma energy spectrum and integrate it into the batch processing template file to generate a simulation spectrum database, which is output to the saturation parameter selection module in the simulation mode;

[0138] The saturation parameter selection module is used to select the saturation parameters to be calculated and their corresponding calculation algorithms in the simulation mode or the measurement mode, and is connected to the saturation simulation database generation module;

[0139] The saturation simulation database generation module is used to automatically execute the saturation parameter calculation algorithm in the simulation mode, save the saturation parameter results calculated based on the simulation spectrum database into a batch processing template file, and use it together with the simulation spectrum database as the saturation simulation database. The calculation results corresponding to the saturation well logging data are output to the saturation calculation module.

[0140] The preprocessing module is used to preprocess the loaded saturation well logging data to obtain the corrected saturation well logging data and output it to the saturation parameter selection module in the measurement mode;

[0141] The saturation calculation module is used to calculate the water saturation, gas saturation or oil saturation at each well depth based on the saturation logging well data.

[0142] Example 3

[0143] This embodiment proposes a terminal, including a memory, a processor, and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the saturation multi-probe instrument core calculation and processing automation method described in Example 1.

[0144] Example 4

[0145] This embodiment provides a computer-readable storage medium having a computer program stored thereon. When the computer program is executed by a processor, the steps of the method for automating the calculation and processing of saturation multi-probe instrument cores described in Example 1 are implemented.

[0146] The system embodiments described above are merely illustrative. Modules described as separate components may or may not be physically separate, and components shown as modules may or may not be physical modules. They may be located in one place or distributed across multiple network modules. Some or all of these modules may be selected based on actual needs to achieve the objectives of this embodiment. Persons of ordinary skill in the art will be able to understand and implement the present invention without inventive effort.

[0147] Through the detailed description of the above embodiments, those skilled in the art can clearly understand that each embodiment can be implemented by means of software plus the necessary general hardware platform, or of course, by means of hardware. Based on this understanding, the above technical solution, in essence, or the portion that contributes to the prior art, can be embodied in the form of a software product, which can be stored in a computer-readable storage medium, including a read-only memory (ROM), a random access memory (RAM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), a one-time programmable read-only memory (OTPROM), an electronically erasable programmable read-only memory (EEPROM), a compact disc read-only memory (CD-ROM) or other optical disc storage, magnetic disk storage, magnetic tape storage, or any other computer-readable medium capable of carrying or storing data.

[0148] Finally, it should be noted that the saturation multi-probe instrument and the calculation and processing automation system and method disclosed in the embodiments of the present invention are only preferred embodiments of the present invention, and are only used to illustrate the technical solutions of the present invention, rather than to limit them. Although the present invention has been described in detail with reference to the aforementioned embodiments, it should be understood by those skilled in the art that the technical solutions described in the aforementioned embodiments can still be modified, or some of the technical features therein can be replaced by equivalents. However, these modifications or replacements do not deviate the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the various embodiments of the present invention.

Claims

1. A method for automating the calculation and processing of saturation multi-probe instrument cores, characterized in that: The following steps are involved: Step 1: Through the user interface, select the instrument parameters of the nuclear physics logging instrument to build a nuclear physics logging instrument model. Then, select the formation information, well parameters, and nuclear physics simulation parameters to automatically generate a nuclear logging simulation batch model. The model is uploaded to the high-performance computing platform for simulation calculation and generates a nuclear physics simulation output file. Step 2: Determine whether the nuclear physics simulation output file includes energy spectrum and time spectrum data. If so, directly extract the energy spectrum and time spectrum data and integrate them into the batch processing template file; otherwise, convert the flux data in the nuclear physics simulation output file into a gamma spectrum and integrate it into the batch processing template file; and then generate a simulation spectrum database in the batch processing template file; The energy spectrum includes inelastic gamma ray spectrum, capture gamma ray spectrum, thermal neutron and epithermal neutron counting spectrum; Step 3: Load the simulated spectrum database and select the saturation parameters to be calculated through the user interface, including C / O, Sigma, FNXS, MultiINE, and MultiCAP; Step 4: Select the calculation algorithm corresponding to the selected saturation parameter through the user interface, automatically generate the saturation parameter results and the corresponding fitting curve plate results, and save them into the batch processing template file, which together with the simulation spectrum database serves as the saturation simulation database; Step 5: Load the saturation logging data, including the total spectrum during the burst period, the captured ray energy spectrum, the time spectrum, the formation porosity and lithology; Step 6: Preprocess the saturation well logging data to obtain corrected saturation well logging data; Step 7: Load the saturation simulation database and the corrected saturation well data, select the saturation parameters to be calculated and their corresponding calculation algorithms through the user interface, generate saturation parameter results for the saturation well data at various well depths, and calculate the water saturation, gas saturation, or oil saturation for the saturation well data at various well depths based on the corresponding fitting curve chart results in the saturation simulation database.

2. The method for automating calculation and processing of saturation multi-probe instrument cores according to claim 1, characterized in that: The formation information described in step 1 includes lithology and pore fluid properties, the well parameters include instrument position, casing parameters, number of tubing layers, cement sheath parameters, and wellbore fluid parameters, and the nuclear physics simulation parameters include the measurement mode, particle number, and thread number of the inelastic ray energy spectrum and the capture ray energy spectrum.

3. The method for automating saturation multi-probe instrument core calculation and processing according to claim 2, characterized in that: In step 2, the detector response function matrix is ​​used to convert the flux data into a gamma spectrum, as follows: Y=XA Where, represents flux data; represents the gamma spectrum; It represents the detector response function matrix, which is calculated by simulating the detector performance of nuclear physics logging instruments using nuclear physics simulation software; n is the number of channels in the energy spectrum.

4. The method for automating calculation and processing of saturation multi-probe instrument cores according to claim 2, characterized in that: The calculation model for FNXS in step 4 is: <h2 style=";text-align:left;direction:ltr">FNXS<h2 style=";text-align:left;direction:ltr"> i <h2 style=";text-align:left;direction:ltr"> =ω6×ln(N<h2 style=";text-align:left;direction:ltr"> I3 <h2 style=";text-align:left;direction:ltr"> )+ω7 Where m represents the total number of lithologies or minerals contained in the formation; FNXS i is the comprehensive FNXS of the formation under the i-th lithology or mineral; N I3 represents the total counts in all windows of the non-elastic ray spectrum of the ultra-long-range detector under the i-th lithology or mineral; ln(·) represents the natural logarithm; ω6 and ω7 represent the inversion coefficients of the comprehensive FNXS of the formation under the i-th lithology or mineral, which are obtained by calculating the relationship between the counts in the non-elastic ray spectrum window of the ultra-long-range detector under the i-th lithology or mineral and the FNXS based on the simulated spectrum database; V i It represents the volume proportion of the i-th lithology or mineral in the formation.

5. The method for automating calculation and processing of saturation multi-probe instrument cores according to claim 2, characterized in that: The specific process of preprocessing in step 6 is: Step 6.1: perform depth averaging and filtering on the total spectrum of the burst period and the captured ray energy spectrum in the saturation logging data, and then perform filtering at each depth; Step 6.2: Based on the total spectrum of the burst period and the captured ray spectrum obtained in step 6.1, obtain the inelastic ray spectrum; Step 6.3: Correct the inelastic ray spectrum and the capture ray spectrum to obtain a corrected inelastic gamma ray spectrum and a corrected capture gamma ray spectrum.

6. The method for automating saturation multi-probe instrument core calculation and processing according to claim 2, characterized in that: The fitting curve plate results in step 4 and step 7 include the water saturation fitting curve plate results, the gas saturation fitting curve plate results, and the oil saturation fitting curve plate results.

7. The method for automating saturation multi-probe instrument core calculation and processing according to claim 6, characterized in that: The specific process of calculating the gas saturation at each well depth using the saturation logging data in step 7 is as follows: The saturation parameters required for calculation include FNXS, MultiINE and MultiCAP. The FNXS saturation parameter result corresponding to the gas-saturated formation in the saturation simulation database is FNXS. g The FNXS saturation parameter corresponding to the water-saturated formation is FNXS w The MultiINE saturation parameter result for the gas-saturated formation is g The MultiINE saturation parameter result for the water-saturated formation is w The MultiCAP saturation parameter result for gas-saturated formations is MultiCAP g The MultiCAP saturation parameter result corresponding to the water-saturated formation is MultiCAP w , and the FNXS saturation parameter result corresponding to the saturation logging data is recorded as FNXS meas , the MultiINE saturation parameter result is MultiINE meas , the MultiCAP saturation parameter result is MultiCAP meas ; Calculate the FNXS gas saturation S corresponding to the saturation logging data g1 、MultiINE gas saturation S g2 and MultiCAP gas saturation S g3 ; Then the gas saturation S is calculated g : S g =ω1×S g1 +ω2×S g2 +ω3×S g3 ω1+ω2+ω3=1 Where ω1, ω2 and ω3 are weight coefficients.

8. A saturation multi-probe instrument core calculation and processing automation system for implementing the saturation multi-probe instrument core calculation and processing automation method according to any one of claims 1 to 7, characterized in that: It includes nuclear physics simulation module, simulation spectrum database generation module, saturation parameter and algorithm selection module, saturation simulation database generation module, preprocessing module and saturation calculation module; among which: The nuclear physics simulation module is used to automatically generate a nuclear well logging simulation batch model, upload it to a high-performance computing platform for simulation to generate a nuclear physics simulation output file, and output it to a simulation spectrum database generation module; The simulation spectrum database generation module is used to extract the energy spectrum and time spectrum data in the nuclear physics simulation output file and integrate them into the batch processing template file, or convert the flux data into the gamma energy spectrum and integrate it into the batch processing template file to generate a simulation spectrum database, which is output to the saturation parameter selection module in the simulation mode; The saturation parameter selection module is used to select the saturation parameters to be calculated and their corresponding calculation algorithms in the simulation mode or the measurement mode, and is connected to the saturation simulation database generation module; The saturation simulation database generation module is used to automatically execute the saturation parameter calculation algorithm in the simulation mode, save the saturation parameter results calculated based on the simulation spectrum database into a batch processing template file, and use it together with the simulation spectrum database as the saturation simulation database. The calculation results corresponding to the saturation well logging data are output to the saturation calculation module. The preprocessing module is used to preprocess the loaded saturation well logging data to obtain the corrected saturation well logging data and output it to the saturation parameter selection module in the measurement mode; The saturation calculation module is used to calculate the water saturation, gas saturation or oil saturation at each well depth based on the saturation logging well data.

9. A terminal comprising a memory, a processor, and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of the saturation multi-probe instrument core calculation and processing automation method according to any one of claims 1 to 7.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method for automating the calculation and processing of saturation multi-probe instrument cores according to any one of claims 1 to 7 are implemented.

Images