A new method for characterizing heterogeneities from a core image using local Hölder exponents

The two-dimensional multifractional Brownian motion (2D-mBm) provides a convenient model for image processing. However, the accurate estimation of its associated local regularity function presents a serious difficulty. Here, we suggest a new estimator of the local Hölder exponent (H) of 2D-mBm paths. The idea consists of extending the algorithm, as suggested by Peltier and Lévy-Véhél (1995) for quantifying regularity of mono-dimensional signals, to the 2D case. The resulting algorithm is first tested on artificial 2D-mBm paths shaped by the kriging method. The good accuracy of the obtained H estimations motivates us to implement it on digitalized image data of a core extracted from an Algerian borehole. The derived regularity map shows a significant correlation with the geological features observed on the investigated core. These lithological discontinuities are marked by local variations of the Hölder exponent value. Besides, it is shown that no clear relationship can be drawn between regularity and digitalized data. To conclude, the proposed method may be of high importance in characterizing heterogeneities from core images using the regularity exponents.     

月表虹湾Hurst指数分析

月球遥感是在登月技术还不完善的情况下,利用遥感数据,通过各种分析来研究月球演化、内部结果和表面特征的方法。粗糙度是地表隆升、沉降、侵蚀和火山喷射物覆盖在行星表面留下的记录,因此对行星表面粗糙度的定量分析可以精确反映其地质演化过程。而Hurst指数则是粗糙度的精确反映。文中利用LOLA制作的DEM,以月球正面虹湾为例,选取水平方向的12条剖面和垂直方向的15条剖面对其各粗糙度参数进行计算,并结合月表成熟度对Hurst指数与底层年龄的相关性进行了分析,得出如下结论:(1)在1km的剖面长度上,虹湾地区均方根高程的平均值在3m左右;在0.2～3km研究尺度范围内,虹湾地区的坡度不超过2°,说明该地表起伏度小。(2)虹湾地区的Hurst指数跨度比较大,尤其在垂直方向(南北)上,在0.4到0.9之间,表明其构造样式丰富。(3)在0.2～3km研究尺度范围内,Hurst指数自北往南增大、自西往东增大,表示在虹湾地区粗糙度自西北向东南方向增大。(4)虹湾地区地质年龄为自西北向东南方向变小,与Hurst指数增大方向一致。(5)在月表年轻撞击坑地区月表成熟度与Hurst指数大致呈负相关关系。     

Three-Dimensional Modeling of Shallow Shear-Wave Velocities for Las Vegas, Nevada, Using Sediment Type

A three-dimensional model of near-surface shear-wave velocity in the deep alluvial basin underlying the metropolitan area of Las Vegas, Nevada (USA), is being developed for earthquake site response projections. The velocity dataset, which includes 230 measurements, is interpolated across the model using depth-dependent correlations of velocity with sediment type. The sediment-type database contains more than 1 400 well and borehole logs. Sediment sequences reported in logs are assigned to one of four units. A characteristic shear-wave velocity profile b developed for each unit by analyzing closely spaced pairs of velocity profiles and well or borehole logs. The resulting velocity model exhibits reasonable values and patterns, although it does not explicitly honor the measured shear-wave velocity profiles. Site response investigations that applied a preliminary version of the velocity model support a two-zone ground-shaking hazard model for the valley. Areas in which clay predominates in the upper 30 m are predicted to have stronger ground motions than the rest of the basin.     

克拉美丽气田石炭系火山岩复杂岩性岩电特征

准噶尔盆地东部克拉美丽气田石炭系火山岩岩石类型复杂多样,且同一岩性由于结构、构造和成分的差异,电性特征差异亦较大,岩性识别困难。本文通过对该区14口井取心段岩电关系研究,认为自然伽马、电阻率、密度三种曲线岩性特征响应明显; 电阻率、中子、密度、声波四条曲线对火山岩岩石构造特征响应明显。并编制了岩性和岩石构造测井识别交会图版11张。进而利用电测资料识别出该区11种火山岩岩石类型: 正长斑岩、二长斑岩、玄武岩、粗面岩、英安岩、流纹岩、霏细岩、沉凝灰岩、熔结凝灰岩、火山角砾岩和熔结火山角砾岩; 识别出5种岩石构造类型: 正长斑岩中气孔及块状构造和玄武岩中杏仁、碎裂及块状构造。通过本区12口钻井取心后验,测井识别结果与钻井岩心分析结果吻合良好,可作为地区性火山岩测井岩性、岩石构造识别模式。     

