Method, device, equipment and medium for predicting imbibition utilization efficiency of blocky shale oil

Abstract

The invention relates to the technical field of unconventional shale oil and gas reservoir exploration, in particular to a method, a device, equipment and a medium for predicting the imbibition utilization efficiency of blocky shale oil. According to the method, under the condition that the rock core is not damaged, the content of the oil in different aperture intervals after the rock core is pressurized and saturated with the oil is rapidly quantified through the two-dimensional nuclear magnetic resonance technology, and then the content of the oil in the different aperture intervals of the block-shaped shale under the stratum condition is calculated by utilizing the respective imbibition and oil discharge efficiency functions and the overall critical use lower limit of the different aperture intervals of the block-shaped shale. The imbibition and oil discharge efficiency of the whole shale core can be rapidly, accurately and quantitatively represented. By means of the method, the shale oil core imbibition utilization efficiency can be rapidly and accurately measured under the condition that the core is not damaged, and effective means are provided for revealing the interaction reaction mechanism of shale reservoir rock and injected media, the shale reservoir rock imbibition oil extraction mechanism, the shale reservoir recovery efficiency improving mechanism and the like. The method disclosed by the invention has important engineering value and scientific significance for promoting effective economic development of shale reservoirs.

Description

Technical Field

[0001] This disclosure relates to the field of unconventional shale oil and gas reservoir exploration technology, and in particular to methods, apparatus, equipment and media for predicting the permeation and utilization efficiency of massive shale oil. Background Technology

[0002] Shale oil and gas resources are widely distributed and have enormous potential, making them an important strategic alternative energy source. However, achieving large-scale, efficient development remains the biggest challenge. Shale oil reservoirs exhibit multi-scale pore-throat structures, with widespread micro- and nano-scale pores and throats. Compared to conventional reservoirs, they possess higher capillary forces and more pronounced adsorption-driven oil recovery. Adsorption-driven oil recovery refers to the process where, when fractures surrounding the matrix rock block are filled with water, the injected water from the fractures or large pores is drawn into the matrix rock block by capillary forces, displacing the crude oil in the matrix and thus achieving a redistribution of oil and water outside the matrix. Adsorption occurs in several key stages of shale oil reservoir development, including volumetric fracturing and water injection, and is one of the important factors affecting shale oil production capacity. Clarifying the characteristics of adsorption and evaluating its utilization efficiency are of great significance for improving shale oil recovery rates.

[0003] Shale-type (massive) lithology is mudstone and shale, with poor bedding. Currently, the evaluation of shale oil adsorption and mobilization efficiency mainly uses laboratory tests for self-adsorption and discharge, which are characterized by three methods: mass method, volumetric method, and nuclear magnetic resonance (NMR) method. The mass method involves processing shale cores and placing them in a pre-prepared adsorption solution. During adsorption, the core mass changes, and the change in mass over time is monitored using an electronic balance. The volumetric method uses a self-made water absorption instrument. The processed core is placed in a pre-prepared adsorption solution, and the volume of oil discharged is recorded over time. The NMR method uses the T2 spectrum obtained during the adsorption experiment to distinguish the fluid mobilization in different pore sizes. The mass and volumetric methods involve relatively cumbersome manual measurement and are prone to significant errors. The NMR method offers precise measurement, but a simple one-dimensional NMR cannot distinguish between oil and water signals, requiring the addition of heavy water or soaking in strong water to eliminate the water signal. In addition, these methods require a relatively long period of percolation experiments (more than one month) to obtain percolation mobilization efficiency. Furthermore, the core samples after chemical percolation experiments are contaminated and damaged to a certain extent, making it impossible to conduct other experiments. Therefore, this disclosure proposes a novel method for rapidly predicting the percolation mobilization efficiency of massive shale oil using two-dimensional nuclear magnetic resonance technology, which addresses the shortcomings of existing shale oil percolation mobilization efficiency evaluation methods.

[0004] In summary, there is an urgent need for a technical solution to rapidly predict the percolation and utilization efficiency of massive shale oil. Summary of the Invention

[0005] To address the aforementioned issues, this disclosure provides a method, apparatus, equipment, and medium for predicting the adsorption and mobilization efficiency of massive shale oil, thereby resolving the problems of excessively long experimental cycles and core contamination damage in existing shale oil adsorption and mobilization efficiency evaluation methods.

[0006] In a first aspect, a method for predicting the adsorption and utilization efficiency of massive shale oil, the method comprising:

[0007] Take several core samples of each physical property grade of blocky shale oil;

[0008] Two-dimensional nuclear magnetic resonance (NMR) spectra were generated for each core sample;

[0009] The two-dimensional NMR spectrum oil signal was divided into bound oil, mesoporous oil region and kerogen region, respectively.

[0010] The first permeation rate was calculated based on the two-dimensional NMR spectrum; the second permeation rate was calculated for the kerogen region.

[0011] Fresh samples of massive shale from the target block were taken, and two-dimensional nuclear magnetic resonance spectra were measured to obtain the saturated oil signal intensity of the fresh samples. The total saturated oil signal intensity of the fresh samples was determined based on the saturated oil signal intensity of the fresh samples.

[0012] Based on the saturated oil signal quantity and the first permeation rate of the fresh sample, the oil signal quantity after permeation in the bound oil and mesoporous oil regions of the fresh sample is determined; based on the saturated oil signal quantity and the second permeation rate of the fresh sample, the oil signal quantity after permeation in the kerogen region of the fresh sample is determined.

