Method for determining uranium discharge capacity of argillaceous hydrocarbon source rock and computer program product

Abstract

The application discloses a method for determining the uranium discharge amount of argillaceous hydrocarbon source rock and a computer program product, and relates to the technical field of uranium discharge amount determination. The method comprises the following steps: obtaining rock attribute data of argillaceous hydrocarbon source rock in a target layer system in a target area, wherein the rock attribute data comprises a total rock volume and a porosity average value; determining a total pore volume of the argillaceous hydrocarbon source rock according to the total rock volume and the porosity average value, and determining a total pore loss volume according to the total pore volume and the porosity average value; determining the mass of the uranium discharged from the argillaceous hydrocarbon source rock according to a corresponding geologic history period of the argillaceous hydrocarbon source rock, the total pore loss volume, and the uranium solubility of the uranium in formation water in the geologic history period, so as to quickly and accurately obtain the uranium discharge amount of the argillaceous hydrocarbon source rock in the process of geologic evolution.

Description

Technical Field

[0001] This invention relates to the field of geological exploration technology, and in particular to a method for determining the uranium content of argillaceous source rocks and a computer program product thereof. Background Technology

[0002] Uranium is both an energy resource and a strategic resource. Sandstone-type uranium deposits, due to their economic and environmental benefits, have become a key focus and hot topic in current exploration and mining. Currently, most large sandstone-type uranium deposits discovered are associated with oil and gas, and the generation and expulsion of hydrocarbons from these deposits are closely linked to uranium mineralization.

[0003] Existing methods for analyzing uranium emissions from argillaceous source rocks mainly focus on analyzing uranium compounds extracted from the source rocks using laboratory analytical techniques to estimate the amount of uranium emitted. These methods include acid leaching, flame absorption spectroscopy, fluorescence spectroscopy, and isotope analysis. However, these methods are often cumbersome and time-consuming, making it difficult to obtain evaluation results quickly and conveniently. They also have certain limitations and can easily lead to inaccurate evaluation results.

[0004] Currently, the evaluation of the resource potential of sandstone-type uranium deposits mainly focuses on uranium sources around the basin periphery, while there is no relatively complete and rapid calculation method or evaluation technology system for uranium released during the diagenesis of argillaceous source rocks. Summary of the Invention

[0005] This invention provides a method and computer program product for determining the uranium discharge of argillaceous source rocks, in order to solve the problem that the existing technology fails to fully consider the impact of pore volume loss generated by argillaceous source rocks during geological evolution on the uranium discharge, which easily leads to inaccurate evaluation results.

[0006] According to one aspect of the present invention, a method for determining the uranium discharge of argillaceous source rocks is provided, wherein the method includes:

[0007] Obtain rock property data of argillaceous source rocks in the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity.

[0008] The total pore volume of the argillaceous source rock is determined based on the total rock volume and the average porosity, and the total pore loss volume is determined based on the total pore volume and the average porosity.

[0009] The mass of uranium discharged from the argillaceous source rock is determined based on the geological history period corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period.

[0010] According to another aspect of the present invention, an apparatus for determining the uranium discharge of argillaceous source rocks is provided, wherein the apparatus comprises:

[0011] The rock property data acquisition module is used to acquire rock property data of argillaceous source rocks of the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity.

[0012] The total pore loss volume determination module is used to determine the total pore volume of the argillaceous source rock based on the total rock volume and the average porosity, and to determine the total pore loss volume based on the total pore volume and the average porosity.

[0013] The uranium discharge determination module is used to determine the mass of uranium discharged from the argillaceous source rock based on the geological history period corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period.

[0014] According to another aspect of the present invention, an electronic device is provided, the electronic device comprising:

[0015] At least one processor; and

[0016] A memory communicatively connected to the at least one processor; wherein,

[0017] The memory stores a computer program that can be executed by the at least one processor, which enables the at least one processor to perform a method for determining the uranium discharge of argillaceous source rocks according to any embodiment of the present invention.

[0018] According to another aspect of the present invention, a computer-readable storage medium is provided, the computer-readable storage medium storing computer instructions, the computer instructions being configured to cause a processor to execute and implement a method for determining the uranium discharge of a argillaceous source rock as described in any embodiment of the present invention.

[0019] According to another aspect of the present invention, embodiments of this disclosure also provide a computer program product, including a computer program that, when executed by a processor, implements a method for determining the uranium discharge of argillaceous source rocks as described in any of the embodiments of this disclosure.

[0020] The technical solution of this invention firstly involves acquiring rock property data of argillaceous source rocks in the target strata within the target area. Since the rock property data includes the total rock volume and average porosity, the physical properties of the argillaceous source rocks, such as grain size distribution, can be effectively understood, thereby reconstructing the sedimentary environment. Next, the total pore volume of the argillaceous source rocks is determined based on the total rock volume and the average porosity, and the total pore loss volume is determined based on the total pore volume and the average porosity, allowing for accurate determination of the pore volume loss of the argillaceous source rocks under compaction. Finally, based on the geological history corresponding to the argillaceous source rocks and the total pore loss... The mass of uranium discharged from the argillaceous source rock is determined by the volume and uranium solubility in formation water during the geological history period. This method comprehensively considers the limitation of fluid flow between rocks by the pore volume loss of the argillaceous source rock under compaction, thereby obtaining the volume of water discharged from the argillaceous source rock during geological history. Combined with uranium solubility, the amount of uranium discharged from the argillaceous source rock can be obtained quickly and accurately. This method enables the calculation of uranium discharge from the argillaceous source rock during geological evolution, effectively evaluating sandstone-type uranium resources in exploration blind areas, and more conveniently and quickly evaluating the resource quantity of oil and gas associated sandstone-type uranium deposits. It has good adaptability, high practical value, and strong operability.

[0021] It should be understood that the description in this section is not intended to identify key or essential features of the embodiments of the present invention, nor is it intended to limit the scope of the invention. Other features of the invention will become readily apparent from the following description. Attached Figure Description

[0022] To more clearly illustrate the technical solutions in the embodiments of the present invention, the accompanying drawings used in the description of the embodiments will be briefly introduced below. Obviously, the accompanying drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0023] Figure 1 This is a flowchart of a method for determining the uranium discharge of argillaceous source rocks according to Embodiment 1 of the present invention;

[0024] Figure 2 This is a flowchart of a method for determining the uranium discharge of argillaceous source rocks according to Embodiment 2 of the present invention;

[0025] Figure 3 This is a flowchart of a method for determining the uranium discharge of argillaceous source rocks according to Embodiment 3 of the present invention;

[0026] Figure 4 This is a schematic diagram of a device for determining the uranium discharge of argillaceous source rocks according to Embodiment 5 of the present invention;

[0027] Figure 5 This is a schematic diagram of the structure of an electronic device for implementing a method for determining the uranium discharge of argillaceous source rocks according to an embodiment of the present invention. Detailed Implementation

[0028] To enable those skilled in the art to better understand the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings of the embodiments of the present invention. Obviously, the described embodiments are only some embodiments of the present invention, and not all embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative effort should fall within the scope of protection of the present invention.

