CN106501206A - A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology - Google Patents

A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology Download PDF

Info

Publication number
CN106501206A
CN106501206A CN201610040780.1A CN201610040780A CN106501206A CN 106501206 A CN106501206 A CN 106501206A CN 201610040780 A CN201610040780 A CN 201610040780A CN 106501206 A CN106501206 A CN 106501206A
Authority
CN
China
Prior art keywords
rock
anisotropism
different
micro
characteristic
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN201610040780.1A
Other languages
Chinese (zh)
Inventor
周汉国
李静
郭建春
张玉
彭成乐
侯江朋
曲险峰
高长海
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Investment And Development Office Shengli Oilfield Branch Of China Petroleum & Chemical Corp
China University of Petroleum East China
Southwest Petroleum University
Original Assignee
Investment And Development Office Shengli Oilfield Branch Of China Petroleum & Chemical Corp
China University of Petroleum East China
Southwest Petroleum University
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Investment And Development Office Shengli Oilfield Branch Of China Petroleum & Chemical Corp, China University of Petroleum East China, Southwest Petroleum University filed Critical Investment And Development Office Shengli Oilfield Branch Of China Petroleum & Chemical Corp
Priority to CN201610040780.1A priority Critical patent/CN106501206A/en
Publication of CN106501206A publication Critical patent/CN106501206A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N21/00Investigating or analysing materials by the use of optical means, i.e. using sub-millimetre waves, infrared, visible or ultraviolet light
    • G01N21/17Systems in which incident light is modified in accordance with the properties of the material investigated
    • G01N21/25Colour; Spectral properties, i.e. comparison of effect of material on the light at two or more different wavelengths or wavelength bands
    • G01N21/31Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry
    • G01N21/35Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light
    • G01N21/3563Investigating relative effect of material at wavelengths characteristic of specific elements or molecules, e.g. atomic absorption spectrometry using infrared light for analysing solids; Preparation of samples therefor

Landscapes

  • Physics & Mathematics (AREA)
  • Spectroscopy & Molecular Physics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigating Or Analysing Materials By Optical Means (AREA)

Abstract

The invention discloses a kind of method of utilization micro ft-ir spectroscopy imaging technique accurate Characterization oil and gas reservoir rock anisotropism, it is that Comprehensive Evaluation of Reservoirs and structure stress scene simulation method modeling provide more accurate effectively information, preferably instructs oilfield prospecting developing.The anisotropism that reservoir rock is characterized at spectrum and two, space aspect using micro ft-ir spectroscopy Image analysis system, different rock mediums are differentiated to the response characteristic of infrared light by rock, using different medium to the absorbance A value of characteristic wavelength as characteristic quantity, by setting its threshold value, different rock mediums are separated, the information such as species, hole, the size in crack and its spatial distribution of rock medium are obtained, the accurate Characterization to oil and gas reservoir rock anisotropism is realized.It is demonstrated experimentally that rock anisotropism characterizing method of the rock anisotropism characterizing method ratio based on micro ft-ir spectroscopy imaging technique based on visible range color digital image treatment technology has higher practical value.

Description