[0013] The oil adsorption and utilization efficiency of the target block's massive shale oil is calculated based on the total saturated oil signal quantity of the fresh sample, the oil signal quantity after infiltration in the bound oil and mesopore oil zones of the fresh sample, and the oil signal quantity after infiltration in the kerogen zone of the fresh sample.

[0014] Furthermore, several core samples were taken from each physical property level of the massive shale oil, including:

[0015] Well logging facies technology was used to identify shale properties at different locations in the reference wellbore, and wire cutting technology was used to drill several core samples of the same size for each property level.

[0016] Furthermore, two-dimensional nuclear magnetic resonance (NMR) spectra were generated for each core sample, including:

[0017] Two-dimensional nuclear magnetic resonance (NMR) was used to measure the T1 and T2 spectra after oil saturation for each core sample to establish a two-dimensional NMR spectrum after oil saturation, and the T1 and T2 spectra after infiltration were measured to establish a two-dimensional NMR spectrum after infiltration.

[0018] Furthermore, the T1 and T2 spectra of the saturated oil were measured to establish a two-dimensional NMR spectrum of the saturated oil, including:

[0019] Each core was placed in a constant temperature chamber and dried for H1 hours to remove residual oil. The core was then placed in a piston container, evacuated for H2 hours to saturate with simulated shale oil, and aged in a constant temperature chamber for H3 hours. The core was then removed and weighed. Two-dimensional nuclear magnetic resonance technology was used to determine the T1 and T2 spectra after oil saturation without damaging the core.

[0020] Furthermore, the determination of the T1 and T2 spectra after saturation of the oil to establish a two-dimensional NMR spectrum after saturation also includes:

[0021] Software quantifies the amount of oil signal W in different pore sizes of each core under saturated oil conditions. xy，饱和 , where x is the T2 spectral coordinate value and y is the T1 spectral coordinate value.

[0022] Furthermore, the T1 and T2 spectra after osmosis were measured to establish a two-dimensional NMR spectrum after osmosis, including:

[0023] After each core was saturated with oil, it was placed in a piston container, and an adsorption medium was added. The piston container was then placed in a constant temperature chamber set to the formation temperature and pressurized to the reservoir pressure to conduct an adsorption experiment. After day 1 of the experiment, the core was removed and weighed, and the T1 and T2 spectra after adsorption were tested using two-dimensional nuclear magnetic resonance technology.

[0024] Furthermore, the determination of post-osmosis T1 and T2 spectra to establish a post-osmosis two-dimensional NMR spectrum also includes:

[0025] Software quantifies the amount of oil signal W in different pore sizes after permeation of each core sample. xy，渗吸后 , where x is the coordinate value of T2 and y is the coordinate value of T1.

[0026] Furthermore, the two-dimensional NMR spectrum oil signal was divided into bound oil, mesoporous oil region, and kerogen region, including:

[0027] The oil signals in the two-dimensional NMR spectra after saturation and after percolation were divided into bound oil, mesoporous oil, and kerogen regions, respectively.

[0028] Furthermore, the oil signals in the two-dimensional NMR spectra after saturation and after percolation were divided into bound oil, mesoporous oil, and kerogen regions, respectively, including:

[0029] The oil content characteristics of the kerogen region differ from those of bound oil and mesoporous oil. The bound oil, mesoporous oil and kerogen regions can be distinguished by T1, T2 and T1 / T2 ratio.

[0030] Furthermore, the first permeation rate is calculated based on the two-dimensional NMR spectrum, including:

[0031] Based on the two-dimensional NMR spectrum of the oil signal, the permeation rate I was calculated. xy The formula is as follows:

[0032] Ixy =(W xy，饱和 -W xy，渗吸后 ) / W xy，饱和

[0033] Among them, W xy，饱和 W represents the oil signal quantity after saturation. xy，渗吸后 This represents the oil signal quantity after absorption, where x is the T2 coordinate value and y is the T1 coordinate value.

[0034] Furthermore, the calculation of the first permeation rate based on the two-dimensional NMR spectrum also includes:

[0035] For each core sample, the average of the permeability values ​​at the same T1 and T2 coordinates is taken as the first permeability, as shown in the following formula:

[0036]

[0037] Among them, I xy,最终 This represents the final absorption rate at points T2 = x and T1 = y; I xy，n Let T2 = x and T1 = y represent the permeation rate at point T2 = x and T1 = y in the nth core sample; n represents the total number of samples.

[0038] Furthermore, the second osmotic rate is calculated for the kerogen zone, including:

[0039] Based on two-dimensional NMR spectroscopy of oil signals, software quantization is used to determine the total oil signal quantity W in the kerogen region under saturated oil conditions. 干酪根，饱和 And the total oil signal W in the kerogen region after infiltration 干酪根，渗吸后 The osmotic mobilization efficiency I in the kerogen region is calculated using the following formula. xy，干酪根 :

[0040]

[0041] Furthermore, the calculation of the second osmotic rate for the kerogen zone also includes:

[0042] The average permeation mobilization efficiency of the kerogen zone in each core is taken as the second permeation rate, as shown in the following formula:

[0043]

[0044] Among them, I 干酪根,最终 This indicates the final permeation rate value of the kerogen display zone; I 干酪根,n This represents the percolation rate of the kerogen region in the nth sample; n represents the total number of samples.

[0045] Furthermore, fresh samples of massive shale from the target block were taken, and two-dimensional nuclear magnetic resonance spectra were measured to obtain the saturated oil signal quantity of the fresh samples, including:

[0046] Fresh samples of massive shale from the target block were taken, and T1 and T2 spectra after saturation with oil were measured to establish a two-dimensional NMR spectrum of the fresh samples after saturation with oil.