[0029] It should be noted that the terms "first," "second," "target," etc., used in the specification, claims, and accompanying drawings of this invention are used to distinguish similar objects and are not necessarily used to describe a specific order or sequence. It should be understood that such data can be interchanged where appropriate so that embodiments of the invention described herein can be implemented in orders other than those illustrated or described herein. Furthermore, the terms "comprising" and "having," and any variations thereof, are intended to cover a non-exclusive inclusion; for example, a process, method, system, product, or apparatus that comprises a series of steps or units is not necessarily limited to those steps or units explicitly listed, but may include other steps or units not explicitly listed or inherent to such processes, methods, products, or apparatus.

[0030] It should be noted that the terms "a" and "a plurality of" used in this disclosure are illustrative rather than restrictive, and those skilled in the art should understand that, unless otherwise expressly indicated in the context, they should be understood as "one or more".

[0031] The names of messages or information exchanged between multiple devices in the embodiments of this disclosure are for illustrative purposes only and are not intended to limit the scope of such messages or information.

[0032] It is understood that before using the technical solutions disclosed in the various embodiments of this disclosure, users should be informed of the types, scope of use, and usage scenarios of the personal information involved in this disclosure in an appropriate manner in accordance with relevant laws and regulations, and user authorization should be obtained.

[0033] For example, upon receiving a user's active request, a prompt message is sent to the user to explicitly inform them that the requested operation will require the acquisition and use of the user's personal information. This allows the user to independently choose whether to provide personal information to the software or hardware, such as the electronic device, application, server, or storage medium performing the operations of this disclosed technical solution, based on the prompt message.

[0034] As an optional but non-limiting implementation, in response to a user's active request, sending a prompt message to the user can be done via a pop-up window, where the prompt message can be presented in text format. Furthermore, the pop-up window can also include a selection control allowing the user to choose "agree" or "disagree" to provide personal information to the electronic device.

[0035] It is understood that the above notification and user authorization process are merely illustrative and do not constitute a limitation on the implementation of this disclosure. Other methods that comply with relevant laws and regulations may also be applied to the implementation of this disclosure.

[0036] It is understood that the data involved in this technical solution (including but not limited to the data itself, the acquisition or use of the data) shall comply with the requirements of relevant laws, regulations and related provisions.

[0037] Example 1

[0038] Figure 1 This invention provides a flowchart of a method for determining the uranium emission of argillaceous source rocks, as described in Embodiment 1. This embodiment addresses the limitations of traditional methods for analyzing uranium emission in argillaceous source rocks, which are cumbersome, time-consuming, and fail to consider the impact of pore volume loss during geological evolution on uranium emission. This method can be executed by the uranium emission determination device provided in this embodiment, which can be implemented in hardware and / or software, optionally through electronic devices such as mobile terminals, PCs, or servers. Figure 1 As shown, the method may specifically include:

[0039] S110. Obtain rock property data of argillaceous source rocks of the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity.

[0040] In this embodiment of the invention, the target area can be understood as the geological region where the argillaceous source rocks are located. The target stratigraphy typically refers to a specific stratum or series of strata within the target area that are considered to have potential oil and gas resources. Specifically, the target stratigraphy can be mudstone or shale layers rich in organic matter, possessing a certain thickness, and in a thermal evolution stage suitable for oil and gas generation. The target stratigraphy is usually determined based on multiple geological data. Selecting a specific target stratigraphy for exploration and evaluation requires comprehensive consideration of multiple geological factors, combined with geological information such as geophysical logging data, seismic data, and drilling results. The target stratigraphy may have the following characteristics: high organic carbon content, appropriate maturity, and good preservation conditions. High organic carbon content is an important indicator of source rock potential; generally, mudstone or shale with an organic carbon content exceeding 2% is considered to have good oil generation potential. Appropriate maturity means that the source rock needs to reach a certain burial depth and temperature conditions for the organic matter within it to be converted into oil and natural gas. Good preservation conditions mean that the source rock has not undergone severe post-processing alteration or oxidation, maintaining its original hydrocarbon generation capacity.

[0041] The total rock volume can be understood as the total volume of argillaceous source rocks in the target strata, specifically the sum of the rock volume including all pore spaces and solid particles. The total rock volume reflects the actual space occupied by the rock in its natural state. The average porosity of argillaceous source rocks reflects the proportion of pore space to the total rock volume.

[0042] As an optional technical solution in this embodiment of the invention, optionally, obtaining the total volume of argillaceous source rocks in the target strata within the target area includes: obtaining the distribution area and isothickness map of the argillaceous source rocks in the target strata within the target area, and determining the total volume of the argillaceous source rocks based on the distribution area and the isothickness map. Calculating the total volume of argillaceous source rocks helps to better understand the spatial distribution of source rocks and their relationship with other geological elements, thereby optimizing exploration and development strategies.

[0043] Isothepth maps, specifically source rock isothighting maps, are geological maps typically used to show the thickness variations of source rocks at a specific depth or stratum underground. An isothighting map is a contour map created from the actual thickness of strata. By connecting points with the same thickness value to form isothighting lines, it visually displays the morphology and thickness distribution characteristics of a geological body in three-dimensional space. In the oil and gas industry, isothighting maps are frequently used to determine the distribution and thickness of source rocks and reservoirs to aid exploration decisions. Isothighting maps of argillaceous source rocks allow for the visual identification of the thickness of source rock layers in different areas, effectively providing insight into the specific distribution of underground source rocks.

[0044] In some embodiments of the present invention, the thickness of the argillaceous source rock can be obtained from the isothighting map of the argillaceous source rock in the target strata within the target area. The total volume of the argillaceous source rock is obtained by multiplying the distribution area of ​​the argillaceous source rock in the target strata within the target area by the integral of the rock mass thickness. For example, the target area is divided into small units (such as a grid), and the area dS and the corresponding thickness h(x,y) of each unit are integrated to obtain the total rock volume V: V=∫∫h(x,y)dS.

[0045] In other embodiments of the present invention, the average thickness of the argillaceous source rock can be read from the isothighting map, or the average value (average thickness) can be calculated from the thickness data of multiple sampling points, and then the distribution area of ​​the argillaceous source rock is multiplied by the average thickness to obtain the total volume of the rock.

[0046] If the thickness of the argillaceous source rock varies significantly (e.g., due to faults or sedimentary differences), the region can be divided into multiple sub-regions. The volume of each sub-region can be calculated separately and then summed to obtain the total volume of the rock.

[0047] As an optional technical solution in this embodiment of the invention, optionally, obtaining the average porosity of the argillaceous source rocks of the target strata within the target area includes: obtaining the sonic transit time or density logging curves of the wells drilled within the target area, and determining the average porosity of the argillaceous source rocks of the target strata within the target area based on the sonic transit time or density logging curves. Using this technical solution, the average porosity of the argillaceous source rocks can be determined simply, quickly, and accurately.

[0048] Optionally, the porosity of argillaceous source rocks can be estimated using sonic transit time (also known as sonic velocity or sonic propagation time). Sonic transit time refers to the time required for a sound wave to travel through a rock sample and is directly related to the porosity of the argillaceous source rock. Generally, the greater the porosity, the greater the sonic transit time, because the sound wave propagates much slower in pores than in solid minerals. Optionally, the sonic transit time of the rock, as well as that of pure rock (without pores) and fluid (usually water) consisting entirely of pore space, can be obtained from sonic logging data. The porosity of the argillaceous source rock can then be calculated based on sonic transit time using empirical formulas or theoretical models (such as the Wyllie time-averaged equation). This method comprehensively considers the difference in sonic propagation velocity between the framework rock and the fluid when estimating porosity. In practical applications, other factors may need to be considered, and it may be necessary to combine this method with other logging techniques to obtain more accurate results.

