应用伽马测井曲线小波变换计算粒度参数

现代地质 ›› 2012, Vol. 26 ›› Issue (4): 778-783.

应用伽马测井曲线小波变换计算粒度参数

杨宁^1,2, 王贵文^1,2, 赖锦^1,2, 李鉴伦^1,2, 苍丹^1,2, 蒋其君^1,2

1中国石油大学油气资源与探测国家重点实验室，北京102249；
2中国石油大学地球科学学院，北京102249

收稿日期:2012-01-14 修回日期:2012-04-09 出版日期:2012-08-18 发布日期:2012-08-24
作者简介:杨宁，男，博士研究生，1982年出生，层序地层学与测井地质学专业，主要从事层序地层学及测井地质学的研究。Email： B100010011@student.cup.edu.cn。
基金资助:
国家大型油气田及煤层气开发科技重大专项（2011ZX05020-008）。

Application of Gamma Curves Wavelet Transform to Calculate Grain Size Parameters

YANG Ning^-1,2, WANG Gui-Wen^-1,2, LAI Jin^-1,2, LI Jian-Lun^-1,2, CANG Dan^-1,2, JIANG Ji-Jun^-1,2

1State Key Laboratory of Petroleum Resource and Prospecting, China University of Petroleum, Beijing102249, China;
2College of Geosciences, China University of Petroleum, Beijing102249, China

Received:2012-01-14 Revised:2012-04-09 Online:2012-08-18 Published:2012-08-24

摘要/Abstract

摘要：

伽马测井曲线包含了储层的粒度信息，采用二进制小波变换对伽马测井曲线处理可以提取粒度参数。通过分析不同粒度区间单位体积含量(以UddenWentworth方案为标准)与伽马测井曲线变换后信号的幅度、频率之间关系，考虑测井仪器探头测量速率造成的采样间隔的影响，建立求解不同粒度区间体积含量的递归方程组。将获得不同粒度区间的绝对体积转化为相对体积含量，按照Trask或FolkWard公式，计算粒度中值、平均粒度、分选系数、偏度、峰度。选择塔里木盆地大北地区DB101等4口井51个取样点进行方法适用性和可靠性验证，此方法的计算结果与筛选法相比，各项参数误差均在允许接受范围内。与国内外同类方法相比，此方法在信息来源、耗时、成本方面的要求均有所降低。

关键词: 粒度参数, 小波变换, 伽马测井曲线, 岩心, 致密储层

Abstract:

The information on grain size of reservoirs from gamma logging curves can be obtained and the grain size parameter can be extracted based on the dyadic wavelet transform. Considering the influence of sampling intervals caused by the measurement rate of logging tools, an algorithm for solving the system of recurrence equations concerning volume fraction in different kinds of grain size distributions was presented through analyzing the content by volume in different kinds of grain size distribution and the relationship between range and frequency of signal which was extracted after gamma logging curves changed. On the basis of this algorithm, the absolute volume content in different kinds of grain size distribution can be transformed into the relative volume fraction. Then the grain size parameters can be calculated including the median grain diameter, the average particle size, the sorting coefficient, the skewness and the kurtosis according to Trask or FolkWard formulae. Four wells like DB101 in Dabei area of Tarim Basin were selected to verify the applicability and reliability of this method. Moreover, compared with the sieve method, all kinds of parameter errors are within an acceptable range. The advantage of this method is that it can reduce the demand for the source of information, timing consuming and cost in contrast to the similar methods at home and abroad.

Key words: grain size parameter, wavelet transform, gamma logging curve, core, tight reservoir

杨宁, 王贵文, 赖锦, 李鉴伦, 苍丹, 蒋其君. 应用伽马测井曲线小波变换计算粒度参数[J]. 现代地质, 2012, 26(4): 778-783.

YANG Ning, WANG Gui-Wen, LAI Jin, LI Jian-Lun, CANG Dan, JIANG Ji-Jun. Application of Gamma Curves Wavelet Transform to Calculate Grain Size Parameters[J]. Geoscience, 2012, 26(4): 778-783.

参考文献

