塔里木盆地南缘8.5 ka以来的黄土岩石磁学特征及其磁化率变化机制

黄土岩石磁学参数是古气候研究中的重要指标,其中磁化率应用最为广泛,并在黄土高原地区取得重大进展,其受控于成壤作用的变化机制也被普遍接受.然而在黄土高原外缘的新疆地区,磁化率的变化机制仍不明确,导致磁化率的古气候意义在该区存在较大争议.本文选取塔里木盆地南缘具有精确年代控制的典型黄土剖面（羊场剖面）开展岩石磁学和高分辨率磁化率研究,利用交叉小波分析方法并结合剖面粒度、矿物及元素特征对该地区磁化率变化机制进行初步探讨.结果显示,羊场剖面的岩石磁学性质主要由粗颗粒软磁性矿物所控制,同时也表现出一定的顺磁性特征.根据载磁矿物和磁化率变化特征可将剖面进一步划分为两个阶段：阶段Ⅰ（8.5~2.5 ka）,载磁矿物以亚铁磁性的磁铁矿为主,磁化率值整体较高;阶段Ⅱ（2.5~0.2 ka）,亚铁磁性矿物依然占据主导地位,但硬磁性矿物和以黄铁矿为代表的顺磁性矿物相对增多,磁化率值显著降低.相关性研究和交叉小波分析表明：阶段Ⅰ磁化率与粗颗粒组分的变化具有一致性,符合"风速论"模式;阶段Ⅱ磁化率不仅与粗颗粒组分具有明显的正相关关系,而且与指示成壤作用强度的频率磁化率百分含量呈现出显著的负相关关系,暗示了阶段Ⅱ的磁化率变化可能受到"风速论"和"还原性成壤"模式的共同影响.本文拓宽了对新疆地区黄土岩石磁学特征及其磁化率变化机制的深入理解,也为利用磁化率恢复新疆及中亚地区全新世以来的古气候变化历史提供了新的线索.     

单裂缝不同传播方向Krauklis波振幅变化的实验研究

裂缝广泛分布于地壳岩石中并具有多尺度的特点,研究地震波在裂缝介质中的传播规律具有重要的科学意义和应用价值.含流体的裂缝中传播的Krauklis波（也曾被称为慢流体波、裂缝波、斯通利导波等）近年来逐渐成为地球物理学界的研究热点之一.本研究利用地震物理模型技术研究横波在单裂缝样品中传播时随传播方向变化的响应特征,在实验室用环氧树脂为基质构建模型样品,一个为不含裂缝的参考样品,一个为含直径30 mm、厚度0.2 mm的薄币状单裂缝的模型,利用0.5 MHz换能器测试了横波在模型中传播时透射信号的特征.实验结果显示在单裂缝样品中透射横波发生了波形的畸变和振幅降低,振幅谱显示透射信号的主频明显降低,在初至横波之后为具有一定能量的尾波.对尾波进行时频分析,发现在方位角较小时尾波能量较弱而方位角较大时尾波能量较强,尾波中信号的主频在0.1 MHz左右,与前人理论和实验研究结果一致.

     

A triple porosity scheme for fluid/solid substitution: theory and experiment

Quantifying the effects of pore-filling materials on elastic properties of porous rocks is of considerable interest in geophysical practice. For rocks saturated with fluids, the Gassmann equation is proved effective in estimating the exact change in seismic velocity or rock moduli upon the changes in properties of pore infill. For solid substance or viscoelastic materials, however, the Gassmann theory is not applicable as the rigidity of the pore fill (either elastic or viscoelastic) prevents pressure communication in the pore space, which is a key assumption of the Gassmann equation. In this paper, we explored the elastic properties of a sandstone sample saturated with fluid and solid substance under different confining pressures. This sandstone sample is saturated with octadecane, which is a hydrocarbon with a melting point of 28°C, making it convenient to use in the lab in both solid and fluid forms. Ultrasonically measured velocities of the dry rock exhibit strong pressure dependency, which is largely reduced for the filling of solid octadecane. Predictions by the Gassmann theory for the elastic moduli of the sandstone saturated with liquid octadecane are consistent with ultrasonic measurements, but underestimate the elastic moduli of the sandstone saturated with solid octadecane. Our analysis shows that the difference between the elastic moduli of the dry and solid-octadecane-saturated sandstone is controlled by the squirt flow between stiff, compliant, and the so-called intermediate pores (with an aspect ratio larger than that of compliant pore but much less than that of stiff pores). Therefore, we developed a triple porosity model to quantify the combined squirt flow effects of compliant and intermediate pores saturated with solid or viscoelastic infill. Full saturation of remaining stiff pores with solid or viscoelastic materials is then considered by the lower embedded bound theory. The proposed model gave a reasonable fit to the ultrasonic measurements of the elastic moduli of the sandstone saturated with liquid or solid octadecane. Comparison of the predictions by the new model to other solid substitution schemes implied that accounting for the combined effects of compliant and intermediate pores is necessary to explain the solid squirt effects.     

Acoustic and petrophysical properties of mechanically compacted overconsolidated sands: Part 2 – Rock physics modelling and applications