[0049] Alternatively, rock porosity can also be calculated using density logging curves. Density logging tools infer porosity based on the measured bulk density of the formation. Specifically, porosity is calculated based on the density differences between rock solid particles, pore space, and the fluids filling the pores. For example, the average porosity of the argillaceous source rocks in a target area is determined based on the rock skeleton density, pore fluid density, and density logging values. The density logging value is typically the bulk density of the formation measured by a density logging instrument. The density logging value reflects the density value resulting from the combined effects of the formation rock skeleton and fluids.

[0050] During the long-term evolution of argillaceous source rocks, as they are continuously deposited and buried deep underground, volume loss typically occurs due to compaction and diagenetic changes. Specifically, as argillaceous source rocks gradually deposit from a few meters to several thousand meters deep, the pressure from the overlying strata reduces the porosity between grains, resulting in a certain degree of volume reduction. With continued deposition, the stratum thickness increases, leading not only to increased ground pressure but also to elevated geothermal temperatures. These combined effects reduce the pore volume of the source rocks, and this loss typically limits the effective expulsion of uranium, thus increasing the probability of uranium retention within the source rocks. As the porosity of argillaceous source rocks decreases due to compaction, the flow of fluids (including water, oil, and gas) is significantly restricted, directly impacting the migration paths and speeds of elements such as uranium and their compounds. Specifically, lower porosity means less space for uranium to exist in a dissolved state and migrate with fluids, thus hindering fluid flow and altering the local chemical environment, such as affecting pH and redox potential (Eh). Meanwhile, rising geothermal temperatures alter the form and solubility of uranium, affecting its retention and migration behavior within source rocks. These factors significantly impact uranium solubility and stability. Therefore, as compaction and diagenesis in argillaceous source rocks reduce porosity, it not only inhibits the migration of elements like uranium but may also promote their enrichment and stabilization under specific geological conditions, further influencing their distribution characteristics and resource potential assessment within strata.

[0051] The average porosity of argillaceous source rocks can be used to measure their storage capacity. Higher porosity usually means more space for storing organic matter and generating oil and gas, thus affecting the source rock's ability to serve as a source of oil and gas. The porosity characteristics of argillaceous source rocks help to understand groundwater flow patterns and changes in reserves, as well as potential geological evolution such as surface subsidence, making resource assessment more accurate and reducing exploration and development costs.

[0052] S120. Determine the total pore volume of the argillaceous source rock based on the total rock volume and the average porosity, and determine the total pore loss volume based on the total pore volume and the average porosity.

[0053] Specifically, determining the total pore volume of the argillaceous source rock based on the total rock volume and the average porosity includes: multiplying the total rock volume and the average porosity to obtain the total pore volume of the argillaceous source rock.

[0054] Optionally, determining the total pore loss volume based on the total pore volume and the average porosity includes: determining the total pore loss volume of the argillaceous source rock based on the total pore volume, the average porosity, and the porosity constant of the most closely packed spherulites of equal diameter.

[0055] For example, the total pore loss volume of the argillaceous source rock can be determined based on the total pore volume and the average porosity using the following formula:

[0056]

[0057] Among them, V l V represents the total volume of pore loss in argillaceous source rocks. p This represents the total pore volume of the argillaceous source rock. It represents the porosity constant of the most compactly packed spherical particles of equal diameter; This represents the average porosity of the argillaceous source rock section.

[0058] During geological history, as argillaceous source rocks are deposited and gradually buried deep underground, their porosity decreases significantly due to compaction and diagenetic changes. Using the formula described above, the difference between the original and present porosity of argillaceous source rocks from the surface to multiple depths can be quantitatively analyzed, allowing for the calculation of the total volume of porosity loss throughout the process. This enables a precise assessment of the fluid discharge volume resulting from porosity reduction, quantifying porosity loss, clarifying the specific changes in porosity at different depths, calculating the total volume of porosity loss, and estimating the accompanying fluid (such as water, oil, and gas) discharge based on porosity loss. This allows for a more accurate assessment of potential mineral resources.

[0059] S130. Determine the mass of uranium discharged from the argillaceous source rock based on the geological history period corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period.

[0060] Specifically, the mass of uranium expelled from argillaceous source rocks during geological history can be quantitatively estimated by multiplying the pore loss volume by the uranium solubility. In practical applications, corrections need to be made based on geological dynamic processes (such as compaction and hydrocarbon expulsion efficiency) and hydrochemical conditions.

[0061] This technical solution calculates the total pore loss volume of argillaceous source rocks to obtain the volume of formation water discharged by the argillaceous source rocks during a predetermined geological period. Combined with the uranium solubility in formation water during the geological period, the mass of uranium metal discharged by the argillaceous source rocks during the predetermined geological period is calculated. This overcomes the shortcomings and deficiencies of existing methods that do not consider the influence of evolutionary processes and directly evaluate based on abundance. It is extremely effective for evaluating sandstone-type uranium resources in exploration blind areas.

[0062] The technical solution of this invention firstly involves acquiring rock property data of argillaceous source rocks in the target strata within the target area. Since the rock property data includes the total rock volume and average porosity, the physical properties such as grain size distribution of the argillaceous source rocks can be effectively understood, thereby reconstructing the sedimentary environment. Next, by determining the total pore volume of the argillaceous source rocks based on the total rock volume and the average porosity, and by determining the total pore loss volume based on the total pore volume and the average porosity, the pore volume loss of the argillaceous source rocks under compaction can be accurately determined. Finally, by considering the geological history of the argillaceous source rocks and the... The total volume of pore loss and the uranium solubility in formation water during the geological history period are used to determine the mass of uranium discharged from the argillaceous source rock. This method comprehensively considers the limitation of fluid flow between rocks by the pore volume loss of the argillaceous source rock under compaction, thereby obtaining the volume of water discharged from the argillaceous source rock during the geological history period. Combined with uranium solubility, the amount of uranium discharged from the argillaceous source rock can be obtained quickly and accurately. This method realizes the calculation of the amount of uranium discharged from the argillaceous source rock during geological processes, which can effectively evaluate sandstone-type uranium resources in exploration blind areas and more conveniently and quickly evaluate the resource quantity of oil and gas associated sandstone-type uranium deposits. It has good adaptability and high practical value.

[0063] Example 2

[0064] Figure 2 This is a flowchart illustrating a method for determining the amount of uranium discharged from a argillaceous source rock, as provided in Embodiment 2 of the present invention. Based on the above embodiments, it further describes the process of determining the mass of uranium discharged from the argillaceous source rock based on the geological history period corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period. Specific implementation details can be found in the description of this embodiment. Technical features that are the same as or similar to those in the foregoing embodiments will not be repeated here. Figure 2 As shown, the method may specifically include:

[0065] S210. Obtain rock property data of argillaceous source rocks of the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity.