大型坳陷湖盆层序地层格架的研究方法体系——以鄂尔多斯盆地中生界延长组为例

摘 要 大型坳陷湖盆层序地层学研究有很大进展,但鄂尔多斯大型湖盆层序地层研究进展缓慢,主要原因是黄土塬地震资料品质差,钻井资料多解性强,横向等时对比困难,尚未找到合适的研究方法体系。本文通过三年攻关研究,总结出适于这种湖盆层序地层格架建立的方法体系,笔者称之为“六定”方法体系:“露头剖面定标准、小波变换定级次、湖平面变化定旋回、测井曲线定样式、地震解释定格架、三者结合定方案”。通过“六定”方法体系的系统应用,建立了可以在全盆地范围内进行地层划分与对比的三级层序地层格架。     

Fracture Length Estimation from Borehole Image Logs

We describe a method to estimate fracture length for circular fractures from borehole image logs. The relative frequency of fractures, which have complete circumference trace on image logs is related to fracture length. A simple functional relationship can be derived for the relative frequency of complete fracture traces in terms of average fracture inclination to borehole, borehole diameter and fracture length. This formulation however, tends to underestimate fracture length because a constant length is assumed. A more accurate length estimate can be obtained by assuming that fracture length is linearly correlated to fracture aperture or spacing. Cumulative frequency distribution of fracture aperture and spacing can be obtained from borehole image logs. The problem then transforms itself to finding the scaling factor between fracture length and aperture (or spacing) from the relative frequency of fractures with complete traces. The product of the scaling factor and average fracture spacing (or aperture) gives the average fracture length.     

用测井曲线识别地质标志层位

利用动态规划算法进行有序序列匹配来识别地质标志层位,在选取特征模式参数时,不仅考虑了曲线形态的特征,而且考虑了地层岩性的特征。对特征模式向量进行了量纲一致性处理。识别时既可以对全孔所有层识别,也可以只做部分层位识别,还可以分段识别。对山东唐口地区的几口井做了识别处理,均取得较满意的结果。      

Détection des transitions lithologiques par l'analyse de la composante fractale des diagraphies par transformée continue en ondelettes

The frequency analysis of many log data permits to verify that their stochastic component show ‘power-law-type’ spectral densities, characteristic of $1/f$ noise. They can be modelled by fractional Brownian motions. Continuous Wavelet Transformation (CWT) provides us with very efficient methods to determine the local spectral exponents of these scaling laws. These new attributes are related to the local fractality of these signals. We first present some theoretical results and an application to a fractional Brownian motion. The second application concerns a dataset recorded in the MAR203 borehole. We show that clustering of these new pseudo-logs leads to a good resolution between different lithofacies. To cite this article: N. Zaourar et al., C. R. Geoscience 338 (2006).     

松辽盆地科学钻探工程松科二井东孔上白垩统地球物理测井科学成果

松辽盆地科学钻探工程是围绕地球深部资源、古环境与古气候等一系列重要地球科学问题而实施的重大科学工程。作为该工程的主体钻孔, 松科二井东孔设计为松辽盆地最深钻孔, 设计深度6 400.0 m, 计划穿透白垩系地层、钻至盆地基底。由于上白垩统井段取心量非常少, 因此地球物理测井的作用尤为重要。该孔正式开钻后, 采用先进的成像测井设备对一开和二开钻井井段(上白垩统井段)先后实施4次裸眼井综合测井和套管井测井作业, 并系统地开展测井综合评价, 取得以下科学成果: (1)原位获取了常规测井、特殊测井和固井质量检查测井资料共20种, 资料丰富、质量优良, 建立了上白垩统连续的岩石物理参数标准剖面; (2)完成了上白垩统地层划分、岩性识别、沉积和构造精细解释; (3)识别出泉头组油气储层、嫩江组和青山口组烃源岩层、嫩江组高放射性异常层, 以及上白垩统地热开发潜力层; (4)探索了青山口组古气候分析的有利测井指标; (5)及时提供准确的井径、钻孔轨迹和固井质量等关键钻探工程信息。研究成果对于松辽盆地地球物理勘探科学“标尺”建立、古环境与古气候研究、钻探工程施工起到至关重要的作用, 为整个钻探工程科学目标的实现奠定良好基础。     

High-resolution reflection seismic imaging of the upper crust at Laxemar, southeastern Sweden

