3D MODELING METHOD FOR PORE-FILLING HYDRATE SEDIMENT BASED ON CT IMAGE

Abstract

The present invention belongs to the technical field of petroleum exploitation engineering, and discloses a 3D modeling method for pore-filling hydrate sediment based on a CT image. Indoor remolding rock cores or in situ site rock cores without hydrate can be scanned by CT; a sediment matrix image stack and a pore image stack are obtained by gray threshold segmentation; then, a series of pore-filling hydrate image stacks with different saturations are constructed through image morphological processing of the pore image stack such as erosion, dilation and image subtraction operation; and a series of digital rock core image stacks of the pore-filling hydrate sediment with different saturations are formed through image subtraction operation and splicing operation to provide a relatively real 3D model for the numerical simulation work of the basic physical properties of a reservoir of natural gas hydrate.

Claims

1. A 3D modeling method for pore-filling hydrate sediment based on CT image, comprising steps of: step 1, scanning remolding or in situ rock cores without natural gas hydrate by CT to obtain digital rock core image stacks; step 2, adjusting the gray threshold of the digital rock core image stacks, conducting binarization segmentation to obtain a sediment matrix and a pore, and respectively saving as the image stacks; step 3, firstly eroding a pore image stack obtained in step 2 at x pixel and then dilating at x pixel to obtain a pore-filling hydrate image stack; step 4, performing image subtraction; and subtracting the pore-filling hydrate image stack obtained in step 3 from the pore image stack obtained in step 2 to obtain a new pore image stack corresponding to the pore-filling hydrate image stack obtained in step 3; step 5, splicing and combining the sediment matrix image stack obtained in step 2, the pore-filling hydrate image stack obtained in step 3 and the new pore image stack obtained in step 4 to form a digital rock core image stack with the sediment matrix, the pore-filling hydrate and the pore, which is the digital rock core image stack of the pore-filling hydrate sediment; step 6, repeatedly executing step 3 to step 5, and adjusting x value to obtain the digital rock core image stacks of the pore-filling hydrate sediment with different hydrate saturations.

Description

DESCRIPTION OF DRAWINGS

[0014] FIG. 1 shows the binarization results of sediment matrices and pores.

[0015] FIG. 2 is a schematic diagram of image erosion at 4 pixels.

[0016] FIG. 3 is a schematic diagram of image dilation at 4 pixels.

[0017] FIG. 4 shows the result of image subtraction.

[0018] FIG. 5 shows a digital rock core image stack of pore-filling hydrate sediment.

DETAILED DESCRIPTION

[0019] Specific embodiments of the present invention are further described below in combination with accompanying drawings and the technical solution.

Embodiments

[0020] A 3D modeling method for pore-filling hydrate sediment based on a CT image comprises the following steps:

[0021] step 1, scanning remolding rock cores (particle size distribution: 0.01-1 mm; median particle size: 0.15 mm; porosity: 41%) without natural gas hydrate by CT to obtain digital rock core image stacks (resolution: 1024*1024; voxel size: 0.004 mm);

[0022] step 2, adjusting the gray threshold of the digital rock core image stacks, conducting binarization segmentation to obtain a sediment matrix and a pore, and respectively saving as the image stacks, as shown in FIG. 1;

[0023] step 3, firstly eroding a pore image stack obtained in step 2 at 4 pixels, as shown in FIG. 2, and then dilating at 4 pixels, as shown in FIG. 3, to obtain a pore-filling hydrate image stack;

[0024] step 4, performing image subtraction; and subtracting the pore-filling hydrate image stack obtained in step 3 from the pore image stack obtained in step 2 to obtain a new pore image stack corresponding to the pore-filling hydrate image stack obtained in step 3, as shown in FIG. 4;

[0025] step 5, splicing and combining the sediment matrix image stack obtained in step 2, the pore-filling hydrate image stack obtained in step 3 and the new pore image stack obtained in step 4 to form a digital rock core image stack with the sediment matrix, the pore-filling hydrate and the pore, which is the digital rock core image stack of the pore-filling hydrate sediment (saturation: 28.5%), as shown in FIG. 5.