[0066] S220. Determine the total pore volume of the argillaceous source rock based on the total rock volume and the average porosity, and determine the total pore loss volume based on the total pore volume and the average porosity.

[0067] S230. Divide the geological history corresponding to the argillaceous source rock into multiple time periods, determine the average solubility of uranium in formation water in each time period, and determine the single-time pore loss volume in each time period based on the total pore loss volume.

[0068] In the evolution of argillaceous source rocks, a predetermined geological history period is selected. To more accurately understand and assess the behavior of uranium and its solubility changes, the geological history corresponding to the argillaceous source rocks is divided into multiple periods. The average solubility of uranium in formation water is calculated for each period, and the porosity loss volume within each period is divided based on the total porosity loss volume calculated above. Compared to a single analysis of the entire geological history, analyzing the geological history by subdividing it into multiple periods can effectively avoid the accumulation of errors caused by long time spans. The conditions (such as temperature and pressure) of each period are relatively stable, making the calculation results more accurate. Geological conditions and physicochemical environments may change significantly in different periods. Segmented analysis can more realistically reflect the changing trend of uranium solubility. Through detailed studies of each period, the complex relationship between porosity changes and uranium solubility, and how these relationships evolve over time, can be better revealed.

[0069] Based on the above scheme, optionally, determining the average solubility of uranium in formation water during each time period includes: constructing a first solubility model of uranium and a formation pressure model corresponding to the argillaceous source rock, and determining the average solubility of uranium in formation water during each time period according to the first solubility model and the formation pressure model.

[0070] To accurately assess the solubility and variation of uranium in argillaceous source rocks, a uranium solubility model is typically constructed. First, the type of uranium-bearing mineral, formation water pH, and geothermal temperature are determined, and a model is constructed to represent the changes in uranium solubility with geothermal temperature and formation water pH. Next, various laboratory methods are used to obtain uranium solubility data under different set formation temperatures and pH values. For example, batch equilibrium experiments involve mixing known amounts of uranium-bearing minerals with simulated formation water, reacting the mixture under specific temperature and pH conditions, and then measuring the uranium concentration in the solution when chemical equilibrium is reached. This leads to uranium solubility variation models corresponding to multiple formation temperatures and pH values; for example, a model based on the same temperature (25°C). At a pH of 6.0, the corresponding uranium solubility is 0.03 mg / L. At a temperature of 50°C and a pH of 7.0, the corresponding uranium solubility is 0.05 mg / L. Furthermore, during geological evolution, the pH of formation water typically exhibits relatively small changes. When constructing the uranium solubility model, it can be assumed that the pH value is a constant or fluctuates within a narrow range. This simplifies the model's complexity, making it easier to process and apply. Optionally, the first solubility model is used to indicate the correspondence between formation temperature and uranium solubility.

[0071] In deep geological environments, formation pressure and temperature often increase proportionally with depth. In such environments, the solubility of uranium-bearing minerals is not influenced by a single factor but is the result of the interaction of multiple factors. Therefore, it is necessary to construct a model that comprehensively considers the effects of formation pressure and temperature on uranium solubility. This model will help to more accurately predict the behavior and solubility of uranium at different geological depths. Optionally, the formation pressure model is used to indicate the correspondence between formation temperature and formation pressure at multiple time points throughout geological history.

[0072] Optionally, a formation pressure model can be constructed using basin simulation software (PetroMod, BaseinMod) combined with fluid inclusion analysis. Fluid inclusions are small volumes of fluid or gas trapped within minerals during crystallization; they preserve information about the physicochemical conditions of the environment at that time, including temperature, pressure, and fluid composition. Analyzing these inclusions can yield crucial data on the properties of fluids under paleogeological conditions.

[0073] Optionally, a formation pressure model can be constructed using basin simulation software combined with fluid inclusion analysis. This includes: collecting basic geological data of the target basin, including but not limited to drilling data, seismic profiles, and stratigraphic columns; collecting rock samples containing fluid inclusions from the target area, typically selecting minerals formed during key geological events; observing the morphology and distribution of inclusions using a microscope and determining their homogenization temperature and freezing point using cryo-thermal stage technology to infer the pressure and temperature conditions at the time of capture; analyzing the chemical composition of the fluids within the inclusions using methods such as laser Raman spectroscopy and secondary ion mass spectrometry to further understand the geological environment at the time, supplementing the data obtained; creating a three-dimensional structural model of the basin using PetroMod or BasicMod software, defining the properties of each stratigraphic unit (such as thickness, lithology, porosity, etc.), and setting the thermal history model based on known geothermal gradients, thermal conductivity, and other parameters; using actual pressure and temperature data obtained from fluid inclusion analysis to calibrate the model, ensuring that the simulation results accurately reflect the actual situation. By comparing the simulation results with actual measurements (such as downhole pressure tests), necessary adjustments are made to the model to improve its accuracy. Run basin simulation software to simulate the evolution of the basin throughout geological history, thereby understanding the impact of sedimentation, compaction, diagenesis, and tectonic activity on formation pressure and temperature during geological evolution.

[0074] Optionally, during the construction of the formation pressure model, the uplift period portion of the geological history can be removed and denoted as P, representing the pore loss volume V of the argillaceous source rock. l The corresponding changes in formation pressure are analyzed. During the uplift phase, due to the uplift of the lithosphere, rocks that were originally in a deep, high-temperature, and high-pressure environment are brought to shallower or near-surface locations. This leads to changes in the geothermal gradient, typically manifested as a faster rate of geothermal rise, deviating from the normal thermal evolution path. Simultaneously, the uplift process is accompanied by a reduction in overlying strata, resulting in a decrease in formation pressure. This significantly impacts studies of pore pressure and fluid activity, potentially masking formation pressure characteristics formed during normal sedimentation. In such cases, radiometric dating can be used to determine the timeframe of the uplift event. By excluding data from the uplift phase and analyzing only the relatively stable subsidence stage, the simulation results can accurately reflect the pressure and temperature conditions during the basin's subsidence phase. This allows for a more accurate capture of the basin's natural evolution over long timescales, avoiding misleading results caused by abnormal conditions during the uplift phase.

[0075] Based on the above scheme, optionally, determining the single-period pore loss volume within each time period based on the total pore loss volume includes: determining the average formation pressure corresponding to each time period according to the formation pressure model, and determining the pressure ratio of the average formation pressure to the total formation pressure within the geological history period; determining the single-period pore loss volume within each time period based on the total pore loss volume and the pressure ratio.

[0076] Optionally, the pore loss volume V of the argillaceous source rock can be obtained based on the formation pressure model. l The corresponding formation pressure variation P is divided into n equal parts according to the geological history period T corresponding to P, with each part having a geological time length of Δt, in Ma (millions of years). Among them, Δt=Δt1=Δt2=…=Δt n Where n is the sequence number corresponding to Δt, which is a dimensionless quantity; the average formation pressure corresponding to each equal geological time segment is... Recorded as The unit is MPa (megapascal). Calculate the values ​​of each... The proportional relationship between them can be represented by r, where r is a dimensionless vector; based on the pressure ratio r between the average formation pressure and the total formation pressure over geological history, the pore loss volume V of argillaceous source rocks is calculated. l Constructing a unit geological time length Δt m The corresponding pore volume loss of the argillaceous source rock, i.e. and in The unit is m 3 .