Part one of this paper reported results from experimental compaction measurements of unconsolidated natural sand samples with different mineralogical compositions and textures. The experimental setup was designed with several cycles of stress loading and unloading applied to the samples. The setup was aimed to simulate a stress condition where sediments underwent episodes of compaction, uplift and erosion. P-wave and S-wave velocities and corresponding petrophysical (porosity and density) properties were reported. In this second part of the paper, rock physics modelling utilizing existing rock physics models to evaluate the model validity for measured data from part one were presented. The results show that a friable sand model, which was established for normally compacted sediments is also capable of describing overconsolidated sediments. The velocity–porosity data plotted along the friable sand lines not only describe sorting deterioration, as has been traditionally explained by other studies, but also variations in pre-consolidation stress or degree of stress release. The deviation of the overconsolidated sands away from the normal compaction trend on the V_P/V_S and acoustic impedance space shows that various stress paths can be predicted on this domain when utilizing rock physics templates. Fluid saturation sensitivity is found to be lower in overconsolidated sands compared to normally consolidated sands. The sensitivity decreases with increasing pre-consolidation stress. This means detectability for four-dimensional fluid saturation changes can be affected if sediments were pre-stressed and unloaded. Well log data from the Barents Sea show similar patterns to the experimental sand data. The findings allow the development of better rock physics diagnostics of unloaded sediments, and the understanding of expected 4D seismic response during time-lapse seismic monitoring of uplifted basins. The studied outcomes also reveal an insight into the friable sand model that its diagnostic value is not only for describing sorting microtextures, but also pre-consolidation stress history. The outcome extends the model application for pre-consolidation stress estimation, for any unconsolidated sands experiencing similar unloading stress conditions to this study.     

CO2 saturation estimates at Sleipner (North Sea) from seismic tomography and rock physics inversion

CO₂ saturations are estimated at Sleipner using a two-step imaging workflow. The workflow combines seismic tomography (full-waveform inversion) and rock physics inversion and is applied to a two-dimensional seismic line located near the injection point at Sleipner. We use baseline data (1994 vintage, before CO₂ injection) and monitor data that was acquired after 12 years of CO₂ injection (2008 vintage). P-wave velocity models are generated using the Full waveform inversion technology and then, we invert selected rock physics parameters using an rock physics inversion methodology. Full waveform inversion provides high-resolution P-wave velocity models both for baseline and monitor data. The physical relations between rock physics properties and acoustic wave velocities in the Utsira unconsolidated sandstone (reservoir formation) are defined using a dynamic rock physics model based on well-known Biot–Gassmann theories. For data prior to injection, rock frame properties (porosity, bulk and shear dry moduli) are estimated using rock physics inversion that allows deriving physically consistent properties with related uncertainty. We show that the uncertainty related to limited input data (only P-wave velocity) is not an issue because the mean values of parameters are correct. These rock frame properties are then used as a priori constraint in the monitor case. For monitor data, the Full waveform inversion results show nicely resolved thin layers of CO₂–brine saturated sandstones under intra-reservoir shale layers. The CO₂ saturation estimation is carried out by plugging an effective fluid phase in the rock physics model. Calculating the effective fluid bulk modulus of the brine–CO₂ mixture (using Brie equation in our study) is shown to be the key factor to link P-wave velocity to CO₂ saturation. The inversion tests are done with several values of Brie/patchiness exponent and show that the CO₂ saturation estimates are varying between 0.30 and 0.90 depending on the rock physics model and the location in the reservoir. The uncertainty in CO₂ saturation estimation is usually lower than 0.20. When the patchiness exponent is considered as unknown, the inversion is less constrained and we end up with values of exponent varying between 5 and 20 and up to 33 in specific reservoir areas. These estimations tend to show that the CO₂–brine mixing is between uniform and patchy mixing and variable throughout the reservoir.     

A semi-empirical approach to model pressure dependence of elastic moduli in granular media accounting for variations of coordination-number and Poisson-ratio

The effective medium theory based on the Hertz–Mindlin contact law is the most popular theory to relate dynamic elastic moduli (or elastic velocities) and confining pressure in dry granular media. However, many experimental results proved that the effective medium theory predicts pressure trends lower than experimental ones and over-predicts the shear modulus. To mitigate these mispredictions, several evolutions of the effective medium theory have been presented in the literature. Among these, the model named modified grain contact theory is an empirical approach in which three parametric curves are included in the effective medium theory model. Fitting the parameters of these curves permits to adjust the pressure trends of the Poisson ratio and the bulk modulus. In this paper, we present two variations of the modified grain contact theory model. First, we propose a minor modification in the fitting function for the porosity dependence of the calibration parameters that accounts for non-linearity in the vicinity of the critical porosity. Second, we propose a major modification that reduces the three-step modified grain contact theory model to a two-step model, by skipping the calibration parameter–porosity fit in the model and directly modelling the calibration parameter–pressure relation. In addition to an increased simplicity (the fitting parameters are reduced from 10 to 6), avoiding the porosity fit permits us to apply the model to laboratory data that are not provided with accurate porosity measurements. For this second model, we also estimate the uncertainty of the fitting parameters and the elastic velocities. We tested this model on dry core measurements from literature and we verified that it returns elastic velocity trends as good as the original modified grain contact theory model with a reduced number of fitting parameters. Possible developments of the new model to add predictive power are also discussed.     