[0047] The signal intensity of the saturated oil in the fresh sample was determined based on the two-dimensional NMR spectrum of the saturated oil.

[0048] Furthermore, fresh samples of massive shale from the target block were taken, and T1 and T2 spectra after saturation with oil were measured to establish a two-dimensional NMR spectrum of the fresh samples after saturation with oil, including:

[0049] Fresh samples of blocky shale from the target area were obtained by drilling using wire cutting technology. The T1 and T2 spectra of the fresh samples were measured using two-dimensional nuclear magnetic resonance technology without damaging the core.

[0050] Furthermore, the total saturated oil signal quantity of the fresh sample is determined based on the saturated oil signal quantity of the fresh sample, including:

[0051] Software quantifies the amount of oil signal W in different pore sizes of each fresh core sample. xy，新鲜 Then, the oil signal data of all fresh sample oil signal points are summed to obtain the total saturated oil signal W of the fresh sample. 总和，新鲜 , where x is the coordinate value of T2 and y is the coordinate value of T1.

[0052] Furthermore, based on the saturated oil signal intensity and the first permeation rate of the fresh sample, the oil signal intensity after permeation into the bound oil and mesopore oil regions of the fresh sample is determined, including:

[0053] Calculate the oil signal quantity W after seepage at different T1 and T2 coordinates. xy，渗吸后 As shown in the following formula:

[0054] W xy，渗吸后 =W xy，新鲜 -W xy，新鲜 ×I xy,最终

[0055] Among them, I xy,最终 W represents the final absorption rate at points T2 = x and T1 = y. xy，新鲜 W represents the oil signal quantity at points T2 = x and T1 = y for fresh samples; xy，渗吸后 This represents the oil signal quantity at points T2=x and T1=y after percolation.

[0056] W, the bound oil and the mesopore oil zone xy，渗吸后 The oil signal data after all percolation in the data volume are summed to obtain the oil signal W of the bound oil and mesopore oil regions of the fresh sample. 束缚油、中孔油，渗吸后 .

[0057] Furthermore, based on the saturated oil signal intensity and the second percolation rate of the fresh sample, the oil signal intensity after percolation in the kerogen region of the fresh sample is determined, including:

[0058] The oil signal data from the kerogen region of the fresh sample were summed to obtain W. 干酪根，新鲜 The oil signal quantity W after osmosis in the kerogen zone is calculated using the following formula. 干酪根，渗吸后 :

[0059] W 干酪根，渗吸后＝ W 干酪根，新鲜 -W 干酪根，新鲜 ×I 干酪根,最终

[0060] Among them, I 干酪根,最终 W represents the final permeation rate value of the kerogen display zone. 干酪根，新鲜 This represents the total oil signal quantity in the kerogen region of a fresh sample.

[0061] Furthermore, based on the total saturated oil signal quantity of fresh samples, the oil signal quantity after permeation in the bound oil and mesoporous oil zones of fresh samples, and the oil signal quantity after permeation in the kerogen zone of fresh samples, the permeation and utilization efficiency of massive shale oil in the target block is calculated, including:

[0062] The formula for calculating the adsorption mobilization efficiency I of fresh shale oil core samples is as follows:

[0063]

[0064] Among them, W 总和，新鲜 W represents the sum of the signal quantities of the fresh sample oil; 干酪根，渗吸后 W represents the total oil signal after permeation in the kerogen oil region; 束缚油、中孔油，渗吸后 This represents the sum of the oil signal quantities after the absorption of bound oil and mesopore oil; I represents the absorption efficiency.

[0065] Secondly, an apparatus for predicting the permeation and utilization efficiency of massive shale oil includes: a control sample permeation rate measurement unit and a fresh sample permeation and utilization efficiency measurement unit.

[0066] The reference sample permeability measurement unit is used to take several core samples of each physical property grade of blocky shale oil.

[0067] The reference sample permeation rate measurement unit is also used to establish a two-dimensional nuclear magnetic resonance spectrum for each core sample;

[0068] The reference sample permeation rate measurement unit is also used to divide the two-dimensional NMR spectrum oil signal into bound oil, mesoporous oil region and kerogen region, respectively.

[0069] The control sample permeation rate measurement unit is also used to calculate the first permeation rate based on the two-dimensional NMR spectrum; and to calculate the second permeation rate for the kerogen region;

[0070] The fresh sample infiltration and kinetic efficiency measurement unit is used to take fresh samples of blocky shale from the target area, measure the two-dimensional nuclear magnetic resonance spectrum, and obtain the saturated oil signal of the fresh sample; the total saturated oil signal of the fresh sample is determined based on the saturated oil signal of the fresh sample.

[0071] The fresh sample percolation efficiency measurement unit is also used to determine the oil signal after percolation in the bound oil and mesoporous oil regions of the fresh sample based on the saturated oil signal and the first percolation rate; and to determine the oil signal after percolation in the kerogen region of the fresh sample based on the saturated oil signal and the second percolation rate.

[0072] The fresh sample percolation efficiency measurement unit is also used to calculate the percolation efficiency of the target block of massive shale oil based on the total saturated oil signal of the fresh sample, the oil signal after percolation of the bound oil and mesopore oil zones of the fresh sample, and the oil signal after percolation of the kerogen zone of the fresh sample.

[0073] Thirdly, an electronic device includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus;

[0074] Memory, which stores computer programs;

[0075] When a processor executes a computer program stored in memory, it implements the above-described method for predicting the percolation and mobilization efficiency of massive shale oil.

[0076] Fourthly, a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method for predicting the percolation and mobilization efficiency of massive shale oil.