[0077] In this technical solution, the total pore loss volume of the argillaceous source rock is known. The pore loss volume is directly proportional to the change in formation pressure. The pressure ratio of the average formation pressure corresponding to each time period to the total formation pressure throughout the entire geological history is multiplied by the total pore loss volume to obtain the single-time pore loss volume within each time period. By dividing the pore loss volume into different geological stages, a more detailed understanding of the porosity changes of the argillaceous source rock in different geological evolution stages and its mapping relationship with formation pressure can be obtained. By calculating the pore loss volume in segments, the specific changes in porosity in different geological stages can be captured more accurately, avoiding the accumulation of errors caused by long time spans and improving accuracy.

[0078] S240. Determine the mass of uranium discharged from the argillaceous source rock based on the single-period pore loss volume of the multiple time periods and the average solubility corresponding to the time period.

[0079] Specifically, the uranium mass can be obtained by multiplying the single-period pore loss volume and the average solubility within each time period and then summing them; the uranium mass can be multiplied by the conversion constant coefficient between uranium mass and uranium mass to obtain the mass of uranium discharged from the argillaceous source rock.

[0080] For example, the mass of metallic uranium discharged from the argillaceous source rock can be determined based on the following formula, which considers the single-period pore loss volume and the average solubility corresponding to the multiple said time periods. U (Unit: 10) 3 kg):

[0081]

[0082] Among them, M U This indicates the mass of metallic uranium extruded from argillaceous source rocks, in units of 10⁻⁶. 3 kg; Δt represents the set unit geological time length, unit: Ma; n represents the total number of Δt within the total geological time length T, which is a dimensionless quantity; m represents the sequence number corresponding to Δt, which is a dimensionless quantity; This represents the pore loss volume of the argillaceous source rock within a unit geological time length Δt, with sequence number m (unit: m). 3 ); This represents the average solubility of uranium (unit: mol / m) within a unit geological time length Δt of sequence number m. 3 ); a represents the conversion constant between the mass of uranium and the mass of uranium.

[0083] This technical solution calculates the pore loss volume of argillaceous source rocks for each time period within a preset geological history, thereby obtaining the volume of formation water discharged due to the pore loss of argillaceous source rocks in each time period. It then calculates the uranium solubility for each time period, thus obtaining the mass of metallic uranium discharged from the formation water in each time period. By constructing a mapping relationship between the pore volume loss of argillaceous source rocks and the amount of uranium discharged, the amount of uranium discharged from argillaceous source rocks can be obtained quickly and accurately, demonstrating strong practicality.

[0084] The technical solution of this application, firstly, divides the geological history of the argillaceous source rock into multiple time periods, determines the average solubility of uranium in formation water within each time period, and determines the single-time pore loss volume within each time period based on the total pore loss volume. This allows for a comprehensive and in-depth understanding of the porosity and uranium solubility changes of the argillaceous source rock in different geological historical stages, providing data support for subsequently constructing a mapping relationship between porosity loss and uranium discharge in the argillaceous source rock. Next, based on the single-time pore loss volume and the average solubility corresponding to each time period, the mass of uranium discharged from the argillaceous source rock is determined. Due to changes in environmental conditions (such as temperature and pH), the solubility of uranium may vary in different time periods. By measuring and calculating in different time periods, these changes can be captured more accurately, thereby achieving dynamic monitoring of uranium discharge and obtaining the uranium discharge from the argillaceous source rock throughout the entire geological history. This provides a more convenient and rapid assessment of the resource quantity of oil and gas associated sandstone-type uranium deposits, exhibiting good adaptability, high practical value, and strong operability.

[0085] Example 3

[0086] Figure 3 This is a flowchart illustrating a method for determining uranium discharge from argillaceous source rocks according to Embodiment 3 of the present invention. Based on the above embodiments, it further describes the process of determining the average solubility of uranium in formation water during each time period according to the first solubility model and the formation pressure model. Specific implementation details can be found in the description of this embodiment. Technical features that are the same as or similar to those in the foregoing embodiments will not be repeated here. Figure 3 As shown, the method in this embodiment may specifically include:

[0087] S310. Obtain rock property data of argillaceous source rocks of the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity.

[0088] S320. Determine the total pore volume of the argillaceous source rock based on the total rock volume and the average porosity, and determine the total pore loss volume based on the total pore volume and the average porosity.

[0089] S330. Construct a first solubility model for uranium and a formation pressure model corresponding to the argillaceous source rock. Construct a second solubility model based on the first solubility model and the formation pressure model. The second solubility model is used to indicate the correspondence between the average formation pressure and the solubility of uranium.

[0090] Specifically, the first solubility model is used to obtain the correspondence between formation temperature and uranium solubility, and the formation pressure model is used to obtain the mapping relationship between formation pressure and formation temperature. The second solubility model is obtained based on the formation temperature and uranium solubility corresponding to the same formation pressure. The second solubility model is used to indicate the correspondence between average formation pressure and uranium solubility. By using formation temperature as an intermediate variable, the correspondence between average formation pressure and uranium solubility is established. Based on the established second solubility model, the changes in uranium solubility under different formation pressure conditions can be understood more accurately.

[0091] S340. Divide the geological history period corresponding to the argillaceous source rock into multiple time periods. For each time period, determine the average solubility of uranium in formation water within that time period based on the average formation pressure corresponding to that time period and the second solubility model.

[0092] Optionally, for each period in geological history, based on the length of each equal geological time segment of the geological history period and the correspondence between average formation pressure and uranium solubility, the average solubility of uranium in formation water corresponding to the average formation pressure of argillaceous source rocks in each period can be obtained, and the unit geological time length Δt can be converted into the average solubility of uranium in formation water during that period. m The average solubility of uranium within is expressed as Right now The unit is mol / m 3 By dividing geological time periods into equal segments based on formation pressure, the average solubility of uranium in formation water within each period can be obtained. By comprehensively considering the changes in average formation pressure and uranium solubility of argillaceous source rocks in each period, the solubility changes of uranium at each stage of geological history can be reconstructed in greater detail. Through the fine division of different geological periods, the changing patterns of uranium solubility over time can be captured more accurately.

[0093] S350. Determine the single-period pore loss volume within each time period based on the total pore loss volume.

[0094] S360. Determine the mass of uranium discharged from the argillaceous source rock based on the single-period pore loss volume of the multiple time periods and the average solubility corresponding to the time period.

[0095] The technical solution of this invention can reconstruct the dynamics of geological processes and capture the pressure-solubility response characteristics at different stages by constructing a first solubility model for uranium and a formation pressure model corresponding to the argillaceous source rock. By constructing a second solubility model based on the first and formation pressure models, and since the second solubility model indicates the correspondence between average formation pressure and uranium solubility, the specific path of pressure's effect on solubility can be quantitatively analyzed. The geological history corresponding to the argillaceous source rock is divided into multiple time periods. For each time period, the average solubility of uranium in formation water is determined based on the average formation pressure and the second solubility model. The single-time pore loss volume within each time period is determined based on the total pore loss volume. This decomposes the macroscopic geological process into quantifiable stages, achieving a refined analysis of the coupling relationship between uranium solubility, pressure, and pore loss. This method reveals the temporal patterns of uranium migration and significantly improves the accuracy and reliability of uranium discharge.