华北盆地强震孕育的动力学机制研究

强震活跃区通常出现在地壳变形剧烈的地区.但是,华北盆地地壳变形速率小,而破坏性强震却十分频繁.文中利用有限元方法对华北盆地发震断层下方可能存在的速度异常体进行了力学分析,计算结果显示：无论这个异常区是低速异常还是高速异常,在水平构造力的作用下（无论是挤压还是引张力作用）,应力及应变能密度都在空间某些区域积聚.根据本文的工作及前人的研究结果,认为华北盆地强震孕育的物理机制如下：华北克拉通在减薄、破坏过程中,浸入的岩浆或上涌的地幔物质经过漫长的地质时期演化,成为低速异常体.当断裂带下面存在这种低速异常体时,水平向构造力的作用,会使异常体上方应力集中、应变能积聚.当异常体上方的应力水平达到某种限度,介质发生破坏.在破坏过程中,异常体释放的能量,使断裂带上的应变能集聚速度加快;另一方面,异常体的破坏使断层的闭锁面积减小.随着外部构造应力的持续加载,当发震断层上的应力及能量积累达到一定极限,断层产生错动,发生强震.同时,文中结果还显示,断层在从闭锁、积累能量到解锁、发生强震的过程中,其下方的低速异常体起着重要的推进作用.     

Thermal properties as controls on rock surface temperature maxima,and possible implications for rock weathering

Four rock types (basalt, sandstone, granite, and chalk) are examined with respect to the maximum surface temperatures which they experience when subjected to similar conditions of exposure. Rock temperature measurements are reported for an urban environment and for two experimental situations in which an infrared lamp is used to simulate heating under cold and hot conditions. Differences in rock temperatures are discussed with reference to thermal rock properties (albedo, specific heat capacity, and thermal conductivity). Some natural situations are suggested in which thermal rock properties could conceivably play a role in determining the extent to which rocks would be affected by particular weathering processes.     

Origin of layering in cumulate gabbros in the Oman ophiolite: Insights from magnetic susceptibility measurements in the Wadi Sadm area

Abstract   Magnetic susceptibility and the anisotropy of magnetic susceptibility were measured on an 800-cm-thick succession of cumulate gabbro in the Sadm area of the Oman ophiolite. The section contained three distinct cumulate units. The susceptibility tends to decrease upward in each from a melanocratic layer (several tens of centimeters thick) to a leucocratic layer (a few meters thick). The susceptibility decreases in accordance with the decreasing number of magnetite grains, which are the alteration product mainly of olivine minerals. This suggests the cyclic downward accumulation of olivine in the cumulate gabbro. The apparent strain deduced from the patterns of magnetic and grain fabrics was the result mostly of simple shear, so that the layering of gabbro is understood to be formed primarily by a crystal cumulus process followed by simple shear deformation.     

Isolating quartz-dominated OSL signal of rock slice by using pulsed stimulation: Implications for dating burial age of cobbles

Although the OSL signal sourced from quartz is expected to be more stable and bleached more rapidly than the IRSL signal sourced from feldspar in a general sense, the former is much less investigated than the latter for rock surface related luminescence dating. It is mainly due to the difficulty in isolating quartz-dominated OSL signals from rock slices, and the low sensitivity and non-fast component dominated OSL signals for quartz in most of rocks. In merit of the sub-conduction band transition of trapped electrons of feldspar under both green light stimulation and IR stimulation at elevated temperature (Jain and Ankjӕrgaard, 2011), it is expected that the contribution of feldspars to the green light stimulated luminescence could be substantially reduced by a prior IR stimulation at elevated temperature. Meanwhile, more fast-component dominated quartz OSL signal (if there is any) could be obtained by green light stimulation (Bailey et al., 2011). Therefore, in this study, we investigated the feasibility of using the post-IR pulsed green stimulation, which is performed at 25 °C (PGLSL₂₅) following two IR stimulations (IR₅₀IR₂₂₅), to isolate the quartz-dominated OSL signals from rock slices of granite cobbles for burial dating. The decay characteristics of stimulation curve, characteristic saturation dose and thermal stability of the pIR₅₀IR₂₂₅ PGLSL₂₅ signal suggest that this signal is quartz-dominated, but still not fast-component dominated. We tentatively validated the equivalent dose (D_e) measurement protocol by dose recovery experiment. The luminescence-depth profiles show that the bleaching depths of pIR₅₀IR₂₂₅ PGLSL₂₅ signals are slightly smaller than or close to that of the IR₅₀ signals, while they are much shallower for the pIR₅₀IR₂₂₅ signals. The pIR₅₀IR₂₂₅ PGLSL₂₅ procedure enables multiple D_e values determined from luminescence signals, with different bleachabilities and stabilities, from both feldspars and quartz in one measurement, which is of potential for buried age dating of cobbles.