［1］赵澄林，朱筱敏.沉积岩石学
[M].北京：石油工业出版社，2001：56-60.
［2］Krumbein W C，Pettijohn. F J. Manual of Sediment Petrography
[M]. New York：Appleton Century Crofts,Inc,1938：1-50.
［3］Gibbs R J. The accuracy of particle size analyses utilizing settling tubes
[J].Journal of Sedimentary Petrology，1972，42：41-145.
［4］Agrawal Y C，McCave I N，Riley J B. Laser diffraction size analysis
[M］//Syvitski J M P. Principles，Methods and Application of Particle Size Analysis. Cambridge：Cambridge University Press，1991：119-128.
［5］Verdu J M，Batalla R J，MartinezCasasnovas J A. High resolution grainsize characterization of gravel bars using imagery analysis and geostatistics
[J]. Geomorphology，2005，72(1/4)：73-93.
［6］Patrick L Barnard，David M Rubin. Field test comparison of an autocorrelation technique for determining grain size using a digital beachball camera versus traditional methods
[J]. Sedimentary Geology，2007，201：180-195.
［7］Sarpun I H，Selami K K. Mean grain size determination in marbles by ultrasonic velocity techniques
[J].NDT & E International，2005，38：21-22.
［8］杨姣.延长油田延长组主力油层沉积微相测井识别技术研究
[D].西安：长安大学，2009：11-13.
［9］Daniel Buscombe. Estimation of distributions and associated parameters from digital images of sediment
[J].Sedimentary Geology，2008，210：1-10.
［10］Ruggiero P，Kaminsky G M，Gelfenbaum G，et al. Seasonal to interannual morphodynamics along a highenergy littoral cell
[J]. Journal of Coastal Research，2005，21(3)：553-578.
［11］Rubin D M. A simple autocorrelation algorithm for determining grain size from digital images of sediment
[J].Journal of Sedimentary Research，2004，74(1)：160-165.
［12］Carbonneau P E，Bergeron N，Lane S N. Automated grain size measurements from airborne remote sensing for long profile measurements of fluvial grain sizes
[J]. Water Resources Research，2005，41：W11426.
［13］Verdu J M, Batalla R J, MartinezCasasnovas J A. Estimating grainsize distributions of a gravel riverbed at reach scale from detailed aerial photos, geostatistics and digital image processing(Isabena River, Spain)［R］. Poster: Braided Rivers Conference,2001.
［14］Sheen D R，Johnston L P. Statistical and spatial properties of forest clutter measured with polarimetric Synthetic Aperture Radar (SAR)
[J]. IEEE Transactions on Geoscience and Remote Sensing，1992，30：578-588.
［15］Ulaby F T，Kouyate F，Brisco B. Textural information in SAR images
[J］.IEEE Transaction on Geoscience and Remote Sensing，1986，24：35-245.
［16］Carbonneau P E，Lane S N，Bergeron N E. Catchment scale mapping of surface grain size in gravelbed rivers using airborne digital imagery
[J]. Water Resources Research，2004,40：W07202.
［17］操应长，葸克来，王健，等.砂岩机械压实与物性演化成岩模拟初探［J］.现代地质，2011，25（6）：1152-1158.
［18］Carr J R. Spectral and textural classification of single and multiple band digital images
[J]. Computer and Geosciences，1996，22：849-865.
［19］Haralick R，Shanmugam K，Dinstein I. Textural features for image classification
[J]. IEEE Transactions on Systems, Man and Cybernetics，1973，3：610-621.

[1]	祁诗阳, 沈宸, 刘小标, 张梦娇, 马新明. 基于小波变换的郑州降雨量与太阳黑子活动关系分析[J]. 现代地质, 2023, 37(01): 184-196.
[2]	王志昊, 赵建华, 蒲秀刚, 刘可禹, 李俊乾, 程斌. 页岩油岩心样品洗油实验效率对比分析[J]. 现代地质, 2022, 36(05): 1304-1312.
[3]	赵建鹏, 崔利凯, 陈惠, 李宁, 王自亮, 马瑶, 杜贵超. 基于CT扫描数字岩心的岩石微观结构定量表征方法[J]. 现代地质, 2020, 34(06): 1205-1213.
[4]	陈松, 陈剑, 乔春生, 马俊学, 刘超. 图解法与矩值法在金沙江上游雪隆囊古滑坡堰塞湖溃坝堆积物粒度分析中的应用[J]. 现代地质, 2017, 31(06): 1278-1283.
[5]	马俊学, 陈剑, 崔之久, 刘超, 陈松. 堰塞湖溃坝堆积物的粒度特征分析:以岷江上游叠溪古滑坡堰塞湖为例[J]. 现代地质, 2017, 31(06): 1261-1268.
[6]	陈晨, 卢双舫, 李俊乾, 张鹏飞, 张萍. 不同岩相泥页岩数字岩心构建方法研究——以东营凹陷为例[J]. 现代地质, 2017, 31(05): 1069-1078.
[7]	庄永涛, 刘鹏程, 郝明强. 高温高盐高黏油藏缔合聚合物-表面活性剂二元驱技术[J]. 现代地质, 2017, 31(04): 814-821.
[8]	卢海峰，王金艳，王瑞. 河北省衡水地区第四纪中—晚期标准地层剖面的建立及其地质含义[J]. 现代地质, 2014, 28(5): 962-970.
[9]	屈乐，孙卫，杜环虹，张创，姜黎明，魏虎. 基于CT扫描的三维数字岩心孔隙结构表征方法及应用—以莫北油田116井区三工河组为例[J]. 现代地质, 2014, 28(1): 190-196.
[10]	肖晖1，郭建春1，任山2. 气井水力压裂有效缝长评估及影响因素研究[J]. 现代地质, 2013, 27(6): 1434-1439.
[11]	高小跃，刘洛夫，姜振学，尚晓庆，苏天喜，戴琦雯. 库车坳陷东部YN5井下侏罗统致密储层特征及主控因素[J]. 现代地质, 2013, 27(5): 1099-1019.
[12]	武粤, 孟小红, 刘国峰. 南海北部缘重力异常的多尺度分析及其构造讨论[J]. 现代地质, 2012, 26(6): 1162-1167.
[13]	李灿苹, 刘学伟, 赵桂艳. 利用小波变换分析随机-均匀介质尺度特征[J]. 现代地质, 2010, 24(2): 392-396.
[14]	张凌, 蒋国盛, 宁伏龙, 涂运中, 吴翔, 窦斌. 天然气水合物保真取心装置内部密封技术分析[J]. 现代地质, 2009, 23(6): 1147-1152.
[15]	曹殿华王安建王高尚李瑞萍李以科. 勘查地球化学异常多尺度分析方法——以赣东北德兴矿集区为例[J]. 现代地质, 2008, 22(6): 1028-1033.