[0096] Example 4

[0097] As an optional example of the present invention, the method for determining the uranium discharge from argillaceous source rocks may specifically include:

[0098] (1) Determine the target area;

[0099] (2) Obtain the distribution area and isothickness map of the target strata within the target area, and determine the total volume of the rock based on the distribution area and isothickness map;

[0100] (3) Determine the current average porosity of the argillaceous source rock section based on the sonic transit time or density logging curves of the wells drilled in the target area. Porosity is expressed as a percentage (%).

[0101] (4) Determine the total pore volume V of the argillaceous source rock based on the total rock volume and the average porosity. p (Unit: m) 3 );

[0102] (5) Determine the total pore loss volume V based on the total pore volume and the average porosity using the following formula (Formula 1). l (Unit: m) 3 ):

[0103]

[0104] in:

[0105] V l Total pore loss volume of argillaceous source rocks, in cubic meters (m³). 3 ;

[0106] V p The total pore volume of present-day argillaceous source rocks, in cubic meters (m³). 3 ;

[0107] The porosity constant for the closest packing of equal-diameter spherical particles is taken as 25.95%.

[0108] The average porosity of the current argillaceous source rock section, in %;

[0109] (6) Construct a first solubility model for uranium solubility as a function of formation temperature and a formation pressure model for formation pressure as a function of formation temperature during geological history.

[0110] (7) Based on the characteristics of formation pressure history, the part in the formation pressure model that is in the uplift period is removed, and it can be set as P, which is the pore loss volume V of the argillaceous source rock. l The corresponding changes in formation pressure;

[0111] (8) Divide the geological history period T corresponding to P into n equal parts, with each part having a geological time length of Δt, in Ma (millions of years). Among them, Δt=Δt1=Δt2=…=Δt n Where n is the sequence number corresponding to Δt, which is a dimensionless quantity; the average formation pressure corresponding to each equal geological time segment is... Recorded as The unit is MPa (megapascal). Calculate the values ​​of each... and The proportional relationship between the sums can be represented by r, where r is a dimensionless vector.

[0112] (9) Based on r, the pore loss volume V of the argillaceous source rock l The pore volume loss of argillaceous source rocks corresponding to a unit geological time length Δt is constructed, i.e. and in The unit is m 3 ;

[0113] (10) Using formation temperature as a medium, a second solubility model is constructed based on the first solubility model and the formation pressure model to indicate the correspondence between mean formation pressure and uranium solubility, and the solubility can be expressed as a unit geological time length Δt. m The average solubility of uranium within is expressed as Right now The unit is mol / m 3 ;

[0114] (11) Calculate the mass of metallic uranium discharged from the argillaceous source rock M using the following formula (Equation 2).U :

[0115]

[0116] in:

[0117] M U Mass of metallic uranium extruded from argillaceous source rocks, unit: 10⁻⁶ 3 kg;

[0118] △t: The set unit geological time length, in Ma;

[0119] n: The number of Δt within the total geological time length T, which is a dimensionless quantity;

[0120] m: The sequence number corresponding to △t, which is a dimensionless quantity;

[0121] The pore loss volume of argillaceous source rocks within a unit geological time length Δt for sequence number m (unit: m³) 3 );

[0122] Average solubility of uranium (unit: mol / m) within a unit geological time length Δt for sequence number m. 3 );

[0123] 238.08291×10 -6 : The conversion constant between uranium mass and uranium mass.

[0124] The technical solution of this invention can accurately depict the environmental conditions of each period in geological history, and specifically analyze the mapping relationship between formation pressure and uranium solubility in each period. This can provide a new perspective for oil and gas exploration, help determine the optimal drilling location, optimize mining strategies, and contribute to breakthroughs in sandstone-type uranium deposit exploration. Based on the pore volume loss of argillaceous source rocks, and using uranium solubility and formation pressure models, the amount of uranium excreted from argillaceous source rocks during geological processes can be calculated. This is a supplement to the resource evaluation technology of oil and gas associated sandstone-type uranium deposits, overcoming the defects and shortcomings of existing methods that directly evaluate based on abundance. It is extremely effective for evaluating sandstone-type uranium resources in exploration blind areas. This method takes a new perspective on the resource evaluation of oil and gas associated sandstone-type uranium deposits. It has a simple structure, reasonable design, and can more conveniently and quickly evaluate the resource quantity of oil and gas associated sandstone-type uranium deposits. It has good adaptability, high practical value, and strong operability.

[0125] Example 5

[0126] Figure 4This is a schematic diagram of a device for determining the uranium emission amount of argillaceous source rocks according to Embodiment 5 of the present invention. This device is used to execute the method for determining the uranium emission amount of argillaceous source rocks provided in any of the above embodiments. This device and the methods for determining the uranium emission amount of argillaceous source rocks in the above embodiments belong to the same inventive concept. Details not described in detail in the embodiments of the device for determining the uranium emission amount of argillaceous source rocks can be referred to the embodiments of the methods for determining the uranium emission amount of argillaceous source rocks described above. Figure 4 As shown, the device includes: a rock property data acquisition module 410, a total pore loss volume determination module 420, and a uranium discharge quantity determination module 430.

[0127] The rock attribute data acquisition module 410 is used to acquire rock attribute data of argillaceous source rocks in the target strata within the target area, wherein the rock attribute data includes the total volume of the rock and the average porosity; the total pore loss volume determination module 420 is used to determine the total pore volume of the argillaceous source rock based on the total volume of the rock and the average porosity, and to determine the total pore loss volume based on the total pore volume and the average porosity; the uranium discharge volume determination module 430 is used to determine the mass of uranium discharged from the argillaceous source rock based on the geological history period corresponding to the argillaceous source rock, the total pore loss volume, and the uranium solubility in formation water during the geological history period.

[0128] The technical solution of this invention embodiment firstly acquires rock attribute data of argillaceous source rocks in the target strata within the target area through a rock attribute data acquisition module 410. Since the rock attribute data includes the total rock volume and the average porosity, it can effectively understand the physical properties of the argillaceous source rocks, such as grain size distribution, thereby reconstructing the sedimentary environment. Next, a total pore loss volume determination module 420 determines the total pore volume of the argillaceous source rocks based on the total rock volume and the average porosity, and determines the total pore loss volume based on the total pore volume and the average porosity, which can accurately determine the pore volume loss of the argillaceous source rocks under compaction. Finally, a uranium discharge volume determination module 430 determines the total pore volume of the argillaceous source rocks based on the total pore volume and the average porosity, thereby accurately determining the pore volume loss of the argillaceous source rocks under compaction. The geological history corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period determine the mass of uranium discharged from the argillaceous source rock. This method comprehensively considers the limitation of fluid flow between rocks by the pore volume loss of the argillaceous source rock under compaction, thereby obtaining the volume of water discharged from the argillaceous source rock during the geological history period. Combined with uranium solubility, the amount of uranium discharged from the argillaceous source rock can be obtained quickly and accurately. This method realizes the calculation of the amount of uranium discharged from the argillaceous source rock during geological processes, which can effectively evaluate sandstone-type uranium resources in exploration blind areas and more conveniently and quickly evaluate the resource quantity of oil and gas associated sandstone-type uranium deposits. It has good adaptability and high practical value.