A major cost in exploring the upper 1–2 km of crystalline crust with reflection seismics is the drilling required for explosive sources. By reducing the charge size to a minimum, shallow inexpensive shotholes can be drilled with handheld equipment. Here, we present results from a full-scale test using small charges for high-resolution seismic surveying over a nuclear waste disposal study site (not an actual site). Two 2–2.5-km-long crossing profiles were acquired in December 1999 with 10-m shot and geophone spacing in the Laxemar area, near Oskarshamn in southeastern Sweden. After standard processing, including dip moveout (DMO), several subhorizontal to moderately dipping reflections are imaged. Many of the dipping ones can be correlated to fracture zones observed in a ca. 1700-m-deep borehole where the profiles cross and/or to fracture zones mapped on the surface. The imaged fracture zones form a complex 3D pattern illustrating the necessity of having 3D control before interpreting seismic reflection data. Analyses of sonic and density logs from the borehole show that greenstones have significantly higher impedances than the more dominant granite found in the borehole (granite/greenstone reflection coefficient is +0.065). These greenstones may contribute to the reflectivity when associated with fracture zones. In some cases, where they are present as larger subhorizontal lenses, they may be the dominant source of reflectivity. A set of north-dipping (10°) reflectors at 3–3.5-km depth can be correlated to a similar set observed below the island of Ävrö about 3 km to the east.     

Fluid Property Discrimination in Dolostone Reservoirs Using Well Logs

The Ordovician Majiagou Formation is one of the main gas-producing strata in the Ordos Basin,China.The identification of hydrocarbon-bearing intervals via conventional well logs is a challenging task.This study describes the litholog of Ma 5(Member 5 of Majiagou Formation)dolostones,and then analyzes the responses of various conventional well logs to the presences of natural gas.The lithology of the gas bearing layers is dominantly of the dolomicrite to fine to medium crystalline dolomite.Natural gas can be produced from the low resistivity layers,and the dry layers are characterized by high resistivities.Neutron-density crossovers are not sensitive to the presences of natural gas.In addition,there are no significant increases in sonic transit times in natural gas bearing layers.NMR(nuclear magnetic resonance)logs,DSI(Dipole Sonic Imager)logs and borehole image logs(XRMI)are introduced to discriminate the fluid property in Majiagou dolostone reservoirs.The gas bearing intervals have broad NMR T_2(transverse relaxation time)spectrum with tail distributions as well as large T_(2gm)(T_2 logarithmic mean values)values,and the T_2 spectrum commonly display polymodal behaviors.In contrast,the dry layers and water layers have low T_(2gm) values and very narrow T_2 spectrum without tails.The gas bearing layers are characterized by low V_p/V_s ratios,low Poisson’s ratios and low P-wave impedances,therefore the fluid property can be discriminated using DSI logs,and the interpretation results show good matches with the gas test data.The apparent formation water resistivity(AFWR)spectrum can be derived from XRMI image logs by using the Archie’s formula in the flushed zone.The gas bearing layers have broad apparent formation water resistivity spectrum and tail distributions compared with the dry and water layers,and also the interpretation results from the image logs exhibit good agreement with the gas test data.The fluid property in Majiagou dolostone reservoirs can be discriminated through NMR logs,DSI logs and borehole image logs.This study helps establish a predictable model for fluid property in dolostones,and have implications in dolostone reservoirs with similar geological backgrounds worldwide.     

The Little Ice Age signature and subsequent warming seen in borehole temperature logs versus solar forcing model

The ‘low’ in the transient temperature versus depth borehole profiles around 120 m seen from deep temperature logs in the Canadian Prairies (southern Alberta–southern Saskatchewan), as well as in some of the European data, has been interpreted to be related to the Little Ice Age (LIA). Data point to the lowest ground surface and subsurface temperatures occurring in the very late eighteenth to nineteenth centuries. Inversion of these logs shows that surface temperature lows were followed by a recent warming period. Further, the synthetic profiles built on the basis of solar forcing history, stretching as far back as the beginning of the seventeenth century, suggest that the LIA signatures interpreted from the inversion of the borehole temperature logs would be difficult to be explained by known published models of past solar irradiation despite large range of assumed sensitivities for the couplings assumed, and that further forcing needs to be considered.     

裂缝性碳酸盐岩储层声波时差曲线的波动和增幅分析

研究了裂缝性碳酸盐岩储层声波时差曲线的波动和增幅现象,认为是裂缝对声波的衰减导致到达远接收器的初至波幅度降低,因而形成3～4μs/ft的时差增幅,从而有效解释了声波时差曲线的波动和增幅现象。利用成像资料和常规测井资料相结合,研究增幅发生条件和规律,发现低角度、较大开度的裂缝导致时差增幅。     