[0077] This disclosure has at least the following beneficial effects:

[0078] This disclosure rapidly quantifies the oil content in different pore size ranges of a pressurized and saturated shale core (fresh sample) without damaging the core under two-dimensional nuclear magnetic resonance (NMR) technology. Furthermore, by utilizing the adsorption and drainage efficiency functions for different pore size ranges of massive shale under formation conditions established in this disclosure, as well as the overall critical utilization lower limit, the adsorption and drainage efficiency of the entire shale core can be rapidly and accurately quantitatively characterized. This method enables the rapid, accurate, and non-destructive determination of shale oil core adsorption and utilization efficiency. It provides an effective means to elucidate the interaction mechanism between shale reservoir rocks and the injection medium, the adsorption and recovery mechanism of shale reservoir rocks, and the enhanced oil recovery mechanism of shale reservoirs. This method has significant engineering value and scientific significance for promoting the effective and economical development of shale reservoirs.

[0079] Other features and advantages of this disclosure will be set forth in the following description and will be apparent in part from the description or may be learned by practicing the disclosure. The objects and other advantages of this disclosure may be realized and obtained by means of the structures pointed out in the description and the accompanying drawings. Attached Figure Description

[0080] To more clearly illustrate the technical solutions in the embodiments of this disclosure or the prior art, the drawings used in the description of the embodiments or the prior art will be briefly introduced below. Obviously, the drawings described below are some embodiments of this disclosure. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0081] Figure 1 This is a schematic diagram of the prediction method flow according to an embodiment of the present disclosure;

[0082] Figure 2 A schematic diagram showing the well logging facies identification results of a representative well in a block of massive shale reservoirs.

[0083] Figure 3 A schematic diagram showing the results of two-dimensional nuclear magnetic resonance quantification of oil content in different pore size ranges for pressurized saturated oil in a certain rock core;

[0084] Figure 4 A schematic diagram showing the results of two-dimensional nuclear magnetic resonance quantification of the remaining oil content in different pore size ranges after permeation of a certain rock core;

[0085] Figure 5 This is a schematic diagram of the predictive device structure according to an embodiment of the present disclosure;

[0086] Figure 6 This is a schematic diagram of the electronic device structure. Detailed Implementation

[0087] To make the objectives, technical solutions, and advantages of the embodiments of this disclosure clearer, the technical solutions of the embodiments of this disclosure will be clearly and completely described below with reference to the accompanying drawings. Obviously, the described embodiments are only some embodiments of this disclosure, and not all embodiments. Based on the embodiments of this disclosure, all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of this disclosure.

[0088] like Figure 1 As shown, a method for predicting the adsorption and mobilization efficiency of massive shale oil is provided, the method comprising:

[0089] S101, take several cores of each physical property grade of blocky shale oil;

[0090] S102, establish two-dimensional nuclear magnetic resonance spectra for each core sample;

[0091] S103, the two-dimensional NMR spectrum oil signal is divided into bound oil, mesoporous oil region and kerogen region respectively;

[0092] S104, calculate the first permeation rate based on the two-dimensional NMR spectrum; calculate the second permeation rate for the kerogen region;

[0093] S105. Take fresh samples of massive shale from the target block, measure the two-dimensional nuclear magnetic resonance spectrum, and obtain the saturated oil signal of the fresh sample; determine the total saturated oil signal of the fresh sample based on the saturated oil signal of the fresh sample.

[0094] S106, Based on the saturated oil signal quantity and the first permeation rate of the fresh sample, determine the oil signal quantity after permeation of the bound oil and mesoporous oil regions of the fresh sample; based on the saturated oil signal quantity and the second permeation rate of the fresh sample, determine the oil signal quantity after permeation of the kerogen region of the fresh sample.

[0095] S107. Based on the total saturated oil signal of fresh samples, the oil signal after permeation in the bound oil and mesopore oil zones of fresh samples, and the oil signal after permeation in the kerogen zone of fresh samples, the permeation utilization efficiency of the massive shale oil in the target block is calculated.

[0096] The specific implementation details are as follows:

[0097] I. Rock Sample Preparation

[0098] Step 1: The shale-type (massive) lithology is mudstone and shale with poor bedding. Based on the physical properties of the massive shale reservoir in a target block, well logging facies technology is used to identify the shale physical properties at different locations along the wellbore (e.g., Figure 2 As shown, several core samples (core diameter 2.5cm, core length approximately 8cm) of each physical property level were drilled using wire cutting technology. The porosity and permeability of the core samples under overburden conditions were accurately measured using a high-precision overburden porosimeter.

[0099] II. Conducting experiments combining the percolating agent percolation with two-dimensional nuclear magnetic resonance imaging.

[0100] Step 2: Each core sample was placed in a constant temperature chamber and dried at 110℃ for 24 hours to remove residual oil, effectively "washing" the oil. The core was then placed in a piston container and vacuum-sealed for 48 hours until saturated with simulated shale oil (5MPa injection pressure). After aging in a constant temperature chamber for 48 hours, the core was removed and weighed. Two-dimensional nuclear magnetic resonance (NMR) technology was used to determine the T1 and T2 spectra after oil saturation without damaging the core (e.g., ...). Figure 3 (As shown). The software quantifies the oil signal quantity W in different pore sizes (x is the T2 coordinate value, y is the T1 coordinate value) of each core under saturated oil conditions. xy，饱和 .

[0101] Step 3: Place each core sample into a piston container, add the permeation medium, place the piston container in a constant temperature chamber set to the formation temperature and pressurized to the reservoir pressure to conduct the permeation experiment. After 30 days of experimentation, remove the core samples, weigh them, and use two-dimensional nuclear magnetic resonance (NMR) technology to test the T1 and T2 spectra after permeation (e.g., ...). Figure 4 (As shown). The software quantifies the oil signal quantity W in different pore sizes of each core sample after permeation (x is the T2 coordinate value, y is the T1 coordinate value). xy，渗吸后 .

[0102] III. Experimental Data Processing and Data Volume Establishment

[0103] The two-dimensional nuclear magnetic resonance (NMR) spectra of shale-type (massive) shale oil show oil signal regions in bound oil and mesoporous oil areas, and kerogen areas (such as...). Figure 3 As shown, the oil content characteristics of the kerogen display area differ significantly from those of bound oil and mesoporous oil. The movable oil in this area is adsorbed oil. The adsorption amount is positively correlated with the total amount of kerogen; therefore, its permeation rate should be considered based on the overall signal volume of the kerogen. Consequently, the data processing mode for the kerogen oil region differs significantly from that for bound oil and mesoporous oil regions, and they should be processed separately. Therefore, two data processing modes are required. Correspondingly, two functional relationships between pore size and permeation efficiency need to be established for the oil display area.

[0104] 1) Establishment of a two-dimensional data volume for single-sample percolation rate:

[0105] A. Regarding the bound oil ③ and the medium-hole oil zone ④:

[0106] Step 4: Apply I xy =(W xy，饱和 -W xy，渗吸后 ) / W xy，饱和 The basic formula is used to calculate the permeation rate (I) at each oil signal location (different T1, T2 coordinates) in the two-dimensional nuclear magnetic resonance image for the bound oil ③ and mesoporous oil region ④. xy This generates a two-dimensional data volume of permeability for each core sample at different T1 and T2 coordinates.

[0107] B. Regarding the kerogen zone ⑤:

[0108] Step 5: Quantify the total oil signal W in the kerogen region (⑤) under saturated oil conditions using nuclear magnetic resonance data processing software. 干酪根，饱和 And the total oil signal W in the kerogen region (⑤) after infiltration. 干酪根，渗吸后 The osmotic kinetic efficiency I in the kerogen region is calculated using formula (1). xy，干酪根 .

[0109]

[0110] 2) Establishment of a two-dimensional data volume function model for the permeability and absorption efficiency of block-type (massive) shale oil:

[0111] Considering the different pore size ranges of shale cores with varying physical properties, a two-dimensional data volume of permeability was constructed based on different T1 and T2 coordinates of each shale reservoir core. This data was then integrated to form a two-dimensional data volume function model of permeability efficiency for block-type (massive) shale oil.

[0112] A. Regarding the bound oil ③ and the medium-hole oil zone ④:

[0113] Step 6: Using data processing software, average the permeability values ​​at the same T1 and T2 coordinates. Permeability values ​​at different T1 and T2 coordinates are directly merged into the final two-dimensional permeability data volume. This yields a two-dimensional data volume of the permeability efficiency of the shale-type (massive) shale oil bound oil ③ and the mesoporous oil zone ④ in this block. As shown in formula (2):

[0114]

[0115] Among them, I xy,最终 —The final absorption rate at point T2 = x, T1 = y; x n —The nth value; n—The total number of valid data samples.

[0116] B. Regarding the kerogen zone ⑤:

[0117] Step 7: Calculate the average value of the permeation efficiency of the kerogen zone of each core according to formula (3), and obtain a two-dimensional data volume function model of the permeation efficiency of the shale oil kerogen zone of this block.

[0118]

[0119] Among them, I 干酪根,最终 —The final permeation rate value of the kerogen display zone; I 干酪根,n —The percolation rate of the kerogen region in the nth sample; n—The total number of valid samples.

[0120] As the number of actual experimental samples increases, the two-dimensional data volume of percolation efficiency can be continuously updated.

[0121] IV. Determining the absorption and utilization efficiency of fresh shale (massive) oil.

[0122] Step 8: Drill fresh samples of massive shale from the target block using wire cutting technology. Measure the T1 and T2 spectra of the fresh samples using two-dimensional nuclear magnetic resonance (NMR) without damaging the core. Quantify the oil signal quantity W in different pore sizes (x = T2 coordinate, y = T1) of each fresh core sample using software. xy，新鲜 Then, the oil signal data of all fresh sample oil signal points are summed to obtain W. 总和，新鲜 .

[0123] A. Regarding the bound oil ③ and the medium-hole oil zone ④:

[0124] Based on the two-dimensional data volume of the shale-type (massive) shale oil bound oil and the permeability of the mesoporous oil zone in this block, the oil signal quantity W after permeation at different T1 and T2 coordinates was calculated. xy，渗吸后 As shown in formula (4):

[0125] W xy，渗吸后 =W xy，新鲜 -W xy，新鲜 ×I xy,最终 ................(4)

[0126] Among them, I xy,最终 —The final absorption rate at point T2 = x, T1 = y; W xy，新鲜 —Oil signal quantity at points T2=x, T1=y for fresh samples; W xy，渗吸后 —Oil signal quantity at points T2=x and T1=y after percolation.

[0127] W, the bound oil and the mesopore oil zone xy，渗吸后 Sum all the oil signal data after absorption in the data body to obtain W. 束缚油、中孔油，渗吸后 .

[0128] B. Regarding the kerogen zone ⑤:

[0129] Step 9: Sum the oil signal data of the kerogen region of the fresh sample to obtain W. 干酪根，新鲜 Based on the two-dimensional data volume function model of the kerogen zone permeation efficiency of shale oil in this area, the oil signal quantity W after permeation in the kerogen zone is calculated using formula (5). 干酪根，渗吸后 .