[0129] Based on the above scheme, optionally, the rock property data acquisition module 410 includes a rock total volume determination submodule. This rock total volume determination submodule is used to acquire the distribution area and isothickness map of the argillaceous source rocks of the target strata within the target area, and to determine the total rock volume of the argillaceous source rocks based on the distribution area and the isothickness map.

[0130] Based on the above scheme, optionally, the rock property data acquisition module 410 includes a porosity average value determination submodule. This porosity average value determination submodule is used to acquire the sonic transit time or density logging curves of the wells drilled within the target area, and determine the average porosity of the argillaceous source rocks of the target strata within the target area based on the sonic transit time or density logging curves.

[0131] Based on the above scheme, optionally, the total pore loss volume determination module is specifically used to determine the total pore loss volume based on the total pore volume and the average porosity using the following formula:

[0132]

[0133] Among them, V l V represents the total volume of pore loss in argillaceous source rocks. p This represents the total pore volume of the argillaceous source rock. It represents the porosity constant of the most compactly packed spherical particles of equal diameter; This represents the average porosity of the argillaceous source rock section.

[0134] Based on the above scheme, optionally, the uranium discharge determination module includes: a single-period pore loss volume determination submodule and a uranium discharge volume determination submodule. The single-period pore loss volume determination submodule is used to divide the geological history corresponding to the argillaceous source rock into multiple periods, determine the average solubility of uranium in formation water within each period, and determine the single-period pore loss volume within each period based on the total pore loss volume; the uranium discharge volume determination submodule is used to determine the mass of uranium discharged from the argillaceous source rock based on the single-period pore loss volumes of the multiple periods and the average solubility corresponding to the periods.

[0135] Based on the above scheme, optionally, the single-period pore loss volume determination submodule includes: a unit for determining the average uranium solubility in each period. This unit is used to construct a first uranium solubility model and a formation pressure model corresponding to the argillaceous source rock, and to determine the average solubility of uranium in formation water for each period based on the first solubility model and the formation pressure model. The first solubility model indicates the correspondence between formation temperature and uranium solubility, and the formation pressure model indicates the correspondence between formation temperature and formation pressure at multiple time points within geological history.

[0136] Based on the above scheme, optionally, the single-period pore loss volume determination submodule includes: a pressure ratio determination unit and a single-period pore loss volume determination unit. The pressure ratio determination unit is used to determine the average formation pressure corresponding to each of the time periods according to the formation pressure model, and to determine the pressure ratio between the average formation pressure and the total formation pressure during the geological history period; the single-period pore loss volume determination unit is used to determine the single-period pore loss volume within each of the time periods based on the total pore loss volume and the pressure ratio.

[0137] Based on the above scheme, optionally, the uranium discharge determination submodule includes: a uranium mass acquisition unit and a uranium discharge quantity acquisition unit. The uranium mass acquisition unit is used to multiply the single-period pore loss volume and the average solubility within each time period and then sum them to obtain the uranium mass; the uranium discharge quantity acquisition unit is used to multiply the uranium mass by a conversion constant coefficient between uranium mass and uranium mass to obtain the mass of uranium discharged from the argillaceous source rock.

[0138] Based on the above scheme, optionally, the unit for determining the average solubility of uranium in each time period is specifically used to: construct a second solubility model based on the first solubility model and the formation pressure model; wherein the second solubility model is used to indicate the correspondence between the average formation pressure and the uranium solubility; and, for each time period, determine the average solubility of uranium in formation water within the time period based on the average formation pressure corresponding to the time period and the second solubility model.

[0139] The uranium emission determination device for argillaceous source rocks provided in Embodiment 4 of the present invention can execute the uranium emission determination method for argillaceous source rocks provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the method.

[0140] Example 6

[0141] Figure 5A schematic diagram of an electronic device 10 that can be used to implement embodiments of the present invention is shown. The electronic device is intended to represent various forms of digital computers, such as laptop computers, desktop computers, workstations, personal digital assistants, servers, blade servers, mainframe computers, and other suitable computers. The electronic device can also represent various forms of mobile devices, such as personal digital processors, cellular phones, smartphones, wearable devices (e.g., helmets, glasses, watches, etc.), and other similar computing devices. The components shown herein, their connections and relationships, and their functions are merely illustrative and are not intended to limit the implementation of the invention described and / or claimed herein.

[0142] like Figure 5 As shown, the electronic device 10 includes at least one processor 11 and a memory, such as a read-only memory (ROM) 12 or a random access memory (RAM) 13, communicatively connected to the at least one processor 11. The memory stores computer programs executable by the at least one processor. The processor 11 can perform various appropriate actions and processes based on the computer program stored in the ROM 12 or loaded from storage unit 18 into the RAM 13. The RAM 13 may also store various programs and data required for the operation of the electronic device 10. The processor 11, ROM 12, and RAM 13 are interconnected via a bus 14. An input / output (I / O) interface 15 is also connected to the bus 14.

[0143] Multiple components in electronic device 10 are connected to I / O interface 15, including: input unit 16, such as keyboard, mouse, etc.; output unit 17, such as various types of displays, speakers, etc.; storage unit 18, such as disk, optical disk, etc.; and communication unit 19, such as network card, modem, wireless transceiver, etc. Communication unit 19 allows electronic device 10 to exchange information / data with other devices through computer networks such as the Internet and / or various telecommunications networks.

[0144] Processor 11 can be a variety of general-purpose and / or special-purpose processing components with processing and computing capabilities. Some examples of processor 11 include, but are not limited to, a central processing unit (CPU), a graphics processing unit (GPU), various special-purpose artificial intelligence (AI) computing chips, various processors running machine learning model algorithms, a digital signal processor (DSP), and any suitable processor, controller, microcontroller, etc. Processor 11 performs the various methods and processes described above, such as a method for determining the uranium discharge of argillaceous source rocks.

[0145] In some embodiments, a method for determining the uranium emissions of a argillaceous source rock can be implemented as a computer program tangibly contained in a computer-readable storage medium, such as storage unit 18. In some embodiments, part or all of the computer program can be loaded and / or installed on electronic device 10 via ROM 12 and / or communication unit 19. When the computer program is loaded into RAM 13 and executed by processor 11, one or more steps of the method for determining the uranium emissions of a argillaceous source rock described above can be performed. Alternatively, in other embodiments, processor 11 can be configured to perform a method for determining the uranium emissions of a argillaceous source rock by any other suitable means (e.g., by means of firmware).

[0146] Various embodiments of the systems and techniques described above herein can be implemented in digital electronic circuit systems, integrated circuit systems, field-programmable gate arrays (FPGAs), application-specific integrated circuits (ASICs), application-specific standard products (ASSPs), systems-on-a-chip (SoCs), payload-programmable logic devices (CPLDs), computer hardware, firmware, software, and / or combinations thereof. These various embodiments may include implementations in one or more computer programs that can be executed and / or interpreted on a programmable system including at least one programmable processor, which may be a dedicated or general-purpose programmable processor, capable of receiving data and instructions from a storage system, at least one input device, and at least one output device, and transmitting data and instructions to the storage system, the at least one input device, and the at least one output device.