Bias in estimating fractal dimension with the rescaled-range (R/S) technique

Fractal geostatistics are being applied to subsurface geological data as a way of predicting the spatial distribution of hydrocarbon reservoir properties. The fractal dimension is the controlling parameter in stochastic methods to produce random fields of porosity and permeability. Rescaled range (R/S)analysis has become a popular way of estimating the fractal dimension, via determination of the Hurst exponent (H). A systematic investigation has been undertaken of the bias to be expected due to a range of factors commonly inherent in borehole data, particularly downhole wireline logs. The results are integrated with a review of previous work in this area. Small datasets. overlapping samples, drift and nonstationariry of means can produce a very large bias, and convergence of estimates of H around 0.85–0.90 regardless of original fractal dimension. Nonstationarity can also account for H>1, which has been reported in the literature but which is theoretically impossible for fractal time series. These results call into question the validity of fractal stochastic models built using fractal dimensions estimated with the R/Smethod.     

Bias in estimating fractal dimension with the rescaled-range (R/S) technique

Fractal geostatistics are being applied to subsurface geological data as a way of predicting the spatial distribution of hydrocarbon reservoir properties. The fractal dimension is the controlling parameter in stochastic methods to produce random fields of porosity and permeability. Rescaled range (R/S)analysis has become a popular way of estimating the fractal dimension, via determination of the Hurst exponent (H). A systematic investigation has been undertaken of the bias to be expected due to a range of factors commonly inherent in borehole data, particularly downhole wireline logs. The results are integrated with a review of previous work in this area. Small datasets. overlapping samples, drift and nonstationariry of means can produce a very large bias, and convergence of estimates of H around 0.85–0.90 regardless of original fractal dimension. Nonstationarity can also account for H>1, which has been reported in the literature but which is theoretically impossible for fractal time series. These results call into question the validity of fractal stochastic models built using fractal dimensions estimated with the R/Smethod.     

Plate motion predictability using Hurst exponent applied to the Maitri-antarctica GPS network

The Global Positioning System (GPS) data collected during years 1997–2007 at Maitri along with the International GNSS Service (IGS) stations viz. Casey, Davis, Yaragadee, Tidbinbilla, Seychelles, Coco, Diego Garcia, Kerguelen and Hartbeesthoek in and around Indian Ocean were processed using Bernese. 4.2 software. The baseline length and their changes were estimated and spectral analysis was carried out subsequently. The value of the Hurst exponent, as estimated from the slope of spectral fall-off of each time series of baseline length, was used to interpret the characteristics of the plate motions in the region. The Hurst exponent is found to be 0.65 (>0.5) for the baseline length between stations Kerguelen and Coco clearly showing a constant increasing trend of the changes in baseline length. The time series of baseline lengths between stations Kerguelen and Diego Garcia; and Tidbinbilla and Yaragadee shows an increasing trend and the Hurst exponent is close to 0.5. This may suggest that the diffuse plate boundary zone and the Australian plate are moving away from the Antarctican plate. The baseline length between Kerguelen and Davis, Casey, Maitri and Seychelles is less than 0.5, which indicates that the baseline length of these stations shows a decreasing trend.     

Exploratory Factor Analysis of Wireline Logs Using a Float-Encoded Genetic Algorithm

In the paper, a novel inversion approach is used for the solution of the problem of factor analysis. The float-encoded genetic algorithm as a global optimization method is implemented to extract factor variables using open-hole logging data. The suggested statistical workflow is used to give a reliable estimate for not only the factors but also the related petrophysical properties in hydrocarbon formations. In the first step, the factor loadings and scores are estimated by Jöreskog’s fast approximate method, which are gradually improved by the genetic algorithm. The forward problem is solved to calculate wireline logs directly from the factor scores. In each generation, the observed and calculated well logs are compared to update the factor population. During the genetic algorithm run, the average fitness of factor populations is maximized to give the best fit between the observed and theoretical data. By using the empirical relation between the first factor and formation shaliness, the shale volume is estimated along the borehole. Permeability as a derived quantity also correlates with the first factor, which allows its determination from an independent source. The estimation results agree well with those of independent deterministic modeling and core measurements. Case studies from Hungary and the USA demonstrate the feasibility of the global optimization based factor analysis, which provides a useful tool for improved reservoir characterization.