A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology
Technical field:
The invention belongs to petroleum geology exploration technical field;Specifically related to a kind of anisotropism characterizing method of oil and gas reservoir rock, particularly a kind of utilization micro ft-ir spectroscopy imaging technique is come the Comprehensive Analysis Technique method of accurate Characterization oil and gas reservoir rock anisotropism.
Background technology:
Multi- scenarios method is acted on while causing Large-scale Stratum structure change, on micro-scale, oil and gas reservoir rock is able to be formed or is changed so as to present strong anisotropism.The anisotropism of research oil and gas reservoir rock, significant to Comprehensive Evaluation of Reservoirs.Nineteen nineties, forefathers are randomly generated the different constituents of rock according to Principle of Statistics and probability distribution function, form rock two-dimensional section or said three-dimensional body, carry out the anisotropism of study of rocks.Although the method can characterize the anisotropism of rock to a certain extent, differ larger with rock truth.In recent years, with the development of digital image processing techniques, numerous scholars start with digital picture to characterize the anisotropism of rock.Image is made up of pixel, and wherein each pixel has a unique numerical value to represent the color of the pixel, i.e., it is single uniform that each pixel can only represent a kind of material at color, namely each pixel or chemical composition.Due to containing different mineral grains in rock, these mineral grains typically have different colors.Therefore, it can, during Digital Image Processing, the pixel with same pixel value or pixel value within a certain range is divided into a class, represent a kind of material or chemical composition, it is also possible to for representing the structure and defect of material internal.But, anisotropism characterizing method based on the rock of Digital Image Processing, simply presenting different colors in the picture according to rock medium will be separated for the constituent of rock, it is impossible to avoid due to the drawbacks of different constituent performances causes error for same color.
At present, micro ft-ir spectroscopy has become most widely used spectral analysis technique, is referred to as the giant for analyzing.Micro ft-ir spectroscopy has the features such as fast analyze speed, efficiency high, low cost, nondestructive analysis.First, the spatial resolution (single pixel dimension) of micro ft-ir spectroscopy imaging is 25 microns or 50 microns, this concrete dimensional properties heterogeneous with oil and gas reservoir rock matches, you can to think that material at each pixel or chemical composition are single uniform.Secondly, spectrum analysis is that there are different microscopic characteristics, microscopic characteristics to determine difference response of the material to infrared light again, be embodied according to different material:There is different material different spectrum picture and characteristic absorption wavelength, micro ft-ir spectroscopy can characterize material property on a microscopic scale, carry out accurately Identification to material.In sum, although the spatial resolution of micro ft-ir spectroscopy imaging is slightly poorer than in visible spectrum wave band area, but the characteristic spectrum resolution ratio in rock and sensitivity and is substantially better than visible spectrum.Therefore, characterized based on the rock anisotropism that micro ft-ir spectroscopy is imaged and characterized more accurately, more effectively, with more use value than the rock anisotropism based on digital image processing techniques.
The present invention is based on micro ft-ir spectroscopy imaging technique, obtains the accurate physical property characteristic of unhomogeneous rock, is that anisotropism, seepage characteristic, stress characteristics, destruction property and the analysis Reservoir Fracture for studying oil and gas reservoir rock forms the accurate foundation of the offer such as rule.
Content of the invention:
It is contemplated that carrying out accurate Characterization using micro ft-ir spectroscopy imaging technique to the anisotropism of oil and gas reservoir rock, from rock core micro ft-ir spectroscopy image, the information such as species and spatial distribution and hole, the size in crack and the distribution of the rock medium are particularly extracted, more accurate, effective information is provided for the modeling of the numerical simulations such as the anisotropism sign of oil and gas reservoir rock, material properties of rock, microscopic seepage characteristic.
The present invention solves the scheme taken of its technical problem:Carry out white light transmittant first under the microscope, obtain the panchromatic smooth black white image of rock core mating plate, and measured zone is selected with reference to the image;Then the measured zone that selectes is measured using micro ft-ir spectroscopy Image analysis system, obtain the size in crack and distribution on the infrared spectroscopic imaging and the rock core mating plate in the region, and spectral signature point analysis is carried out to the region, the species of rock medium is recognized with the characteristic wavelength of different rock mediums according to the sharp peaks characteristic of characteristic point infrared spectrum curve;Degree of absorption also dependent on rock core mating plate to infrared light, i.e. absorbance A recognize the information of relevant information and different rock medium contents in hole present in rock core mating plate and crack etc. simultaneously;Finally using rock core mating plate to the algebraic expression of the degree of absorption (absorbance A value) of one or more characteristic wavelength infrared lights as characteristic quantity, by setting its threshold value, will be separated for different rock mediums, realize the accurate Characterization to oil and gas reservoir rock anisotropism.
In the present invention, the anisotropism of heterogeneous reservoir rock sample is analyzed in spectrum and two, space aspect using micro ft-ir spectroscopy Image analysis system and is characterized, in particular with infrared microscopy technology, different rock mediums are differentiated to the response characteristic of infrared light not only by rock sample, can also using different medium to the degree of absorption (absorbance A) of characteristic wavelength as partition threshold, different rock mediums are separated, the information such as the spatial distribution of different medium are obtained, so as to more accurately characterize the anisotropism of oil and gas reservoir rock.It is demonstrated experimentally that rock anisotropism characterizing method of the rock anisotropism characterizing method ratio based on micro ft-ir spectroscopy imaging technique based on visible range color digital image treatment technology has more use value.
Description of the drawings:
Fig. 1 is the implementing procedure block diagram of the reservoir rock anisotropism characterizing method based on micro ft-ir spectroscopy imaging technique proposed by the invention.
Fig. 2 is the panchromatic smooth black-and-white photograph of rock core mating plate and measured zone schematic diagram.
It is 4000cm to wave-number range that Fig. 3 is measured zone-1-1400cm-1Infrared light mean light absorbency figure.
Fig. 4 is the infrared spectrum curve figure of 6 characteristic points.
It is 2530cm to wave number that Fig. 5 is measured zone-1Infrared light unicast number absorbance figure.
Fig. 6 is the change curve of the absorbance A value of each point on y=-710 straight lines in measured zone.
Fig. 7 is the 4 kinds of different medium area schematics divided according to characteristic quantity A values.
Fig. 8 is the vector expression figure of 4 kinds of different medium area schematics, provides reliable data foundation for numerical simulation modeling.
Specific embodiment:
The present invention is described in detail with reference to embodiment.
1st, prepared by rock core mating plate
Experiment is made mating plate (thick 30 microns) with core sample by national oil Gas Industry standard SY/T5913-2004 " rock flaking method ", is bonded in glass slide.
2nd, the panchromatic smooth black white image of rock core mating plate is obtained
The rock core mating plate for having prepared is carried out white light transmittant under fluorescence microscope and retains digital image recording, obtain the panchromatic smooth black white image of rock core mating plate, and measured zone is selected with reference to the image, as shown in Figure 2.The region be can be seen that from selected region obvious black region and white portion, illustrates in the region that black region is main absorbance region, and white portion is larger due to white light transmittant rate, may be to be not filled by crack, pore region.
3rd, survey region micro ft-ir spectroscopy measurement
Micro ft-ir spectroscopy experiment instrument is Fourier transform infrared/near infrared imaging system, using 25 μm of resolution ratio, infrared transmission is carried out to the measured zone that chooses on visible ray black white image, the infrared spectroscopic imaging with different absorbances is obtained, as shown in Figure 3.There are Fig. 3 different absorbances (to be represented with different colors, absorbance intensity is as shown in color-bar on the right of Fig. 3), different selected characteristic points 1,2,3,4,5,6 according to absorbance obtain the infrared spectrum curve of 6 characteristic points, as shown in Figure 4 (figure is transmittance graph figure, and its size of transmissivity T and absorbance A is contrary).
4th, characteristic point curve of spectrum analysis
As shown in Figure 4, the trend of 6 curves of spectrum is unanimous on the whole, in 3200cm-1-3650cm-1There is obvious absworption peak (characteristic absorption peak of the wave band for hydroxyl) in wave band, but zones of different is different to the infrared Absorption intensity of the wave band, and the hydroxyl containing obvious uneven distribution in the rock core mating plate is described;As Spectrum4 with Spectrum5 variation tendencies, transmissivity T is very big (absorbance A is little), only one of which more significantly hydroxyl group absorption wave band 3200cm-1-3650cm-1, illustrate that region rock medium representated by characteristic point 4 and 5 is few, be hole and crack region (blueness, purple, black region in Fig. 3) on the rock mating plate;, as Spectrum3 variation tendencies, simply transmissivity T is not of uniform size for Spectrum1, Spectrum2, except 3200cm-1-3650cm-1This obvious absorption bands, in 2899cm-1Wave number and 2530cm-1There are stronger absworption peak, wherein 2530cm at wave number-1Wave number is the characteristic absorption peak of limestone, illustrate that the Region Medium representated by characteristic point 1,2,3,6 is mainly dolomite, but representated by 1,2,3,6 characteristic points, the absorbance A in region is not of uniform size in figure 3, illustrate that the region dolomite content representated by characteristic point 1,2,3,6 is different, therefore, region representated by characteristic point 1,2,3,6 can be divided into by different areas of dielectric according to the size of absorbance A.
5th, rock medium is extracted
According to characteristic point curve of spectrum analytical conclusions, the rock core mating plate is extracted in " spectral cube " to 2530cm-1The absorbance image of wave number, as shown in Figure 5.In Figure 5, along y=-710 draw a horizontal line, with the straight line each point X-coordinate as abscissa, with each point to 2530cm-1The absorbance A of wave number is ordinate, draws out the change curve of each materials absorbed light degree A on the straight line, as shown in Figure 6.It will be appreciated from fig. 6 that the change of absorbance A is consistent with the change of rock medium, along the line, different dolomite concentration areas A values are different.According to the difference of A values, the rock core mating plate is divided into 4 regions:1. hole, 3. crack area (A≤0.5800), 2. hole and rock medium transitional region (0.5800 < A≤1.1547), dolomite areas of dielectric (1.1547 < A≤1.9162), 4. more complete dolomite areas of dielectric (A > 1.9162), 4 kind regions are individually represented with different colors, as shown in Figure 7.
6th, image vector
4 kinds of regions that Fig. 7 is represented are present in the form of a bitmap, and this form is unfavorable for that Deterministic Finite meta-model is set up in the work such as follow-up numerical simulation.Therefore need bitmap to be converted into polar plot, the data of FEM model are set up in convenient acquisition.Polar plot is to describe figure using straight line and curve, and the element of these figures is some point, line, rectangle, polygon, circle and camber lines etc., and they are to calculate to obtain by mathematical formulae.Each object in polar plot is an entity having a style of one's own, with attributes such as color, shape, profile, size and locations.Image vector program is worked out using software Matlab, vector quantization is carried out to Fig. 7, as a result as shown in Figure 8.
Fig. 8 contains the shape in each piece of region in composition rock section, size and location information data, the data in any one piece of region can be extracted from Fig. 8, information (position) that puts in the information (size, position) in face, the information (length, position) of boundary line and boundary line including the composition region etc..According to the actual requirements, corresponding data are extracted, it is possible to set up the model of accurate Characterization rock anisotropism.For example, the stress characteristics of unhomogeneous rock to be studied, when FEM model is set up, the data boundary (position, length) in the every kind of region generated by Fig. 8 can be just extracted, according to data boundary, every kind of region is generated in corresponding finite element software, the rock section is constituted, then grid division, being further applied load carries out stress analysis.

