A method, device, equipment and medium for obtaining multidimensional nuclear magnetic resonance spectrum
By constructing a three-dimensional digital core model and numerical simulation method, the problem of low acquisition efficiency of multidimensional nuclear magnetic resonance spectra in the existing technology is solved, and the multidimensional nuclear magnetic resonance spectra is obtained efficiently and at low cost is achieved, and the response mechanism in the CO2 storage process in rocks is studied.
Patent Information
- Application Number
- CN202310166057.8
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2023-02-22
- Publication Date
- 2025-08-22
- Estimated Expiration
- 2043-02-22
AI Technical Summary
The prior art is difficult to obtain multidimensional nuclear magnetic resonance spectra efficiently and at low cost, especially in the process of CO2 geological storage in rocks, the response mechanism of saturated oil, water and CO2 is difficult to determine.
By constructing a three-dimensional digital core model, the process of CO2 injection into saturated oil and water rocks was simulated by morphological methods, and the multidimensional nuclear magnetic resonance echo data was obtained using numerical simulation methods, and finally inversion was performed to obtain the multidimensional nuclear magnetic resonance spectrum.
The complexity and cost of multi-dimensional NMR spectrum acquisition is reduced, the acquisition efficiency is improved, and the multi-dimensional NMR response mechanism of saturated oil, water and CO2 rocks is economically and efficiently studied.
Smart Images