[0130] W 干酪根，渗吸后＝ W 干酪根，新鲜 -W 干酪根，新鲜 ×I 干酪根,最终 ...............(5)

[0131] Among them, I 干酪根,最终 —The final permeation rate value of the kerogen display zone, W 干酪根，新鲜 —The sum of oil signal quantities in the kerogen region of the fresh sample.

[0132] Step 10: Obtain W through the above method 总和，新鲜 and W 总和，渗吸后 The data was used to calculate the permeation efficiency I of the shale oil core, as shown in formula (3):

[0133]

[0134] Among them, W 总和，新鲜 —Sum of signal quantities in fresh sample oil; W 干酪根，渗吸后—The sum of oil signal quantities after permeation and adsorption in the kerogen oilfield of shale oil cores; W 束缚油、中孔油，渗吸后 —The sum of the oil signal after the core's bound oil and medium-pore oil zones have been absorbed; I—The core's absorption utilization efficiency.

[0135] Step 11: Based on the established two-dimensional data volume (different T1, T2) of the shale oil permeation efficiency under the target block formation conditions, a method for predicting the permeation efficiency of massive shale under formation conditions using two-dimensional nuclear magnetic resonance technology was developed without damaging the core. Corresponding calculation software was independently developed, thus realizing the prediction of permeation efficiency under massive shale formation conditions without damaging the core.

[0136] like Figure 5 As shown, an apparatus for predicting the permeation efficiency of massive shale oil includes: a control sample permeation rate measuring unit 501 and a fresh sample permeation efficiency measuring unit 502.

[0137] The reference sample permeability measurement unit 501 is used to take several core samples of each physical property grade of blocky shale oil.

[0138] The reference sample permeation rate measurement unit 501 is also used to establish a two-dimensional nuclear magnetic spectrum for each core sample;

[0139] The reference sample permeation rate measurement unit 501 is also used to divide the two-dimensional NMR spectrum oil signal into bound oil, mesoporous oil region and kerogen region, respectively.

[0140] The control sample permeation rate measuring unit 501 is also used to calculate the first permeation rate based on the two-dimensional NMR spectrum; and to calculate the second permeation rate for the kerogen region;

[0141] The fresh sample infiltration and kinetic efficiency measurement unit 502 is used to take fresh samples of blocky shale from the target area, measure the two-dimensional nuclear magnetic resonance spectrum, and obtain the saturated oil signal quantity of the fresh sample; the total saturated oil signal quantity of the fresh sample is determined based on the saturated oil signal quantity of the fresh sample.

[0142] The fresh sample percolation efficiency measurement unit 502 is also used to determine the oil signal quantity after percolation of the bound oil and mesoporous oil regions of the fresh sample based on the saturated oil signal quantity and the first percolation rate; and to determine the oil signal quantity after percolation of the kerogen region of the fresh sample based on the saturated oil signal quantity and the second percolation rate.

[0143] The fresh sample infiltration efficiency measurement unit 502 is also used to calculate the infiltration efficiency of the target block of massive shale oil based on the total saturated oil signal of the fresh sample, the oil signal after infiltration of the bound oil and mesopore oil zones of the fresh sample, and the oil signal after infiltration of the kerogen zone of the fresh sample.

[0144] like Figure 6As shown, this disclosure provides an electronic device, including a processor 601, a communication interface 602, a memory 603, and a communication bus 604, wherein the processor 601, the communication interface 602, and the memory 603 communicate with each other through the communication bus 604;

[0145] Memory 603 stores computer programs;

[0146] The processor 601 implements the above method when executing a computer program stored in the memory 603.

[0147] This disclosure provides a computer-readable storage medium storing a computer program that, when executed by a processor, implements the above-described method.

[0148] The computer-readable storage medium may be included in the device / apparatus described in the above embodiments; or it may exist independently and not assembled into the device / apparatus. The computer-readable storage medium carries one or more programs that, when executed, implement the method according to the embodiments of this disclosure.

[0149] According to embodiments of this disclosure, the computer-readable storage medium can be a non-volatile computer-readable storage medium, such as including, but not limited to: portable computer disks, hard disks, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM or flash memory), portable compact disk read-only memory (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination thereof. In this disclosure, the computer-readable storage medium can be any tangible medium that contains or stores a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0150] Although the present disclosure has been described in detail with reference to the foregoing embodiments, those skilled in the art should understand that modifications can still be made to the technical solutions described in the foregoing embodiments, or equivalent substitutions can be made to some of the technical features; and such modifications or substitutions do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of the embodiments of the present disclosure.

Claims

1. A method for predicting the adsorption and utilization efficiency of massive shale oil, characterized in that, The method includes: Take several core samples of each physical property grade of blocky shale oil; Two-dimensional nuclear magnetic resonance (NMR) spectra were generated for each core sample; The two-dimensional NMR spectrum oil signal was divided into bound oil, mesoporous oil region and kerogen region, respectively. The first permeation rate was calculated based on the two-dimensional NMR spectrum; the second permeation rate was calculated for the kerogen region. Fresh samples of massive shale from the target block were taken, and two-dimensional nuclear magnetic resonance spectra were measured to obtain the saturated oil signal intensity of the fresh samples. The total saturated oil signal intensity of the fresh samples was determined based on the saturated oil signal intensity of the fresh samples. Based on the saturated oil signal quantity and the first permeation rate of the fresh sample, the oil signal quantity after permeation in the bound oil and mesoporous oil regions of the fresh sample is determined; based on the saturated oil signal quantity and the second permeation rate of the fresh sample, the oil signal quantity after permeation in the kerogen region of the fresh sample is determined. The oil adsorption and utilization efficiency of the target block's massive shale oil is calculated based on the total saturated oil signal quantity of the fresh sample, the oil signal quantity after infiltration in the bound oil and mesopore oil zones of the fresh sample, and the oil signal quantity after infiltration in the kerogen zone of the fresh sample.

2. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Take several core samples from each physical property grade of massive shale oil, including: Well logging facies technology was used to identify shale properties at different locations in the reference wellbore, and wire cutting technology was used to drill several core samples of the same size for each property level.

3. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Two-dimensional nuclear magnetic resonance (NMR) spectra were generated for each core sample, including: Two-dimensional nuclear magnetic resonance (NMR) was used to measure the T1 and T2 spectra after oil saturation for each core sample to establish a two-dimensional NMR spectrum after oil saturation, and the T1 and T2 spectra after infiltration were measured to establish a two-dimensional NMR spectrum after infiltration.

4. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 3, characterized in that, The T1 and T2 spectra of saturated oil were measured to establish a two-dimensional NMR spectrum of saturated oil, including: Each core was placed in a constant temperature chamber and dried for H1 hours to remove residual oil. The core was then placed in a piston container, evacuated for H2 hours to saturate with simulated shale oil, and aged in a constant temperature chamber for H3 hours. The core was then removed and weighed. Two-dimensional nuclear magnetic resonance technology was used to determine the T1 and T2 spectra after oil saturation without damaging the core.

5. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 4, characterized in that, The determination of T1 and T2 spectra after saturation of oil to establish a two-dimensional NMR spectrum of saturated oil also includes: Software quantifies the amount of oil signal W in different pore sizes of each core under saturated oil conditions. xy，饱和 , where x is the T2 spectral coordinate value and y is the T1 spectral coordinate value.

6. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 3, characterized in that, Determining the T1 and T2 spectra after osmosis to establish a two-dimensional NMR spectrum after osmosis, including: After each core was saturated with oil, it was placed in a piston container, and an adsorption medium was added. The piston container was then placed in a constant temperature chamber set to the formation temperature and pressurized to the reservoir pressure to conduct an adsorption experiment. After day 1 of the experiment, the core was removed and weighed, and the T1 and T2 spectra after adsorption were tested using two-dimensional nuclear magnetic resonance technology.

7. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 6, characterized in that, The determination of post-osmosis T1 and T2 spectra to establish a post-osmosis two-dimensional NMR spectrum also includes: Software quantifies the amount of oil signal W in different pore sizes after permeation of each core sample. xy，渗吸后 , where x is the coordinate value of T2 and y is the coordinate value of T1.

8. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 3, characterized in that, The two-dimensional NMR spectrum oil signal was divided into three regions: bound oil, mesoporous oil, and kerogen. The oil signals in the two-dimensional NMR spectra after saturation and after percolation were divided into bound oil, mesoporous oil, and kerogen regions, respectively.

9. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 8, characterized in that, The oil signals from the two-dimensional NMR spectra after saturation and after percolation were divided into bound oil, mesoporous oil, and kerogen regions, respectively, including: The oil content characteristics of the kerogen region differ from those of bound oil and mesoporous oil. The bound oil, mesoporous oil and kerogen regions can be distinguished by T1, T2 and T1 / T2 ratio.

10. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, The first permeation rate was calculated based on the two-dimensional NMR spectrum, including: Based on the two-dimensional NMR spectrum of the oil signal, the permeation rate I was calculated. xy The formula is as follows: I xy =(W xy，饱和 -W xy，渗吸后 ) / W xy，饱和 Among them, W xy，饱和 W represents the oil signal quantity after saturation. xy，渗吸后 This represents the oil signal quantity after absorption, where x is the T2 coordinate value and y is the T1 coordinate value.

11. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 10, characterized in that, The calculation of the first permeation rate based on the two-dimensional NMR spectrum also includes: For each core sample, the average of the permeability values ​​at the same T1 and T2 coordinates is taken as the first permeability, as shown in the following formula: Among them, I xy,最终 This represents the final absorption rate at points T2 = x and T1 = y; I xy，n Let T2 = x and T1 = y represent the permeation rate at point T2 = x and T1 = y in the nth core sample; n represents the total number of samples.

12. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, The calculation of the second osmotic rate for the kerogen zone includes: Based on two-dimensional NMR spectroscopy of oil signals, software quantization is used to determine the total oil signal quantity W in the kerogen region under saturated oil conditions. 干酪根，饱和 And the total oil signal W in the kerogen region after infiltration 干酪根，渗吸后 The osmotic mobilization efficiency I in the kerogen region is calculated using the following formula. xy，干酪根 :

13. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 12, characterized in that, The calculation of the second osmotic rate for the kerogen zone also includes: The average permeation mobilization efficiency of the kerogen zone in each core is taken as the second permeation rate, as shown in the following formula: Among them, I 干酪根,最终 This indicates the final permeation rate value of the kerogen display zone; I 干酪根,n This represents the percolation rate of the kerogen region in the nth sample; n represents the total number of samples.

14. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Fresh samples of massive shale from the target block were taken, and two-dimensional nuclear magnetic resonance (NMR) spectra were measured to obtain the saturated oil signal quantity of the fresh samples, including: Fresh samples of massive shale from the target block were taken, and T1 and T2 spectra after saturation with oil were measured to establish a two-dimensional NMR spectrum of the fresh samples after saturation with oil. The signal intensity of the saturated oil in the fresh sample was determined based on the two-dimensional NMR spectrum of the saturated oil.

15. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 14, characterized in that, Fresh samples of massive shale from the target block were collected, and T1 and T2 spectra were measured after saturation with oil to establish a two-dimensional NMR spectrum of the fresh samples after saturation with oil, including: Fresh samples of blocky shale from the target area were obtained by drilling using wire cutting technology. The T1 and T2 spectra of the fresh samples were measured using two-dimensional nuclear magnetic resonance technology without damaging the core.

16. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, The total saturated oil signal quantity of the fresh sample is determined based on the saturated oil signal quantity of the fresh sample, including: Software quantifies the amount of oil signal W in different pore sizes of each fresh core sample. xy，新鲜 Then, the oil signal data of all fresh sample oil signal points are summed to obtain the total saturated oil signal W of the fresh sample. 总和，新鲜 , where x is the coordinate value of T2 and y is the coordinate value of T1.

17. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Based on the saturated oil signal intensity and the first permeation rate of the fresh sample, determine the oil signal intensity after permeation into the bound oil and mesopore oil regions of the fresh sample, including: Calculate the oil signal quantity W after seepage at different T1 and T2 coordinates. xy，渗吸后 As shown in the following formula: IN xy，渗吸后 =In xy，新鲜 -IN xy，新鲜 ×I xy,最终 Among them, I xy,最终 W represents the final absorption rate at points T2 = x and T1 = y. xy，新鲜 W represents the oil signal quantity at points T2 = x and T1 = y for fresh samples; xy，渗吸后 This represents the oil signal quantity at points T2=x and T1=y after percolation. W, the bound oil and the mesopore oil zone xy，渗吸后 The oil signal data after all percolation in the data volume are summed to obtain the oil signal W of the bound oil and mesopore oil regions of the fresh sample. 束缚油、中孔油，渗吸后 .

18. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Based on the saturated oil signal intensity and second permeation rate of the fresh sample, the oil signal intensity after permeation in the kerogen region of the fresh sample is determined, including: The oil signal data from the kerogen region of the fresh sample were summed to obtain W. 干酪根，新鲜 The oil signal quantity W after osmosis in the kerogen zone is calculated using the following formula. 干酪根，渗吸后 : IN 干酪根，渗吸后＝ IN 干酪根，新鲜 -IN 干酪根，新鲜 ×I 干酪根,最终 Among them, I 干酪根,最终 W represents the final permeation rate value of the kerogen display zone. 干酪根，新鲜 This represents the total oil signal quantity in the kerogen region of a fresh sample.

19. The method for predicting the adsorption and utilization efficiency of massive shale oil according to claim 1, characterized in that, Based on the total saturated oil signal quantity of fresh samples, the oil signal quantity after permeation in the bound oil and mesoporous oil zones of fresh samples, and the oil signal quantity after permeation in the kerogen zone of fresh samples, the permeation and utilization efficiency of massive shale oil in the target block is calculated, including: The formula for calculating the adsorption mobilization efficiency I of fresh shale oil core samples is as follows: Among them, W 总和，新鲜 W represents the sum of the signal quantities of the fresh sample oil; 干酪根，渗吸后 W represents the total oil signal after permeation in the kerogen oil region; 束缚油、中孔油，渗吸后 This represents the sum of the oil signal quantities after the absorption of bound oil and mesopore oil; I represents the absorption efficiency.

20. A device for predicting the percolation and mobilization efficiency of massive shale oil, characterized in that, include: Unit for measuring the permeation rate of control samples and unit for measuring the permeation mobilization efficiency of fresh samples; The reference sample permeability measurement unit is used to take several core samples of each physical property grade of blocky shale oil. The reference sample permeation rate measurement unit is also used to establish a two-dimensional nuclear magnetic resonance spectrum for each core sample; The reference sample permeation rate measurement unit is also used to divide the two-dimensional NMR spectrum oil signal into bound oil, mesoporous oil region and kerogen region, respectively. The control sample permeation rate measurement unit is also used to calculate the first permeation rate based on the two-dimensional NMR spectrum; and to calculate the second permeation rate for the kerogen region; The fresh sample infiltration and kinetic efficiency measurement unit is used to take fresh samples of blocky shale from the target area, measure the two-dimensional nuclear magnetic resonance spectrum, and obtain the saturated oil signal of the fresh sample; the total saturated oil signal of the fresh sample is determined based on the saturated oil signal of the fresh sample. The fresh sample percolation efficiency measurement unit is also used to determine the oil signal after percolation in the bound oil and mesoporous oil regions of the fresh sample based on the saturated oil signal and the first percolation rate; and to determine the oil signal after percolation in the kerogen region of the fresh sample based on the saturated oil signal and the second percolation rate. The fresh sample percolation efficiency measurement unit is also used to calculate the percolation efficiency of the target block of massive shale oil based on the total saturated oil signal of the fresh sample, the oil signal after percolation of the bound oil and mesopore oil zones of the fresh sample, and the oil signal after percolation of the kerogen zone of the fresh sample.

21. An electronic device, characterized in that, It includes a processor, a communication interface, a memory, and a communication bus, wherein the processor, the communication interface, and the memory communicate with each other through the communication bus; Memory, which stores computer programs; A processor, when executing a computer program stored in memory, implements a method for predicting the percolation and mobilization efficiency of massive shale oil as described in any one of claims 1-19.

22. A computer-readable storage medium storing a computer program, characterized in that, When the computer program is executed by a processor, it implements a method for predicting the adsorption and utilization efficiency of massive shale oil as described in any one of claims 1-19.

Images