[0147] Computer programs used to implement the methods of the present invention may be written in any combination of one or more programming languages. These computer programs may be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, such that when executed by the processor, the computer programs cause the functions / operations specified in the flowcharts and / or block diagrams to be performed. The computer programs may be executed entirely on a machine, partially on a machine, or as a standalone software package, partially on a machine and partially on a remote machine, or entirely on a remote machine or server.

[0148] In the context of this invention, a computer-readable storage medium can be a tangible medium that may contain or store a computer program for use by or in conjunction with an instruction execution system, apparatus, or device. A computer-readable storage medium may include, but is not limited to, electronic, magnetic, optical, electromagnetic, infrared, or semiconductor systems, apparatus, or devices, or any suitable combination thereof. Alternatively, a computer-readable storage medium may be a machine-readable signal medium. More specific examples of machine-readable storage media include electrical connections based on one or more wires, portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), optical fibers, portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof.

[0149] To provide interaction with a user, the systems and techniques described herein can be implemented on an electronic device having: a display device (e.g., a CRT (cathode ray tube) or LCD (liquid crystal display) monitor) for displaying information to the user; and a keyboard and pointing device (e.g., a mouse or trackball) through which the user provides input to the electronic device. Other types of devices can also be used to provide interaction with the user; for example, feedback provided to the user can be any form of sensory feedback (e.g., visual feedback, auditory feedback, or tactile feedback); and input from the user can be received in any form (including sound input, voice input, or tactile input).

[0150] The systems and technologies described herein can be implemented in computing systems that include backend components (e.g., as data servers), or computing systems that include middleware components (e.g., application servers), or computing systems that include frontend components (e.g., user computers with graphical user interfaces or web browsers through which users can interact with implementations of the systems and technologies described herein), or any combination of such backend, middleware, or frontend components. The components of the system can be interconnected via digital data communication of any form or medium (e.g., communication networks). Examples of communication networks include local area networks (LANs), wide area networks (WANs), blockchain networks, and the Internet.

[0151] A computing system can include clients and servers. Clients and servers are generally located far apart and typically interact through communication networks. The client-server relationship is created by computer programs running on the respective computers and having a client-server relationship with each other. The server can be a cloud server, also known as a cloud computing server or cloud host, which is a hosting product within the cloud computing service system to address the shortcomings of traditional physical hosts and VPS services, such as high management difficulty and weak business scalability.

[0152] In particular, according to embodiments of the present invention, the processes described above with reference to the flowcharts can be implemented as computer software programs. For example, embodiments of the present invention include a computer program product comprising a computer program carried on a non-transitory computer-readable medium, the computer program containing program code for performing the methods shown in the flowcharts. In such embodiments, the computer program can be downloaded and installed from a network via communication unit 19, or installed from storage unit 18, or installed from ROM 12. When the computer program is executed by processor 11, it performs the functions defined in the methods of the embodiments of the present invention.

[0153] It should be understood that the various forms of processes shown above can be used, with steps reordered, added, or deleted. For example, the steps described in this invention can be executed in parallel, sequentially, or in different orders, as long as the desired result of the technical solution of this invention can be achieved, and this is not limited herein.

[0154] The specific embodiments described above do not constitute a limitation on the scope of protection of this invention. Those skilled in the art should understand that various modifications, combinations, sub-combinations, and substitutions can be made according to design requirements and other factors. Any modifications, equivalent substitutions, and improvements made within the spirit and principles of this invention should be included within the scope of protection of this invention.

Claims

1. A method for determining the uranium discharge from argillaceous source rocks, characterized in that, include: Obtain rock property data of argillaceous source rocks in the target strata within the target area, wherein the rock property data includes the total rock volume and the average porosity. The total pore volume of the argillaceous source rock is determined based on the total rock volume and the average porosity, and the total pore loss volume is determined based on the total pore volume and the average porosity. The mass of uranium discharged from the argillaceous source rock is determined based on the geological history period corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history period.

2. The method according to claim 1, characterized in that, The determination of the mass of uranium discharged from the argillaceous source rock based on the geological history corresponding to the argillaceous source rock, the total volume of pore loss, and the uranium solubility in formation water during the geological history includes: The geological history corresponding to the argillaceous source rock is divided into multiple time periods. The average solubility of uranium in formation water in each time period is determined, and the single-time pore loss volume in each time period is determined based on the total pore loss volume. The mass of uranium discharged from the argillaceous source rock is determined based on the single-period pore loss volume and the average solubility corresponding to the period.

3. The method according to claim 2, characterized in that, The determination of the average solubility of uranium in formation water during each of the aforementioned time periods includes: A first solubility model for uranium and a formation pressure model corresponding to the argillaceous source rock are constructed. The average solubility of uranium in formation water during each time period is determined based on the first solubility model and the formation pressure model. The first solubility model is used to indicate the correspondence between formation temperature and uranium solubility, and the formation pressure model is used to indicate the correspondence between formation temperature and formation pressure at multiple time points in geological history.

4. The method according to claim 3, characterized in that, The step of determining the average solubility of uranium in formation water for each time period based on the first solubility model and the formation pressure model includes: A second solubility model is constructed based on the first solubility model and the formation pressure model; wherein, the second solubility model is used to indicate the correspondence between the average formation pressure and the uranium solubility; For each time period, the average solubility of uranium in formation water during that time period is determined based on the average formation pressure corresponding to that time period and the second solubility model.

5. The method according to claim 3, characterized in that, The step of determining the single-period pore loss volume within each time period based on the total pore loss volume includes: The average formation pressure for each time period is determined according to the formation pressure model, and the pressure ratio of the average formation pressure to the total formation pressure during the geological history period is determined. The single-period pore loss volume within each time period is determined based on the total pore loss volume and the pressure ratio.

6. The method according to claim 2, characterized in that, The determination of the uranium discharge mass from the argillaceous source rock based on the single-period pore loss volume and the average solubility corresponding to the multiple time periods includes: The mass of uranium is obtained by multiplying the single-period pore loss volume and the average solubility within each time period and then summing them. Multiplying the mass of the uranium by the conversion constant between the mass of the uranium and the mass of uranium yields the mass of the uranium discharged from the argillaceous source rock.

7. The method according to claim 1, characterized in that, The total pore loss volume is determined based on the total pore volume and the average porosity using the following formula: Among them, V l V represents the total volume of pore loss in argillaceous source rocks. p This represents the total pore volume of the argillaceous source rock. It represents the porosity constant of the most compactly packed spherical particles of equal diameter; This represents the average porosity of the argillaceous source rock section.

8. The method according to claim 1, characterized in that, The acquisition of the total rock volume of the argillaceous source rocks of the target strata within the target area includes: Obtain the distribution area and isothickness map of the argillaceous source rocks of the target strata within the target area, and determine the total volume of the argillaceous source rocks based on the distribution area and the isothickness map.

9. The method according to claim 1, characterized in that, The process of obtaining the average porosity of the argillaceous source rocks in the target strata within the target area includes: Obtain the sonic transit time or density logging curves of the wells drilled in the target area, and determine the average porosity of the argillaceous source rocks of the target strata in the target area based on the sonic transit time or density logging curves.

10. A computer program product, comprising a computer program, characterized in that, When executed by a processor, the computer program implements the method for determining the uranium discharge of argillaceous source rocks as described in any one of claims 1-9.

Images