Figure CN116297615B_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the technical field of CO2 geological storage, and more specifically, to a method, apparatus, device and readable storage medium for acquiring a multidimensional nuclear magnetic resonance spectrum. Background Art
[0002] Geological storage of CO2 is currently one of the most economical and effective methods for controlling CO2 emissions. This involves long-term storage of captured CO2 in geological formations such as coal gas reservoirs, oil and gas fields, and saline aquifers, mitigating its impact on global climate. Abandoned oil and gas fields offer one of the greatest potential locations for CO2 storage, making research on their geological storage crucial.
[0003] Nuclear magnetic resonance (NMR) technology can effectively detect hydrogen-containing fluid signals, such as oil, gas, and water, in rock pores. It can also indirectly assess CO2 stored in formations. However, the signals of oil, gas, and water in rocks often overlap in one-dimensional NMR spectra, making the identification and quantitative evaluation of formation pore fluids difficult. Multidimensional NMR technology can simultaneously identify and quantitatively evaluate formation pore fluids based on the difference between relaxation and diffusion information.
[0004] In order to effectively monitor the amount of CO2 stored in rocks, it is necessary to clarify the response mechanism of multidimensional nuclear magnetic resonance (NMR) of rocks saturated with oil, water, and CO2. This allows the storage amount of CO2 to be determined during the process of CO2 storage in rocks and enables quantitative evaluation of CO2. Determining the response mechanism of multidimensional NMR of rocks saturated with oil, water, and CO2 can be achieved by obtaining multidimensional NMR spectra and analyzing their characteristics. Currently, multidimensional NMR spectra are generally obtained by NMR of rocks saturated with oil, water, and CO2. However, this experiment not only requires continuous injection of CO2 into actual rocks saturated with oil and water, but also requires a lot of related equipment. Therefore, it is difficult to implement, and the acquisition efficiency of multidimensional NMR spectra is low and the cost is high.
[0005] In summary, how to obtain multidimensional nuclear magnetic resonance spectra efficiently and at low cost is a technical problem that needs to be solved urgently by those skilled in the art. Summary of the Invention
[0006] In view of this, the purpose of this application is to provide a method, device, equipment and readable storage medium for obtaining multidimensional nuclear magnetic resonance spectra, which are used to obtain multidimensional nuclear magnetic resonance spectra efficiently and at low cost.
[0007] In order to achieve the above objectives, this application provides the following technical solutions:
[0008] A method for obtaining a multidimensional nuclear magnetic resonance spectrum, comprising:
[0009] Obtaining a rock CT image, and constructing a three-dimensional digital rock core based on the rock CT image;
[0010] Based on the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions;
[0011] According to the three-dimensional digital core model of the rock, a multi-dimensional nuclear magnetic resonance acquisition process is simulated by a numerical simulation method to obtain multi-dimensional nuclear magnetic resonance echo data;
[0012] The multidimensional nuclear magnetic resonance echo data is inverted to obtain a corresponding multidimensional nuclear magnetic resonance spectrum.
[0013] Preferably, based on the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2, including:
[0014] The saturated oil and water are regarded as fluids in the same phase, and the distribution state of CO2 and oil-water mixture in the three-dimensional digital core is simulated using the opening operation in the morphological method to obtain an intermediate model;
[0015] In the intermediate model, CO2 and the skeleton are regarded as the same phase medium, and the opening operation in the morphological method is used to simulate the distribution state of oil and water to obtain the distribution state of the three-phase fluid of CO2, oil and water. The three-phase fluids of CO2, oil and water are replaced by different numbers respectively to obtain the three-dimensional digital core model of the rock.
[0016] Preferably, when the multidimensional nuclear magnetic resonance is T1-D-T2 three-dimensional nuclear magnetic resonance, simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multidimensional nuclear magnetic resonance echo data includes:
[0017] According to the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase signal recovery of water in the rock pores, the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase signal attenuation of water, and the multi-dimensional nuclear magnetic resonance T2 acquisition phase signal attenuation of water are sequentially simulated by the numerical simulation method to obtain multiple echo signal intensities of water in the T2 acquisition phase;
[0018] Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multidimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil, the signal attenuation of the multidimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil, and the signal attenuation of the multidimensional nuclear magnetic resonance T2 acquisition phase of the oil in the rock pores, thereby obtaining multiple echo signal intensities of the oil in the T2 acquisition phase;
[0019] The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0020] Changing the waiting time of the T1 editing stage and / or the time interval between two 180° pulses in the D editing stage, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multidimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multidimensional nuclear magnetic resonance diffusion coefficient D editing stage of water, and the signal attenuation of the multidimensional nuclear magnetic resonance T2 acquisition stage of water in the rock pores by the numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until N groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where N is an integer greater than 1;
[0021] N groups of nuclear magnetic resonance echo data of rocks saturated with oil, water and CO2 are used as the multi-dimensional nuclear magnetic resonance echo data.
[0022] Preferably, when the multidimensional nuclear magnetic resonance is T1-T2 two-dimensional nuclear magnetic resonance, simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multidimensional nuclear magnetic resonance echo data includes:
[0023] Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water to obtain multiple echo signal intensities of water in the T2 acquisition phase;
[0024] Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil to obtain multiple echo signal intensities of the oil in the T2 acquisition phase;
[0025] The intensities of each echo signal of water in the T2 acquisition phase and the intensities of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0026] Changing the waiting time of the T1 editing stage, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition stage of water, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until M groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where M is an integer greater than 1;
[0027] M groups of oil-, water- and CO2-saturated rock nuclear magnetic resonance echo data are used as the multidimensional nuclear magnetic resonance echo data.
[0028] Preferably, when the multidimensional nuclear magnetic resonance is D-T2 two-dimensional nuclear magnetic resonance, simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multidimensional nuclear magnetic resonance echo data includes:
[0029] According to the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance diffusion coefficient D of water in the rock pores is sequentially simulated by the numerical simulation method. The signal attenuation during the editing phase and the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water are simulated to obtain multiple echo signal intensities of water during the T2 acquisition phase.
[0030] Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil in the rock pores to obtain multiple echo signal intensities of the oil in the T2 acquisition phase;
[0031] The intensities of each echo signal of water in the T2 acquisition phase and the intensities of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0032] changing the time interval between the two 180° pulses in the D editing phase, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient of water in the rock pores in the D editing phase and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water by the numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition phase, until P groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where P is an integer greater than 1;
[0033] The P group of oil-, water- and CO2-saturated rock nuclear magnetic resonance echo data is used as the multi-dimensional nuclear magnetic resonance echo data.
[0034] Preferably, after constructing the three-dimensional digital core according to the rock CT image, the method further includes:
[0035] Based on the three-dimensional digital core, a three-dimensional digital core model of rock containing only saturated oil and water is established using the morphological method.
[0036] Preferably, based on the three-dimensional digital core, a three-dimensional digital core model of rock containing only saturated oil and water is established using the morphological method, including:
[0037] Assuming that the pores of the three-dimensional digital core contain only oil and water wetting phase fluids, the pore space is subjected to the opening operation in the morphological method to simulate the non-wetting phase fluids containing oil and water, thereby obtaining a three-dimensional digital core model of the rock containing only oil and water.
[0038] A multi-dimensional nuclear magnetic resonance spectrum acquisition device, comprising:
[0039] A construction module is used to obtain a rock CT image and construct a three-dimensional digital rock core based on the rock CT image;
[0040] The first establishment module is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions using a morphological method based on the three-dimensional digital core;
[0041] A simulation module is used to simulate the acquisition process of multi-dimensional nuclear magnetic resonance by a numerical simulation method based on the three-dimensional digital core model of the rock, and obtain multi-dimensional nuclear magnetic resonance echo data;
[0042] The inversion module is used to invert the multidimensional nuclear magnetic resonance echo data to obtain a corresponding multidimensional nuclear magnetic resonance spectrum.
[0043] A multi-dimensional nuclear magnetic resonance spectrum acquisition device, comprising:
[0044] memory for storing computer programs;
[0045] A processor is configured to implement the steps of any of the above-described methods for acquiring multi-dimensional nuclear magnetic resonance spectra when executing the computer program.
[0046] A readable storage medium stores a computer program, which, when executed by a processor, implements the steps of any of the above-mentioned methods for acquiring multi-dimensional nuclear magnetic resonance spectra.
[0047] The present application provides a method, apparatus, device, and readable storage medium for acquiring a multidimensional nuclear magnetic resonance spectrum, wherein the method comprises: acquiring a rock CT image, and constructing a three-dimensional digital rock core based on the rock CT image; using a morphological method to establish a three-dimensional digital rock core model of a rock saturated with oil, water, and CO2 under different saturation conditions based on the three-dimensional digital rock core; simulating a multidimensional nuclear magnetic resonance acquisition process using a numerical simulation method based on the three-dimensional digital rock core model to acquire multidimensional nuclear magnetic resonance echo data; and inverting the multidimensional nuclear magnetic resonance echo data to obtain a corresponding multidimensional nuclear magnetic resonance spectrum.
[0048] Compared to the existing experimental methods for obtaining multidimensional nuclear magnetic resonance spectra, the technical solution disclosed in this application uses a morphological method to establish a three-dimensional digital core model of rocks saturated with oil, water, and CO2 under different saturation conditions to simulate the process of CO2 injection into oil- and water-saturated rocks. A simulation method is then used to simulate the multidimensional nuclear magnetic resonance acquisition process to obtain multidimensional nuclear magnetic resonance echo data. The echo data is then inverted to obtain a multidimensional nuclear magnetic resonance spectrum. In other words, this application uses numerical simulation to obtain a multidimensional nuclear magnetic resonance spectrum during the CO2 injection process into oil- and water-saturated rocks. Since numerical simulation does not require actual CO2 injection experiments into oil- and water-saturated rocks and does not require various experimental-related equipment, it can reduce the complexity and cost of obtaining multidimensional nuclear magnetic resonance spectra, improve the efficiency of obtaining multidimensional nuclear magnetic resonance spectra, and achieve economical and efficient acquisition of multidimensional nuclear magnetic resonance spectra, thereby achieving economical and efficient research on the multidimensional nuclear magnetic resonance response mechanism of rocks saturated with oil, water, and CO2. BRIEF DESCRIPTION OF THE DRAWINGS
[0049] In order to more clearly illustrate the embodiments of the present application 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 merely embodiments of the present application. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without any creative work.
[0050] Figure 1 A flowchart of a method for acquiring a multidimensional nuclear magnetic resonance spectrum provided in an embodiment of the present application;
[0051] Figure 2 A three-dimensional digital core model of carbonate rock obtained from a CT experiment provided in an embodiment of the present application;
[0052] Figure 3 A digital core model diagram of a carbonate rock with an oil saturation of 51.22% and a water saturation of 48.78% provided in an embodiment of the present application;
[0053] Figure 4 A digital core model diagram of a carbonate rock with an oil saturation of 40.93%, a water saturation of 43.50%, and a CO2 saturation of 15.57% provided in the embodiment of the present application;
[0054] Figure 5 A digital core model diagram of a carbonate rock with an oil saturation of 27.19%, a water saturation of 36.71%, and a CO2 saturation of 36.10% provided in the embodiment of the present application;
[0055] Figure 6A digital core model diagram of a carbonate rock with an oil saturation of 19.92%, a water saturation of 28.86%, and a CO2 saturation of 51.22% provided in the embodiment of the present application;
[0056] Figure 7 For the simulation Figure 3 D-T2 spectrum of carbonate digital core;
[0057] Figure 8 For the simulation Figure 4 D-T2 spectrum of carbonate digital core in;
[0058] Figure 9 For the simulation Figure 5 D-T2 spectrum of carbonate digital core;
[0059] Figure 10 For the simulation Figure 6 D-T2 spectrum of carbonate digital core in;
[0060] Figure 11 A schematic structural diagram of a multi-dimensional nuclear magnetic resonance spectrum acquisition device provided in an embodiment of the present application;
[0061] Figure 12 A schematic diagram of the structure of a multi-dimensional nuclear magnetic resonance spectrum acquisition device provided in an embodiment of the present application. DETAILED DESCRIPTION
[0062] The core of this application is to provide a method, device, equipment and readable storage medium for obtaining multidimensional nuclear magnetic resonance spectra, which are used to obtain multidimensional nuclear magnetic resonance spectra efficiently and at low cost.
[0063] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by ordinary technicians in this field without making creative efforts are within the scope of protection of this application.
[0064] See also Figure 1 , which shows a flow chart of a method for obtaining a multidimensional nuclear magnetic resonance spectrum provided in an embodiment of the present application. The method for obtaining a multidimensional nuclear magnetic resonance spectrum provided in an embodiment of the present application may include:
[0065] S11: Obtain a rock CT image and construct a three-dimensional digital rock core based on the rock CT image.
[0066] In this application, a rock CT (Computed Tomography) experiment is first conducted on the rock to be used for CO2 geological storage to obtain the corresponding rock CT images. Then, a 3D digital rock core (i.e., a 3D rock model) is constructed based on the rock CT images. This 3D digital rock core model, saturated with oil, water, and CO2, is then constructed. The constructed 3D digital rock core only contains the skeleton and pores.
[0067] S12: Based on the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rocks saturated with oil, water and CO2 under different saturation conditions.
[0068] After constructing the 3D digital core, morphological methods can be used to build 3D digital core models of rocks saturated with oil, water, and CO2 at different saturation levels to simulate the injection of CO2 into oil- and water-saturated rocks. During CO2 injection into oil- and water-saturated rocks, the CO2 saturation gradually increases over time.
[0069] That is to say, multiple three-dimensional digital core models of rocks are established. In these multiple three-dimensional digital core models of rocks, the saturations of saturated oil, saturated water and saturated CO2 are different, so as to well simulate the saturation changes of the three-phase fluids of oil, water and CO2 in the saturated oil-water rock during the process of CO2 injection into the oil-water-saturated rock, so that the established three-dimensional digital core model of rock is more in line with the actual situation, thereby improving the accuracy of the multi-dimensional nuclear magnetic resonance spectrum finally obtained.
[0070] S13: Based on the three-dimensional digital core model of the rock, a multi-dimensional nuclear magnetic resonance acquisition process is simulated by a numerical simulation method to obtain multi-dimensional nuclear magnetic resonance echo data.
[0071] On the basis of step S12, according to the established three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions, the multi-dimensional nuclear magnetic resonance acquisition process is simulated by a numerical simulation method, and multi-dimensional nuclear magnetic resonance echo data in the simulated multi-dimensional nuclear magnetic resonance acquisition process is obtained, so as to obtain the corresponding multi-dimensional nuclear magnetic resonance spectrum according to the multi-dimensional nuclear magnetic resonance echo data.
[0072] The multidimensional nuclear magnetic resonance mentioned above may specifically be two-dimensional nuclear magnetic resonance, three-dimensional nuclear magnetic resonance, etc.
[0073] S14: Invert the multi-dimensional nuclear magnetic resonance echo data to obtain a corresponding multi-dimensional nuclear magnetic resonance spectrum.
[0074] After acquiring multidimensional NMR echo data, the multidimensional NMR echo data can be inverted to obtain a corresponding multidimensional NMR spectrum, thereby facilitating the study of the NMR response mechanism of rocks saturated with oil, water, and CO2 based on the multidimensional NMR spectrum. Specifically, after obtaining the corresponding multidimensional NMR spectrum, the multidimensional NMR spectrum characteristics under different saturation conditions and different acquisition parameters (specifically, the acquisition parameters used in the simulated multidimensional NMR acquisition process) are analyzed based on the multidimensional NMR spectrum to obtain the multidimensional NMR response mechanism of rocks saturated with oil, water, and CO2.
[0075] It can be seen from this that compared with the existing experimental method for obtaining multidimensional nuclear magnetic resonance spectra, the present application realizes the acquisition of multidimensional nuclear magnetic resonance spectra during the CO2 injection into oil-water saturated rocks through numerical simulation. Therefore, the present application can reduce the complexity and cost of obtaining multidimensional nuclear magnetic resonance spectra during the CO2 injection into oil-water saturated rocks, and can improve the efficiency of obtaining multidimensional nuclear magnetic resonance spectra during the CO2 injection into oil-water saturated rocks, thereby realizing low-cost and efficient research on the multidimensional nuclear magnetic resonance response mechanism of saturated oil, water and CO2 rocks.
[0076] Compared to the existing experimental methods for obtaining multidimensional nuclear magnetic resonance spectra, the technical solution disclosed in this application uses a morphological method to establish a three-dimensional digital core model of rocks saturated with oil, water, and CO2 under different saturation conditions to simulate the process of CO2 injection into oil- and water-saturated rocks. A simulation method is then used to simulate the multidimensional nuclear magnetic resonance acquisition process to obtain multidimensional nuclear magnetic resonance echo data. The echo data is then inverted to obtain a multidimensional nuclear magnetic resonance spectrum. In other words, this application uses numerical simulation to obtain a multidimensional nuclear magnetic resonance spectrum during the CO2 injection process into oil- and water-saturated rocks. Since numerical simulation does not require actual CO2 injection experiments into oil- and water-saturated rocks and does not require various experimental-related equipment, it can reduce the complexity and cost of obtaining multidimensional nuclear magnetic resonance spectra, improve the efficiency of obtaining multidimensional nuclear magnetic resonance spectra, and achieve economical and efficient acquisition of multidimensional nuclear magnetic resonance spectra, thereby achieving economical and efficient research on the multidimensional nuclear magnetic resonance response mechanism of rocks saturated with oil, water, and CO2.
[0077] The present invention provides a method for acquiring a multidimensional nuclear magnetic resonance spectrum, which uses a morphological method to establish a three-dimensional digital core model of a rock saturated with oil, water, and CO2 based on a three-dimensional digital core. The method may include:
[0078] The saturated oil and water are considered as the same phase fluid, and the distribution of CO2 and oil-water mixture in the 3D digital core is simulated by using the opening operation in the morphological method to obtain the intermediate model.
[0079] In the intermediate model, CO2 and the skeleton are regarded as the same-phase medium, and the opening operation in the morphological method is used to simulate the distribution state of oil and water to obtain the distribution state of the three-phase fluid of CO2, oil and water. The three-phase fluids of CO2, oil and water are replaced by different numbers respectively to obtain the three-dimensional digital core model of the rock.
[0080] In this application, when CO2 is injected into oil- and water-saturated rock, the CO2 occupies the central region of the pores, while the non-wetting phase fluids in the oil and water are distributed in the secondary central region, and the wetting phase fluids in the oil and water are distributed near the pore walls. The fluid that easily adheres to the rock is the wetting phase fluid, while the fluid that does not adhere to the rock is the non-wetting phase fluid. For water-wet rock, water is the wetting phase fluid and oil is the non-wetting phase fluid; for oil-wet rock, oil is the wetting phase fluid and water is the non-wetting phase fluid.
[0081] On the basis of the above, the specific process of establishing a three-dimensional digital core model of a rock saturated with oil, water, and CO2 under different saturation conditions using a morphological method based on the three-dimensional digital core is as follows: first, oil and water are regarded as the same-phase fluid (the same-phase fluid is an oil-water mixture), and the opening operation in the morphological method is used to simulate the distribution state of CO2 and the oil-water mixture in the digital core (specifically, the three-dimensional digital core constructed in step S11) (specifically, the simulation is performed according to the saturation conditions of the three-phase fluids saturated with oil, water, and CO2). Thereafter, CO2 and the skeleton in the digital core are regarded as the same-phase medium, and the opening operation in the morphological method is again used to simulate the distribution state of oil and water in the digital core (i.e., oil and water are separated, and the simulation is performed according to the saturation conditions of saturated oil and water). Finally, the distribution state of the three-phase fluids of CO2, oil, and water in the digital core is obtained. Different fluids are replaced by different numbers to obtain a three-dimensional digital core model of a rock saturated with oil, water, and CO2.
[0082] Specifically, the above process involves establishing a three-dimensional digital core model of rock saturated with oil, water, and CO2 under one saturation condition. Specifically, the above process can be performed after determining the saturations of oil, water, and CO2 (i.e., the distribution states of these three-phase fluids) to establish a three-dimensional digital core model of rock saturated with oil, water, and CO2 under one saturation condition. The process for establishing three-dimensional digital core models of rock saturated with oil, water, and CO2 under other saturation conditions is the same as the above process and will not be repeated here. It should be noted that the three-dimensional digital core models of rock saturated with oil, water, and CO2 under different saturation conditions can be established simultaneously or sequentially.
[0083] The present invention provides a method for acquiring a multidimensional nuclear magnetic resonance spectrum. When the multidimensional nuclear magnetic resonance is a three-dimensional T1-D-T2 nuclear magnetic resonance, the method may include: simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on a three-dimensional digital core model of the rock to acquire multidimensional nuclear magnetic resonance echo data;
[0084] Based on the three-dimensional digital core model of rock, numerical simulation methods are used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of water, and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water, and the multiple echo signal intensities of water in the T2 acquisition phase are obtained.
[0085] Based on the three-dimensional digital core model of the rock, numerical simulation methods are used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil, and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil, and the multiple echo signal intensities of the oil in the T2 acquisition phase are obtained.
[0086] The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0087] Changing the waiting time of the T1 editing stage and / or the time interval between two 180° pulses in the D editing stage, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing stage of water, and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition stage of water in the rock pores by a numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until N sets of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where N is an integer greater than 1;
[0088] N groups of nuclear magnetic resonance echo data of rocks saturated with oil, water and CO2 are used as multi-dimensional nuclear magnetic resonance echo data.
[0089] In the present application, when the multidimensional nuclear magnetic resonance is T1-D-T2 three-dimensional nuclear magnetic resonance, the acquisition process of the multidimensional nuclear magnetic resonance is simulated by a numerical simulation method according to the three-dimensional digital core model of the rock, and the process of obtaining the multidimensional nuclear magnetic resonance echo data can be: according to the three-dimensional digital core model of the rock saturated with oil, water and CO2 under different saturation conditions established in step S12, the multidimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage signal recovery of water in the rock pores is first simulated by a numerical simulation method, and then the multidimensional nuclear magnetic resonance diffusion coefficient D editing stage signal attenuation of water in the rock pores is simulated, and then the multidimensional nuclear magnetic resonance T2 acquisition stage signal attenuation of water in the rock pores is simulated to obtain multiple echo signal intensities of water in the T2 acquisition stage. Similarly, a multidimensional NMR simulation of oil in pores is performed. Specifically, based on the three-dimensional digital core model of the rock saturated with oil, water, and CO2 at different saturation conditions established in step S12, a numerical simulation method is used to first simulate the signal recovery during the T1 editing phase of the multidimensional NMR longitudinal relaxation time of the oil in the rock pores. Then, the signal attenuation during the T2 acquisition phase of the multidimensional NMR diffusion coefficient of the oil in the rock pores is simulated. Finally, the signal attenuation during the T2 acquisition phase of the multidimensional NMR acquisition phase of the oil in the rock pores is simulated to obtain multiple echo signal intensities of the oil in the T2 acquisition phase. Then, the oil and water signals at the corresponding acquisition moments are added together to obtain the NMR echo data of the rock saturated with oil, water, and CO2. Specifically, the intensities of each echo signal of water in the T2 acquisition stage and the intensities of each echo signal of oil in the T2 acquisition stage are added together to obtain a set of nuclear magnetic resonance echo data of rock saturated with oil, water, and CO2. Then, the waiting time of the T1 editing stage and / or the time interval between two 180° pulses in the D editing stage are changed, and the above steps are repeated until N (N is an integer greater than 1) sets of nuclear magnetic resonance echo data of rock saturated with oil, water, and CO2 are obtained. The obtained N sets of nuclear magnetic resonance echo data of rock saturated with oil, water, and CO2 are multidimensional nuclear magnetic resonance echo data.
[0090] Taking the rock as a hydrophilic rock as an example, the specific steps of the above process are as follows:
[0091] ① For the multi-dimensional NMR response simulation of water in pores, first randomly place n particles with an initial signal intensity of 1 in the voxels representing water in the three-dimensional digital core model of the rock (i.e., at the digital position representing water, there is at most one particle in each voxel), and set the diffusion radius r, which is no more than 1 / 3 of the resolution set by the CT experiment. The time step of the simulation is calculated based on the set diffusion radius. , is the diffusion coefficient of water.
[0092] ② These n particles move randomly within each time step. When the particles collide with the rock skeleton during movement, the signal intensity of the particles will weaken, and the signal intensity will decay to (1-p) times the original value after each collision. When the particles collide with other pore fluids (oil or CO2) during movement, the signal intensity of the particles remains unchanged, and the particles return to the position before the collision.
[0093] ③ Define the time of T1 editing stage as waiting time T W , then T W The signal recovery of the i-th particle after time can be expressed as: , where is the longitudinal bulk relaxation time of water, T W The remaining signal intensity of the i-th particle after colliding with the rock skeleton after time T W The number of collisions of the particle in time is x, then ,in, is the longitudinal surface relaxation strength of the rock.
[0094] ④ The time interval between the two 180° pulses in the diffusion coefficient D editing phase is defined as T E1 , then the total signal strength of n particles after the D editing stage is completed can be expressed as: , where is the transverse bulk relaxation time of water, represents the phase shift of the i-th particle due to spin after the D editing stage, is the remaining signal intensity of the i-th particle after it collides with the rock skeleton after the D editing stage. If 2 T E1 The number of collisions of the particle in time is y, then ,in is the longitudinal surface relaxation strength of the rock.
[0095] ⑤ Assume that the T2 acquisition phase contains m 180° pulses, and the time interval between two adjacent 180° pulses is T E , then the signal strength of the jth (1 ≤ j ≤ m) echo in the T2 acquisition phase can be expressed as: , where is the remaining signal strength after the ith particle collides with the rock skeleton in the first j echoes. If the number of collisions of the particle in the first j echoes is z, then In order to simulate the signal intensity ratio of water and oil in relation to their saturation ratio, the simulated echo signal is normalized and then multiplied by the corresponding water porosity value, which is (that is, the multiple echo signal intensities of water in the T2 acquisition stage are): , where is the water-containing porosity.
[0096] ⑥ The multi-dimensional NMR response simulation of oil in pores is consistent with that of water in pores. The difference is that oil, as a non-wetting phase fluid, does not collide with the rock skeleton, so its signal intensity is not weakened by collision. Therefore, the signal intensity of the jth (1 ≤ j ≤ m) echo in the T2 acquisition phase (i.e., the multiple echo signal intensities of oil in the T2 acquisition phase) can be expressed as: , where is the longitudinal bulk relaxation time of the oil, is the transverse bulk relaxation time of the oil, is the oil-bearing porosity, is the number of particles in the voxel representing oil in the three-dimensional digital core model of rock at the initial moment.
[0097] ⑦ Add the oil and water echo signals simulated in steps ⑤ and ⑥ to obtain the fixed T w and T E1 saturated oil, water and CO2 rock nuclear magnetic resonance echo data under conditions (i.e., a set of saturated oil, water and CO2 rock nuclear magnetic resonance echo data, this set of saturated oil, water and CO2 rock nuclear magnetic resonance echo data contains m echo data).
[0098] ⑧Then change T w and / or T E1 , repeat steps ③-⑦ to get different T w 、T E1 The nuclear magnetic resonance echo data of rocks saturated with oil, water and CO2 under the above conditions are obtained, that is, N groups of nuclear magnetic resonance echo data of rocks saturated with oil, water and CO2 are obtained, which are multi-dimensional nuclear magnetic resonance echo data.
[0099] It should be noted that the above description uses water-wetting rock as an example. For oil-wetting rock, oil collides with the rock skeleton, while water does not. This means that the echo signal intensity for oil is determined according to steps 1-5, while the echo signal intensity for water is determined according to step 6. In other words, for rock, the wetting phase fluids (oil and water) collide with the rock skeleton, and multiple echo signal intensities for the wetting phase fluids (oil and water) during the T2 acquisition phase are determined according to steps 1-5. The non-wetting phase fluids (oil and water) do not collide with the rock skeleton, and multiple echo signal intensities for the non-wetting phase fluids (oil and water) during the T2 acquisition phase are determined according to step 6.
[0100] The present invention provides a method for acquiring a multidimensional nuclear magnetic resonance spectrum. When the multidimensional nuclear magnetic resonance is a two-dimensional T1-T2 nuclear magnetic resonance, the method may include: simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on a three-dimensional digital core model of the rock to acquire multidimensional nuclear magnetic resonance echo data;
[0101] Based on the three-dimensional digital core model of rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water are simulated in sequence through numerical simulation methods, and the multiple echo signal intensities of water in the T2 acquisition phase are obtained.
[0102] Based on the three-dimensional digital core model of the rock, the signal recovery during the T1 editing phase of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time of oil in the rock pores and the signal attenuation during the T2 acquisition phase of the multi-dimensional nuclear magnetic resonance of oil are simulated in sequence through numerical simulation methods, and the multiple echo signal intensities of the oil in the T2 acquisition phase are obtained.
[0103] The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0104] Changing the waiting time of the T1 editing stage, returning to the execution of the steps of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition stage of water, and obtaining the multiple echo signal intensities of water in the T2 acquisition stage, in sequence, by a numerical simulation method, until M groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained; M is an integer greater than 1;
[0105] The M groups of oil-, water- and CO2-saturated rock NMR echo data are used as multi-dimensional NMR echo data.
[0106] In the present application, when the multidimensional NMR is T1-T2 two-dimensional NMR, compared to T1-D-T2 three-dimensional NMR, there is no need to simulate the signal attenuation during the D editing phase. That is, based on the three-dimensional digital core model of the rock, the multidimensional NMR acquisition process is simulated by a numerical simulation method. The process of obtaining multidimensional NMR echo data can specifically be as follows: based on the three-dimensional digital core model of the rock saturated with oil, water, and CO2 under different saturation conditions established in step S12, the multidimensional NMR longitudinal relaxation time T1 editing phase signal recovery of water in the rock pores is first simulated by a numerical simulation method, and then the multidimensional NMR T2 acquisition phase signal attenuation of water in the rock pores is simulated to obtain multiple echo signal intensities of water in the T2 acquisition phase. Similarly, a multidimensional NMR simulation of oil in pores is performed. Specifically, based on the three-dimensional digital core model of the rock saturated with oil, water, and CO2 at different saturation conditions established in step S12, a numerical simulation method is used to first simulate the signal recovery during the T1 editing phase of the multidimensional NMR longitudinal relaxation time of oil in the rock pores. Then, the signal attenuation during the T2 acquisition phase of the multidimensional NMR of oil in the rock pores is simulated to obtain the multiple echo signal intensities of the oil during the T2 acquisition phase. Then, the oil and water signals at the corresponding acquisition times are added to obtain the NMR echo data of the rock saturated with oil, water, and CO2. Specifically, the intensities of each echo signal of water in the T2 acquisition stage and the intensities of each echo signal of oil in the T2 acquisition stage are added together to obtain a set of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2. Then, the waiting time of the T1 editing stage is changed, and the above steps are repeated until M (M is an integer greater than 1) sets of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2 are obtained. The obtained M sets of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2 are multidimensional nuclear magnetic resonance echo data.
[0107] Taking the rock as a hydrophilic rock as an example, the specific steps of the above process are similar to those of ①-⑧ above. The difference is that there is no ④ step in the T1-T2 two-dimensional nuclear magnetic resonance simulation process, and only T is changed in ⑧. w , T E1 Considered as 0. If it is an oil-wet rock, the oil collides with the rock skeleton, and the water does not collide with the rock skeleton. Combining the above-mentioned differences, the oil is obtained by following steps ①-⑤ (without step ④) to obtain the echo signal strength (T E1 Considered as 0), water is to obtain the echo signal strength according to step ⑥ (T E1 That is to say, for the rock, the wetting phase fluid in oil and water collides with the rock skeleton, and according to steps ①-⑤ (without step ④), multiple echo signal intensities of the wetting phase fluid in oil and water in the T2 acquisition stage are obtained (T E1The non-wetting phase fluids in oil and water do not collide with the rock skeleton. According to step ⑥, the multiple echo signal intensities of the non-wetting phase fluids in oil and water in the T2 acquisition stage are obtained (T E1 treated as 0).
[0108] The present invention provides a method for acquiring a multidimensional nuclear magnetic resonance spectrum. When the multidimensional nuclear magnetic resonance is a D-T2 two-dimensional nuclear magnetic resonance, the method may include: simulating the acquisition process of the multidimensional nuclear magnetic resonance by a numerical simulation method based on a three-dimensional digital core model of the rock to acquire multidimensional nuclear magnetic resonance echo data;
[0109] Based on the three-dimensional digital core model of rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D of water in rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water are simulated in sequence through numerical simulation methods, and the multiple echo signal intensities of water in the T2 acquisition phase are obtained.
[0110] Based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing stage of the oil in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition stage of the oil are simulated in sequence through numerical simulation methods, and the multiple echo signal intensities of the oil in the T2 acquisition stage are obtained.
[0111] The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0112] changing the time interval between the two 180° pulses in the D editing phase, and returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient of water in the rock pores in the D editing phase and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water in the rock pores by a numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition phase, until P groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained; P is an integer greater than 1;
[0113] The nuclear magnetic resonance echo data of the P group of rocks saturated with oil, water and CO2 are used as multi-dimensional nuclear magnetic resonance echo data.
[0114] In the present application, when the multidimensional NMR is D-T2 two-dimensional NMR, compared to T1-D-T2 three-dimensional NMR, there is no need to simulate the signal attenuation during the T1 editing phase. That is, based on the three-dimensional digital core model of the rock, the multidimensional NMR acquisition process is simulated by a numerical simulation method. The process of obtaining multidimensional NMR echo data can specifically be as follows: based on the three-dimensional digital core model of the rock saturated with oil, water, and CO2 under different saturation conditions established in step S12, the multidimensional NMR diffusion coefficient D editing phase signal attenuation of water in the rock pores is first simulated by a numerical simulation method, and then the multidimensional NMR T2 acquisition phase signal attenuation of water in the rock pores is simulated to obtain multiple echo signal intensities of water in the T2 acquisition phase. Similarly, a multidimensional NMR simulation of oil in pores is performed. Specifically, based on the three-dimensional digital core model of the rock saturated with oil, water, and CO2 at different saturation conditions established in step S12, a numerical simulation method is used to first simulate the signal attenuation during the D editing phase of the multidimensional NMR diffusion coefficient of oil in the rock pores. Then, the signal attenuation during the T2 acquisition phase of the multidimensional NMR of oil in the rock pores is simulated to obtain multiple echo signal intensities of the oil during the T2 acquisition phase. Then, the oil and water signals at corresponding acquisition times are added to obtain the NMR echo data of the rock saturated with oil, water, and CO2. Specifically, the intensities of each echo signal of water in the T2 acquisition stage and the intensities of each echo signal of oil in the T2 acquisition stage are added together to obtain a set of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2. Then, the time interval between the two 180° pulses in the D editing stage is changed, and the above steps are repeated until P (P is an integer greater than 1) sets of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2 are obtained. The obtained P sets of nuclear magnetic resonance echo data of rocks saturated with oil, water, and CO2 are the multidimensional nuclear magnetic resonance echo data.
[0115] Taking the rock as a hydrophilic rock as an example, the specific steps of the above process are similar to those of ①-⑧ above. The difference is that there is no step ③ in the D-T2 two-dimensional nuclear magnetic resonance simulation process, and only the T E1 , T W Considered to be infinite. If it is an oil-wet rock, the oil collides with the rock skeleton, while water does not. Combining the above differences, the oil is obtained by following steps ①-⑤ (without step ③) to obtain the echo signal intensity (T W Considered as infinite), water is to obtain the echo signal strength according to step ⑥ (T W In other words, for the rock, the wetting phase fluids in oil and water collide with the rock skeleton, and according to steps ①-⑤ (without step ③), multiple echo signal intensities of the wetting phase fluids in oil and water in the T2 acquisition stage are obtained (T WConsidered as infinite), the non-wetting phase fluids in oil and water do not collide with the rock skeleton, and the multiple echo signal intensities of the non-wetting phase fluids in oil and water in the T2 acquisition stage are obtained according to step ⑥ (T W is considered infinite).
[0116] The multi-dimensional nuclear magnetic resonance spectrum acquisition method provided in an embodiment of the present application may further include, after constructing a three-dimensional digital rock core based on the rock CT image:
[0117] Based on the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rock containing only saturated oil and water.
[0118] In the present application, after constructing a three-dimensional digital core based on the rock CT image, a three-dimensional digital core model of the rock containing only saturated oil and water can also be established based on the three-dimensional digital core using a morphological method, that is, a three-dimensional digital core model of the rock containing only saturated oil and water without CO2 injection is constructed, that is, a three-dimensional digital core model of the rock in the initial state is constructed, so as to study the multi-dimensional nuclear magnetic resonance response mechanism of the rock saturated with oil, water and CO2 based on the three-dimensional digital core model of the rock containing only saturated oil and water and the three-dimensional digital core model of the rock saturated with oil, water and CO2 under different saturation conditions, thereby improving the accuracy of the response mechanism analysis.
[0119] Specifically, a 3D digital core model of rock saturated only with oil and water can be constructed under one saturation condition. This saturation condition can be the same as the oil and water saturation of the rock to be used for CO2 geological storage. Of course, 3D digital core models of rock saturated only with oil and water can also be constructed under multiple saturation conditions.
[0120] The present invention provides a method for acquiring a multidimensional nuclear magnetic resonance spectrum, which uses a morphological method to establish a three-dimensional digital core model of a rock containing only saturated oil and water based on a three-dimensional digital core. The method may include:
[0121] The pores of the three-dimensional digital core are assumed to contain only oil and water as the wetting phase fluid. The pore space is opened using the morphological method to simulate the non-wetting phase fluid containing oil and water, and a three-dimensional digital core model of the rock containing only oil and water is obtained.
[0122] For water-wet rock, when the rock is saturated with oil and water, the non-wetting oil is primarily distributed in the center of the macropores. For oil-wet rock, when the rock is saturated with oil and water, the non-wetting water is primarily distributed in the center of the macropores. Assuming that the pores of the constructed 3D digital core contain only oil and water wetting phase fluids, a morphological opening operation is performed on the core pore space to simulate the non-wetting phase fluids containing oil and water, thereby obtaining a 3D digital core model containing only oil and water.
[0123] Among them, when constructing a three-dimensional digital core model of rock containing only saturated oil and water under various saturation conditions, the above can specifically be to simulate the process of gradually increasing saturation of non-wetting phase fluid in oil and water to obtain a three-dimensional digital core model of rock containing only saturated oil and water under different saturation conditions.
[0124] In order to verify the effectiveness of the above simulation method, this application takes the numerical simulation of the two-dimensional nuclear magnetic resonance D-T2 response during the injection of CO2 into oil-water saturated rocks as an example to illustrate:
[0125] First, a three-dimensional digital core model was constructed based on the carbonate rock micron CT experimental data. Figure 2 As shown in FIG, it shows a three-dimensional digital core model diagram of carbonate rock obtained by CT experiment provided in the embodiment of the present application. When establishing a digital rock model of saturated two-phase fluid, it is assumed that the non-wetting phase fluid is present in the pore center and the wetting phase fluid is present near the pore wall. Therefore, based on the established three-dimensional digital core, a three-dimensional digital core model of saturated oil-water two-phase fluid is established using a morphological method, as shown in FIG. Figure 3 As shown in FIG, it shows a digital core model diagram of carbonate rock with an oil saturation of 51.22% and a water saturation of 48.78% provided by the embodiment of the present application. When CO2 is injected into oil- and water-saturated rock, CO2 will first occupy the center of the pores, while the oil gradually moves closer to the pore walls. Therefore, based on the three-dimensional digital core saturated with oil and water, the morphological method is used again to establish a three-dimensional digital core model saturated with oil, water, and CO2, as shown in FIG. Figure 4-Figure 6 As shown, Figure 4 The digital core model diagram of the carbonate rock with an oil saturation of 40.93%, a water saturation of 43.50% and a CO2 saturation of 15.57% provided in the embodiment of the present application is provided. Figure 5 The digital core model diagram of the carbonate rock with an oil saturation of 27.19%, a water saturation of 36.71% and a CO2 saturation of 36.10% provided in the embodiment of the present application is provided. Figure 6 A digital core model diagram of carbonate rock with an oil saturation of 19.92%, a water saturation of 28.86%, and a CO2 saturation of 51.22% is provided in the embodiment of the present application.
[0126] The simulation method provided in this application is used to simulate the acquisition process of two-dimensional nuclear magnetic resonance D-T2, and obtain Figure 3-Figure 6 The two-dimensional D-T2 data of the carbonate rock digital core are then inverted to obtain the corresponding D-T2 spectrum, such as Figure 7-10 As shown, Figure 7 For the simulation Figure 3 D-T2 spectrum of carbonate digital core in Figure 8For the simulation Figure 4 D-T2 spectrum of carbonate digital core in Figure 9 For the simulation Figure 5 D-T2 spectrum of carbonate digital core in Figure 10 For the simulation Figure 6 D-T2 spectrum of carbonate digital core in. Figure 7-10 The observed porosity of the digital core under different saturation conditions is calculated using the D-T2 spectra in the figure (obtained by dividing the pore volume occupied by oil and water by the total volume of the core).
[0127] according to Figure 7-10 It can be seen that with the increase of CO2 saturation and the decrease of water saturation, the water signal in the D-T2 spectrum gradually moves toward the short relaxation direction, and the overlap with the oil signal on the T2 spectrum becomes more and more obvious until it completely overlaps. At this time, oil and water cannot be identified and quantitatively evaluated based on the T2 spectrum. However, oil and water are better separated on the D spectrum and are not affected by water saturation. This fully demonstrates that two-dimensional nuclear magnetic resonance technology has obvious advantages over one-dimensional nuclear magnetic resonance in fluid identification and quantitative evaluation.
[0128] Figure 7 The digital core contains only oil and water, so the observed porosity calculated from the D-T2 spectrum is the total porosity, which is 23.14%. Figure 7 The observed porosity calculated from the D-T2 spectrum is subtracted from Figure 8 、 Figure 9 、 Figure 10 The observed porosity calculated from the D-T2 spectrum can be used to obtain the porosity saturated with CO2 in the corresponding core, which are 3.62%, 8.37% and 11.87% respectively. Dividing it by the total porosity, the CO2 saturation can be obtained, which are 15.64%, 36.17% and 51.30% respectively. These are very consistent with the CO2 saturation corresponding to the digital core model, verifying the accuracy of the simulation method proposed in this application.
[0129] The present application also provides a multi-dimensional nuclear magnetic resonance spectrum acquisition device, see Figure 11 , which shows a schematic structural diagram of a multi-dimensional nuclear magnetic resonance spectrum acquisition device provided in an embodiment of the present application, which may include:
[0130] Construction module 101 is used to obtain a rock CT image and construct a three-dimensional digital rock core based on the rock CT image;
[0131] The first establishment module 102 is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions using a morphological method based on the three-dimensional digital core;
[0132] The simulation module 103 is used to simulate the acquisition process of the multi-dimensional nuclear magnetic resonance by a numerical simulation method based on the three-dimensional digital core model of the rock, and obtain the multi-dimensional nuclear magnetic resonance echo data;
[0133] The inversion module 104 is configured to invert the multi-dimensional nuclear magnetic resonance echo data to obtain a corresponding multi-dimensional nuclear magnetic resonance spectrum.
[0134] In an embodiment of the present application, a multi-dimensional nuclear magnetic resonance spectrum acquisition device is provided, wherein the first establishing module 102 may include:
[0135] The first simulation unit is used to treat saturated oil and water as the same phase fluid, and use the opening operation in the morphological method to simulate the distribution state of CO2 and oil-water mixture in the three-dimensional digital core to obtain an intermediate model;
[0136] The second simulation unit is used to regard CO2 and the skeleton as the same-phase medium in the intermediate model, and use the opening operation in the morphological method to simulate the distribution state of oil and water to obtain the distribution state of the three-phase fluid of CO2, oil and water. The three-phase fluids of CO2, oil and water are replaced by different numbers respectively to obtain a three-dimensional digital core model of the rock.
[0137] In an embodiment of the present application, a multi-dimensional nuclear magnetic resonance spectrum acquisition device is provided. When the multi-dimensional nuclear magnetic resonance is T1-D-T2 three-dimensional nuclear magnetic resonance, the simulation module 103 may include:
[0138] The third simulation unit is used to simulate, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of water, and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water in the rock pores by a numerical simulation method, and obtain multiple echo signal intensities of water in the T2 acquisition phase;
[0139] The fourth simulation unit is used to simulate, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil, and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil by a numerical simulation method, and obtain multiple echo signal intensities of the oil in the T2 acquisition phase;
[0140] The first obtaining unit is used to add the respective echo signal intensities of water in the T2 acquisition stage and the respective echo signal intensities of oil in the T2 acquisition stage to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0141] The first changing unit is used to change the waiting time of the T1 editing stage and / or the time interval between two 180° pulses in the D editing stage, and return to execute the steps of simulating, based on the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance longitudinal relaxation time of water in the rock pores in the T1 editing stage signal recovery, the multi-dimensional nuclear magnetic resonance diffusion coefficient of water in the D editing stage signal attenuation, and the multi-dimensional nuclear magnetic resonance T2 acquisition stage signal attenuation of water, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until N groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where N is an integer greater than 1;
[0142] The first is used as a unit for treating N groups of oil-saturated, water-saturated and CO2-saturated rock nuclear magnetic resonance echo data as multi-dimensional nuclear magnetic resonance echo data.
[0143] In an embodiment of the present application, a multi-dimensional nuclear magnetic resonance spectrum acquisition device is provided. When the multi-dimensional nuclear magnetic resonance is a T1-T2 two-dimensional nuclear magnetic resonance, the simulation module 103 may include:
[0144] A fifth simulation unit is used to simulate, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water by a numerical simulation method, and obtain multiple echo signal intensities of water in the T2 acquisition phase;
[0145] The sixth simulation unit is used to simulate, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil by a numerical simulation method, and obtain multiple echo signal intensities of the oil in the T2 acquisition phase;
[0146] The second obtaining unit is used to add the respective echo signal intensities of water in the T2 acquisition stage and the respective echo signal intensities of oil in the T2 acquisition stage to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0147] The second changing unit is used to change the waiting time of the T1 editing stage, and return to execute the steps of simulating, based on the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance longitudinal relaxation time of water in the rock pores in the T1 editing stage signal recovery, the multi-dimensional nuclear magnetic resonance T2 acquisition stage signal attenuation of water, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until M groups of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks are obtained, where M is an integer greater than 1;
[0148] The second is used as a unit for treating M groups of oil-saturated, water-saturated and CO2-saturated rock NMR echo data as multi-dimensional NMR echo data.
[0149] In an embodiment of the present application, a multi-dimensional nuclear magnetic resonance spectrum acquisition device is provided. When the multi-dimensional nuclear magnetic resonance is a D-T2 two-dimensional nuclear magnetic resonance, the simulation module 103 may include:
[0150] The seventh simulation unit is used to simulate the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D of water in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water in the rock pores in sequence through a numerical simulation method based on the three-dimensional digital core model of the rock, and obtain multiple echo signal intensities of water in the T2 acquisition phase;
[0151] An eighth simulation unit is used to simulate, based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D of the oil in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil in succession by a numerical simulation method, and obtain multiple echo signal intensities of the oil in the T2 acquisition phase;
[0152] The third obtaining unit is used to add the respective echo signal intensities of water in the T2 acquisition stage and the respective echo signal intensities of oil in the T2 acquisition stage to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks;
[0153] The third changing unit is configured to change the time interval between the two 180° pulses in the D editing phase, and return to executing the steps of simulating, based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient of water in the rock pores in the D editing phase and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water in the rock pores by a numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition phase, until P groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where P is an integer greater than 1;
[0154] The third unit is used to convert the P group of oil-saturated, water-saturated and CO2-saturated rock nuclear magnetic resonance echo data into multi-dimensional nuclear magnetic resonance echo data.
[0155] The multi-dimensional nuclear magnetic resonance spectrum acquisition device provided in the embodiment of the present application may further include:
[0156] The second establishment module is used to establish a three-dimensional digital core model of the rock containing only saturated oil and water based on the three-dimensional digital core after constructing the three-dimensional digital core according to the rock CT image using a morphological method.
[0157] In an embodiment of the present application, a multi-dimensional nuclear magnetic resonance spectrum acquisition device is provided, wherein the second establishment module may include:
[0158] The ninth simulation unit is used to set the pores of the three-dimensional digital core to contain only oil and water as the wetting phase fluid, and perform the opening operation of the pore space in the morphological method to simulate the non-wetting phase fluid containing oil and water, so as to obtain a three-dimensional digital core model of the rock containing only oil and water.
[0159] The present application also provides a multi-dimensional nuclear magnetic resonance spectrum acquisition device, see Figure 12 , which shows a schematic structural diagram of a multi-dimensional nuclear magnetic resonance spectrum acquisition device provided in an embodiment of the present application, which may include:
[0160] Memory 201, used for storing computer programs;
[0161] The processor 202, when used to execute the computer program stored in the memory 201, can implement the following steps:
[0162] Obtain a rock CT image and construct a three-dimensional digital core based on the rock CT image; based on the three-dimensional digital core, use morphological methods to establish a three-dimensional digital core model of the rock saturated with oil, water, and CO2 under different saturation conditions; based on the three-dimensional digital core model, simulate the multi-dimensional nuclear magnetic resonance acquisition process through numerical simulation methods to obtain multi-dimensional nuclear magnetic resonance echo data; invert the multi-dimensional nuclear magnetic resonance echo data to obtain the corresponding multi-dimensional nuclear magnetic resonance spectrum.
[0163] The present application also provides a readable storage medium, which stores a computer program. When the computer program is executed by a processor, the following steps can be implemented:
[0164] Obtain a rock CT image and construct a three-dimensional digital core based on the rock CT image; based on the three-dimensional digital core, use morphological methods to establish a three-dimensional digital core model of the rock saturated with oil, water, and CO2 under different saturation conditions; based on the three-dimensional digital core model, simulate the multi-dimensional nuclear magnetic resonance acquisition process through numerical simulation methods to obtain multi-dimensional nuclear magnetic resonance echo data; invert the multi-dimensional nuclear magnetic resonance echo data to obtain the corresponding multi-dimensional nuclear magnetic resonance spectrum.
[0165] The readable storage medium may include: a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk, etc., which can store program codes.
[0166] The description of the relevant parts of the multidimensional nuclear magnetic resonance spectrum acquisition device, equipment and readable storage medium provided in this application can be found in the detailed description of the corresponding parts of the multidimensional nuclear magnetic resonance spectrum acquisition method provided in the embodiments of this application, and will not be repeated here.
[0167] It should be noted that, in this article, relational terms such as first and second, etc. are only used to distinguish one entity or operation from another entity or operation, and do not necessarily require or imply the existence of any such actual relationship or order between these entities or operations. Moreover, the terms "include", "comprise" or any other variants thereof are intended to cover non-exclusive inclusion, so that the process, method, article or device comprising a series of elements are inherent to the elements. In the absence of further restrictions, the elements limited by the sentence "comprise one..." do not exclude the presence of other identical elements in the process, method, article or device comprising the elements. In addition, the above-mentioned technical solutions provided in the embodiments of the present application are not described in detail in accordance with the corresponding technical solutions in the prior art to achieve the same principle, so as to avoid excessive elaboration.
[0168] The above description of the disclosed embodiments will enable one skilled in the art to implement or use the present application. Various modifications to these embodiments will be readily apparent to one skilled in the art, and the general principles defined herein may be implemented in other embodiments without departing from the spirit or scope of the present application. Therefore, the present application is not limited to the embodiments shown herein, but is intended to conform to the widest scope consistent with the principles and novel features disclosed herein.
Claims
1. A method for obtaining a multidimensional nuclear magnetic resonance spectrum, characterized in that: include: Obtaining a rock CT image, and constructing a three-dimensional digital rock core based on the rock CT image; Based on the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions; According to the three-dimensional digital core, a morphological method is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2, including: The saturated oil and water are regarded as fluids in the same phase, and the distribution state of CO2 and oil-water mixture in the three-dimensional digital core is simulated using the opening operation in the morphological method to obtain an intermediate model; In the intermediate model, CO2 and the skeleton are regarded as the same phase medium, and the distribution state of oil and water is simulated by using the opening operation in the morphological method to obtain the distribution state of the three-phase fluid of CO2, oil and water. The three-phase fluids of CO2, oil and water are replaced by different numbers respectively to obtain the three-dimensional digital core model of the rock; According to the three-dimensional digital core model of the rock, a multi-dimensional nuclear magnetic resonance acquisition process is simulated by a numerical simulation method to obtain multi-dimensional nuclear magnetic resonance echo data; the multi-dimensional nuclear magnetic resonance includes T1-D-T2 three-dimensional nuclear magnetic resonance, T1-T2 two-dimensional nuclear magnetic resonance, or D-T2 two-dimensional nuclear magnetic resonance; The multidimensional nuclear magnetic resonance echo data is inverted to obtain a corresponding multidimensional nuclear magnetic resonance spectrum.
2. The method for obtaining a multidimensional nuclear magnetic resonance spectrum according to claim 1, wherein When the multi-dimensional nuclear magnetic resonance is a T1-D-T2 three-dimensional nuclear magnetic resonance, the acquisition process of the multi-dimensional nuclear magnetic resonance is simulated by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multi-dimensional nuclear magnetic resonance echo data, including: According to the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase signal recovery of water in the rock pores, the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase signal attenuation of water, and the multi-dimensional nuclear magnetic resonance T2 acquisition phase signal attenuation of water are sequentially simulated by the numerical simulation method to obtain multiple echo signal intensities of water in the T2 acquisition phase; Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multidimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil, the signal attenuation of the multidimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil, and the signal attenuation of the multidimensional nuclear magnetic resonance T2 acquisition phase of the oil in the rock pores, thereby obtaining multiple echo signal intensities of the oil in the T2 acquisition phase; The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks; Changing the waiting time of the T1 editing stage and / or the time interval between two 180° pulses in the D editing stage, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multidimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multidimensional nuclear magnetic resonance diffusion coefficient D editing stage of water, and the signal attenuation of the multidimensional nuclear magnetic resonance T2 acquisition stage of water in the rock pores by the numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until N groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where N is an integer greater than 1; N groups of nuclear magnetic resonance echo data of rocks saturated with oil, water and CO2 are used as the multi-dimensional nuclear magnetic resonance echo data.
3. The method for acquiring multidimensional nuclear magnetic resonance spectrum according to claim 1, wherein When the multi-dimensional nuclear magnetic resonance is a T1-T2 two-dimensional nuclear magnetic resonance, the acquisition process of the multi-dimensional nuclear magnetic resonance is simulated by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multi-dimensional nuclear magnetic resonance echo data, including: Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of water in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water to obtain multiple echo signal intensities of water in the T2 acquisition phase; Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing phase of the oil in the rock pores and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil to obtain multiple echo signal intensities of the oil in the T2 acquisition phase; The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks; Changing the waiting time of the T1 editing stage, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal recovery of the multi-dimensional nuclear magnetic resonance longitudinal relaxation time T1 editing stage of water in the rock pores, the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition stage of water, and obtaining multiple echo signal intensities of water in the T2 acquisition stage, until M groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where M is an integer greater than 1; M groups of oil-, water- and CO2-saturated rock nuclear magnetic resonance echo data are used as the multidimensional nuclear magnetic resonance echo data.
4. The method for acquiring multidimensional nuclear magnetic resonance spectrum according to claim 1, wherein When the multi-dimensional nuclear magnetic resonance is a D-T2 two-dimensional nuclear magnetic resonance, the acquisition process of the multi-dimensional nuclear magnetic resonance is simulated by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multi-dimensional nuclear magnetic resonance echo data, including: According to the three-dimensional digital core model of the rock, the multi-dimensional nuclear magnetic resonance diffusion coefficient D of water in the rock pores is sequentially simulated by the numerical simulation method. The signal attenuation during the editing phase and the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water are simulated to obtain multiple echo signal intensities of water during the T2 acquisition phase. Based on the three-dimensional digital core model of the rock, the numerical simulation method is used to sequentially simulate the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient D editing phase of the oil and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of the oil in the rock pores to obtain multiple echo signal intensities of the oil in the T2 acquisition phase; The intensity of each echo signal of water in the T2 acquisition phase and the intensity of each echo signal of oil in the T2 acquisition phase are added together to obtain a set of nuclear magnetic resonance echo data of saturated oil, water and CO2 rocks; changing the time interval between the two 180° pulses in the D editing phase, returning to the step of simulating, based on the three-dimensional digital core model of the rock, the signal attenuation of the multi-dimensional nuclear magnetic resonance diffusion coefficient of water in the rock pores in the D editing phase and the signal attenuation of the multi-dimensional nuclear magnetic resonance T2 acquisition phase of water by the numerical simulation method, and obtaining multiple echo signal intensities of water in the T2 acquisition phase, until P groups of nuclear magnetic resonance echo data of saturated oil, water, and CO2 rocks are obtained, where P is an integer greater than 1; The P group of oil-, water- and CO2-saturated rock nuclear magnetic resonance echo data is used as the multi-dimensional nuclear magnetic resonance echo data.
5. The method for acquiring multidimensional nuclear magnetic resonance spectrum according to claim 1, wherein: After constructing a three-dimensional digital core according to the rock CT image, the method further includes: Based on the three-dimensional digital core, a three-dimensional digital core model of rock containing only saturated oil and water is established using the morphological method.
6. The method for acquiring multidimensional nuclear magnetic resonance spectrum according to claim 5, wherein: Based on the three-dimensional digital core, a three-dimensional digital core model of rock containing only saturated oil and water is established using the morphological method, including: Assuming that the pores of the three-dimensional digital core contain only oil and water wetting phase fluids, the pore space is subjected to the opening operation in the morphological method to simulate the non-wetting phase fluids containing oil and water, thereby obtaining a three-dimensional digital core model of the rock containing only oil and water.
7. A multi-dimensional nuclear magnetic resonance spectrum acquisition device, characterized in that: include: A construction module is used to obtain a rock CT image and construct a three-dimensional digital rock core based on the rock CT image; The first establishment module is used to establish a three-dimensional digital core model of rock saturated with oil, water and CO2 under different saturation conditions using a morphological method based on the three-dimensional digital core; The first establishment module includes: a first simulation unit for treating saturated oil and water as the same-phase fluids, and using the opening operation in the morphological method to simulate the distribution state of CO2 and the oil-water mixture in the three-dimensional digital core to obtain an intermediate model; a second simulation unit for treating CO2 and the skeleton as the same-phase medium in the intermediate model, and using the opening operation in the morphological method to simulate the distribution state of oil and water to obtain the distribution state of the three-phase fluid of CO2, oil and water, and replacing the three-phase fluids of CO2, oil and water with different numbers to obtain the three-dimensional digital core model of the rock; a simulation module for simulating a multidimensional nuclear magnetic resonance acquisition process by a numerical simulation method based on the three-dimensional digital core model of the rock to obtain multidimensional nuclear magnetic resonance echo data; the multidimensional nuclear magnetic resonance includes T1-D-T2 three-dimensional nuclear magnetic resonance, T1-T2 two-dimensional nuclear magnetic resonance, or D-T2 two-dimensional nuclear magnetic resonance; The inversion module is used to invert the multidimensional nuclear magnetic resonance echo data to obtain a corresponding multidimensional nuclear magnetic resonance spectrum.
8. A multi-dimensional nuclear magnetic resonance spectrum acquisition device, characterized in that: include: memory for storing computer programs; A processor, configured to implement the steps of the method for acquiring a multi-dimensional nuclear magnetic resonance spectrum according to any one of claims 1 to 6 when executing the computer program.
9. A readable storage medium, characterized in that: The readable storage medium stores a computer program, and when the computer program is executed by a processor, the steps of the method for acquiring a multi-dimensional nuclear magnetic resonance spectrum according to any one of claims 1 to 6 are implemented.