Claims (3)

1. a kind of based on micro ft-ir spectroscopy be imaged reservoir rock anisotropism characterizing method, first under the microscope White light transmittant is carried out, and is obtained the panchromatic smooth black white image of rock core mating plate, and measured zone is selected with reference to the image;Then The measured zone that selectes is measured using micro ft-ir spectroscopy Image analysis system, obtain the infrared spectrum in the region The size in crack and distributed intelligence on image and the rock core mating plate;And spectral signature point analysis is carried out to the region, according to The characteristic wavelength of the sharp peaks characteristic of characteristic point infrared spectrum curve and different rock mediums is recognizing the species of rock medium; Hole present in rock core mating plate, crack and different rocks are recognized according to rock core mating plate to the degree of absorption of infrared light simultaneously The information such as medium content;Last degree of absorption (the absorbance by rock core mating plate to one or more characteristic wavelength infrared lights A values) algebraic expression as characteristic quantity, by setting its threshold value, will be separated for different rock mediums, realize to oil The accurate Characterization of gas reservoir rock anisotropism.
2. according to claim 1 based on micro ft-ir spectroscopy be imaged reservoir rock anisotropism characterizing method, It is characterized in that:Infrared spectrum imaging is carried out to rock mating plate using micro ft-ir spectroscopy imaging technique, in infrared spectrum In image, selected characteristic point carries out characteristic point infrared spectrum curve analysis, special according to the peak value of characteristic point infrared spectrum curve Seek peace different rock mediums characteristic wavelength recognizing the species of rock medium.
3. according to claim 1 based on micro ft-ir spectroscopy be imaged reservoir rock anisotropism characterizing method, It is characterized in that:Different to the degree of absorption (absorbance A value) of infrared light according to different rock mediums, by rock core light Piece, leads to as characteristic quantity to the algebraic expression of the degree of absorption (absorbance A value) of one or more characteristic wavelength infrared lights Cross and set its threshold value, will be separated for different rock mediums, realize the accurate Characterization to oil and gas reservoir rock anisotropism.
CN201610040780.1A 2016-01-21 2016-01-21 A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology Pending CN106501206A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610040780.1A CN106501206A (en) 2016-01-21 2016-01-21 A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610040780.1A CN106501206A (en) 2016-01-21 2016-01-21 A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology

Publications (1)

Publication Number Publication Date
CN106501206A true CN106501206A (en) 2017-03-15

Family

ID=58287225

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610040780.1A Pending CN106501206A (en) 2016-01-21 2016-01-21 A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology

Country Status (1)

Country Link
CN (1) CN106501206A (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110487794A (en) * 2019-09-20 2019-11-22 鞍钢集团矿业有限公司 A kind of homogeneity recognition methods of the rock strength based on Digital image technology
CN112881328A (en) * 2021-02-26 2021-06-01 华北水利水电大学 Altered rock strength prediction method based on short-wave infrared spectroscopy technology

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101339131A (en) * 2008-08-13 2009-01-07 中国石油天然气集团公司 Rock core microscopic various light spectrum image-forming information comprehensive processing method
CN101573606A (en) * 2006-08-09 2009-11-04 昆士兰大学 Characterisation of geological materials by thermally induced response

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101573606A (en) * 2006-08-09 2009-11-04 昆士兰大学 Characterisation of geological materials by thermally induced response
CN101339131A (en) * 2008-08-13 2009-01-07 中国石油天然气集团公司 Rock core microscopic various light spectrum image-forming information comprehensive processing method

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
宋宁: "非均匀介质的光学和红外光谱性质——对含油岩屑、岩心研究的应用", 《中国博士学位论文全文数据库 基础科学辑》 *
李静等: "光谱成像信息的数据融合技术在储层表征中的应用", 《光谱学与光谱分析》 *
黄乔松等: "含油岩心显微红外光谱成像方法的研究", 《光谱学与光谱分析》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110487794A (en) * 2019-09-20 2019-11-22 鞍钢集团矿业有限公司 A kind of homogeneity recognition methods of the rock strength based on Digital image technology
CN110487794B (en) * 2019-09-20 2021-10-22 鞍钢集团矿业有限公司 Rock strength homogeneity degree identification method based on digital image technology
CN112881328A (en) * 2021-02-26 2021-06-01 华北水利水电大学 Altered rock strength prediction method based on short-wave infrared spectroscopy technology
CN112881328B (en) * 2021-02-26 2024-03-22 华北水利水电大学 Method for predicting strength of altered rock based on short-wave infrared spectrum technology

Similar Documents

Publication Publication Date Title
CN102854531B (en) Multi-scale rock physical plate method and device for reservoir hydrocarbon detection
Yuan et al. Soil moisture retrieval model for remote sensing using reflected hyperspectral information
Tusa et al. Mineral mapping and vein detection in hyperspectral drill-core scans: Application to porphyry-type mineralization
CN105308588B (en) Effective contour and Jian Men in flow cytometer
Reinhardt et al. Benchmarking conventional and machine learning segmentation techniques for digital rock physics analysis of fractured rocks
CN104885124A (en) Method for producing a three-dimensional characteristic model of a porous material sample for analysis of permeability characteristics
Wang et al. Study on the pore structure and fractal dimension of tight sandstone in coal measures
Wang et al. Determination of NMR T2 cutoff and CT scanning for pore structure evaluation in mixed siliciclastic–carbonate rocks before and after acidification
Visalli et al. Micro-Fabric Analyzer (MFA): A new semiautomated ArcGIS-based edge detector for quantitative microstructural analysis of rock thin-sections
Zhang et al. Permeability prediction of carbonate rocks based on digital image analysis and rock typing using random forest algorithm
Elsayad Spectral phasor analysis for Brillouin microspectroscopy
Pavičić et al. Geometric and fractal characterization of pore systems in the Upper Triassic dolomites based on image processing techniques (example from Žumberak Mts, NW Croatia)
Liu et al. X-ray micro CT based characterization of pore-throat network for marine carbonates from South China Sea
Zhao et al. Research on micro-pore structure and 3D visual characterization of inter-salt shale based on X-CT imaging digital core technology
Frank et al. Experimental reproducibility and natural variability of hydraulic transport properties of fractured sandstone samples
CN106501206A (en) A kind of reservoir rock anisotropism characterizing method based on infrared microscopy technology
Li et al. MV-GPRNet: multi-view subsurface defect detection network for airport runway inspection based on GPR
Mayerich et al. Real-time interactive data mining for chemical imaging information: application to automated histopathology
Kang et al. Construction of complex digital rock physics based on full convolution network
Xue et al. Insights into heterogeneity and representative elementary volume of vuggy dolostones
Xie et al. Investigation of the Oriented Structure Characteristics of Shale Using Fractal and Structural Entropy Theory
Zhang et al. Application of multifractal theory for determination of fluid movability of coal-measure sedimentary rocks using nuclear magnetic resonance (NMR)
Lopes et al. Interactive three-dimensional chemical element maps with laser-induced breakdown spectroscopy and photogrammetry
Lang et al. Three-dimensional finite element simulation and reconstruction of jointed rock models using CT scanning and photogrammetry
Wang et al. Study on multi-scale coral reef pore structure identification and characterization

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20170315

WD01 Invention patent application deemed withdrawn after publication