Geomagnetic field reversal in Cretaceous Tirupati sandstone formation from India

Paleomagnetic study of core samples drilled through 456 cm vertical thickness of the Tirupati sandstones formation of Lower to Middle Cretaceous age from Godavary Valley in India is reported. The results give a strong evidence of a geomagnetic field reversal during the period of magnetization of this formation. This field reversal is opposite in sense to that observed in Deccan traps of Cretaceous to Eocene age by several workers, thus suggesting that there have been at least two field reversals during this period.     

Paleomagnetic field intensity derived from meteorite magnetization

The paleomagnetic field intensity is estimated with the aid of the Koenigsberger-Thellier method for four ordinary chondrites and one carbonaceous chondrite by assuming that the stable NRM component of these meteorites is attributable to the TRM acquired in a low-temperature range (lower than about 400°C) during their extremely-slow cooling process. The results are summarized in Table IV, where the paleomagnetic field intensity ranges from 0.10 to 0.97 Oe.Possible effects of the extremely-slow cooling rate of meteorites and the secondary TRM acquisition of the surface fusion crust upon the original NRM of the meteorite interior are discussed.     

砂岩核磁共振响应模拟及受限扩散

本文运用随机游走方法模拟了砂岩储层中流体的核磁共振(NMR)响应及其受限扩散现象.通过改变数字岩心的分辨率模拟生成不同孔隙尺寸的砂岩,研究了不同孔隙尺寸砂岩饱含水时流体扩散系数随扩散时间的变化关系,同时模拟了砂岩饱和单相流体和两相流体的NMR响应;研究了流体的受限扩散系数与横向弛豫时间T2的关系,分析了表面弛豫率和胶结指数对润湿相流体受限扩散系数线位置的影响,并将其用于解释砂岩储层的D-T2分布.结果表明:孔隙流体的扩散系数会随扩散时间的增加而逐渐减小并趋于定值.随着岩石孔隙尺寸的减小,受限扩散现象越明显,受限扩散对岩石NMR响应的影响也越大.润湿相流体受限扩散系数线的位置受岩石胶结指数和表面弛豫率的影响较大.由于润湿相流体扩散系数减小,导致D-T2分布中润湿相流体信号偏离其自由扩散系数线,需要利用流体的受限扩散系数线准确识别D-T2分布中的润湿相流体.     

核磁共振共轭梯度解谱方法研究

核磁共振（NMR）解谱技术是核磁共振资料应用的基础和关键.本文将核磁共振解谱的混定线性反演问题转化为求目标函数极小值的最优化问题,然后利用共轭梯度算法具有二次终止性、收敛速度快的特点解决上述最优化问题.将该方法应用于无噪声和不同信噪比的理论数据解谱以及岩心NMR实验解谱并与真谱及实验室国外软件解谱结果对比表明：在信噪比SNR≥5时解谱结果和真谱符合得很好,谱线光滑连续;随着信噪比的降低对初始点的要求随之提高;两块岩心解谱结果与实验室结果符合得很好,利用解谱结果计算的核磁孔隙度与实验室氦孔隙度绝对误差分别为0.78%和0.57%.因此,本文方法有效、实用,具有较强的抗噪能力,对初始点的要求不高,能够应用于生产和科研实践中.     

用GRACE卫星跟踪数据反演地球重力场

利用141天GRACE卫星观测资料,包括K波段、星载加速度和卫星轨道数据,反演了80阶地球重力场模型IGGGRACE01S,该模型在半波长为500km的空间分辨率上,确定大地水准面的精度约为0012m,中长波(<80阶)精度优于重力卫星发射以前研制的重力场模型. 与EIGEN_GRACE02S、EIGEN_CHAMP03S和EGM96模型的位系数相比,该模型系数最接近于EIGEN_GRACE02S,与另两个模型差异较大. 比较几种模型确定的全球重力异常和大地水准面起伏,结果发现IGGGRACE01S与EIGEN_GRACE02S模型的计算结果比较接近,与EGM96模型结果差异较大,差别较大地区主要在南极等地区. 对于中国大陆,比较IGGGRACE01S模型（前72阶）计算的重力异常和NIMA重力异常数据（25°×25°网格）,两者之间的标准偏差为48mGal.     

Current horizontal strain field in Chinese mainland derived from GPS data

Introduction In the years when the reliable data could not be obtained and in the analysis of strain property and magnitude in history, the intensity, property and activity pattern of strain field were mainly inferred on the bases of geometric characters of surface traces and behaviors (especially the faults) as well as the characteristics of petrology (XIE, et al, 1993; Molnar, Tapponnier, 1975, 1977; Tapponnier, Molnar, 1977; FU, et al, 2000). However, they are the averaged results accumu…     

核磁共振技术在特低渗砂岩微观孔隙结构评价中的应用

特低渗透砂岩储层微观孔隙结构变化复杂,且随机性强,利用核磁共振技术对取自鄂尔多斯盆地延长组的特低渗砂岩样品进行了测试分析.结果表明,实验样品的T2截止值小,T2谱分布表现出双峰态;可动流体参数差异较大,主要受孔隙(尤其是次生孔隙)的发育和连通程度、微裂缝的发育程度、粘土矿物含量及其赋存形式、孔隙的比表面大小及重结晶程度...     

Recovery of the gravity field from GOCE data by using the invariants of gradient tensor

The gravity field models GUCAS_EGM and GUCAS_EGM_DL are established from GOCE data (GOCE Level 2 Products from Nov. 1 to Dec. 31, 2009) based on the method of the invariants of the gravity gradient tensor, where GUCAS_EGM is derived after GOCE gravity gradient data are filtered with FIR, and GUCAS_EGM_DL is computed with an additional Durbin-Levison arithmetic apart from FIR. Since this method, different from current programs dealing with GOCE data, is introduced for the first time, some new problems are required to be discussed in advance; for example, how to filter GOCE gravity gradient data, how to compute the invariants of the gradient tensor, and how to deal with the pole gap and so on. In addition, by comparing our models with ones recommended by ESA, it can be seen that the variations of GUCAS_EGM and the models recommended by ESA to EGM08 are almost equivalent, and the variation of GUCAS_EGM_DL to EGM08 is obviously less than ones of the recommended models.     

Pore size distributions and hydraulic conductivities of rocks derived from Magnetic Resonance Sounding relaxation data using multi-exponential decay time inversion

In hydrogeology there is a variety of empirical formulae available for determination of hydraulic conductivity of porous media, all based on the analysis of grain size distributions of aquifer materials. Sensitivity of NMR measurements to pore sizes makes it a good indicator of hydraulic conductivity. Analogous to laboratory NMR, Magnetic Resonance Sounding (MRS) relaxation data are of a multi-exponential (ME) nature due to the distribution of different pore sizes in an investigated rock layer. ME relaxation behaviour will also arise due to the superposition of NMR signals which originate from different layers. It has been shown, that both kinds of ME behaviour coexist in MRS and can principally be separated by ME inversion of the field data. Only a few publications exist that have proposed approaches to qualitatively and quantitatively estimate petrophysical parameters such as the hydraulic conductivity from MRS measurements, i.e. MRS porosity and decay times. The so far used relations for the estimation of hydraulic conductivity in hydrogeology and NMR experiments are compared and discussed with respect to their applicability in MRS. Taking into account results from a variety of laboratory NMR and MRS experiments mean rock specific calibration factors are introduced for a data-base-calibrated estimation of hydraulic conductivity when no on-site calibration of MRS is available. Field data have been analysed using conventional and ME inversion using such mean calibration values. The results for conventional and ME inversion agree with estimates obtained from well core analysis for shallow depths but are significantly improved using a ME inversion approach for greater depths.     

Pore size distributions and hydraulic conductivities of rocks derived from Magnetic Resonance Sounding relaxation data using multi-exponential decay time inversion

In hydrogeology there is a variety of empirical formulae available for determination of hydraulic conductivity of porous media, all based on the analysis of grain size distributions of aquifer materials. Sensitivity of NMR measurements to pore sizes makes it a good indicator of hydraulic conductivity. Analogous to laboratory NMR, Magnetic Resonance Sounding (MRS) relaxation data are of a multi-exponential (ME) nature due to the distribution of different pore sizes in an investigated rock layer. ME relaxation behaviour will also arise due to the superposition of NMR signals which originate from different layers. It has been shown, that both kinds of ME behaviour coexist in MRS and can principally be separated by ME inversion of the field data. Only a few publications exist that have proposed approaches to qualitatively and quantitatively estimate petrophysical parameters such as the hydraulic conductivity from MRS measurements, i.e. MRS porosity and decay times. The so far used relations for the estimation of hydraulic conductivity in hydrogeology and NMR experiments are compared and discussed with respect to their applicability in MRS. Taking into account results from a variety of laboratory NMR and MRS experiments mean rock specific calibration factors are introduced for a data-base-calibrated estimation of hydraulic conductivity when no on-site calibration of MRS is available. Field data have been analysed using conventional and ME inversion using such mean calibration values. The results for conventional and ME inversion agree with estimates obtained from well core analysis for shallow depths but are significantly improved using a ME inversion approach for greater depths.     

Controls on the geometry of a Holocene crater row: a field study from southwest Iceland

The Reykjanes Peninsula rift zone in southwest Iceland is a highly oblique segment of the Mid-Atlantic ridge system which accommodates NW–SE extension during rifting episodes that consist of eruptions and normal faulting, and E–W left-lateral shear strain along strike-slip faults during longer amagmatic periods. Dominant tectonic features on the peninsula are a series of generally NE-striking, sub-parallel eruptive fissures and normal faults, and a cross-cutting zone of N–S striking, right-lateral strike-slip faults. The last series of rifting episodes ended in 1227, and a proposed 1,000 year cyclicity predicts the start of a new series of eruptions within the next 200 years. In order to more accurately characterize the nature of eruptions on the Reykjanes Peninsula, we present a new, spatially accurate map of the ∼2,350 year old Sundhnúkur crater row in the western part of the peninsula, which was examined in detail in order to determine the structural controls on crater row geometry and to understand the interactions that take place between eruptive fissures and pre-existing geological structures. Volcanism is sometimes influenced by small perturbations in the surroundings such as gravitational loading, topography, changes in crustal properties or the presence of fault zones, but there are few field examples showing how fissures are influenced by these pre-existing structures. We identify 27 fissure segments, ranging in strike from 006° to 053°, with varying spacing and overlap between them. Significant local variability in strike and stepping sense of segments occurs in proximity to topographic highs as well as within zones of faulting that pre-date the crater row. Strike also varies at the northern end of the crater row as it approaches a region of older crust at the rift margin. Our data support numerical and laboratory modeling results which show that local topography, pre-existing fractures and crustal properties influence the path taken by magma on its way through the shallow crust.     

融合GOCE和GRACE卫星数据的无约束重力场模型Tongji-GOGR2019S

本文在法方程层面融合GOCE卫星的V_xx、V_yy、V_zz和V_xz重力梯度分量观测数据和GRACE卫星观测数据,采用直接法解算了220阶次的重力场模型Tongji-GOGR2019S.首先利用ⅡR带通滤波器在5~41 mHz的重力梯度带宽范围内对约24个月的GOCE重力梯度观测方程进行无相移滤波处理,并组成解算220阶次重力场模型的法方程,各梯度分量根据相对于参考模型统计精度进行定权;然后与13.5 a GRACE数据建立的180阶次Tongji-Grace02s重力场模型的法方程进行叠加,解算了220阶次的无约束纯卫星重力场模型Tongji-GOGR2019S.利用EIGEN-6C4重力场模型、GNSS/水准数据、DTU15重力异常数据以及欧洲区域似大地水准面模型EGG2015等数据对Tongji-GOGR2019S模型精度进行全面的检核评定,结果表明：引入GOCE卫星梯度数据后,高于72阶的位系数精度优于Tongji-Grace02s模型,Tongji-GOGR2019S模型的整体精度接近同阶次的DIR-R6等GOCE卫星第6代模型.     

Active blocks and their deformations derived from the velocity field in North China

Based on the horizontal movement velocities and their error estimations of 144GPS stations in North China de-duced from the results obtained from multiple GPS measurements carried out in the period of 1992-2001,the horizontal deformation in the area is stuklied.The movements,entire deformations,local deformations,activeity patterns and intensities along the boundary zones are derived.And then the risk of strong earthquake in the area is estimated.In the research,the horizontal displacement observations can be considered as the sum of three parts.The first part is the entire motion following Eurasian plate,which can be derived from NUVAL-IA model;the second part is the relative motions and deformations between the sub-tectonic blocks in the studied area,which can be derived from a set of displacement observations determined by the expanded QUAD method in the paper;the third part is the local deformations and errors in the sub-tectonic blocks,which can be described as the inhomoge-neous strains of the block.The method and results are introduced in detail in the paper.     

Gravity field contribution analysis of GOCE gravitational gradient components

A gravity field model is computed from the four accurate gravitational gradient components of GOCE (Gravity field and steady-state Ocean Circulation Explorer), combined with the analysis of the kinematic orbits, and some moderate constraint (or stabilization) in the polar areas where no observation from GOCE is available due to the orbit geometry. The normal matrix of each component is computed individually in order to study its contribution to the combined solution. The results show that the contribution of V_zz is the largest, with an average value of 32.74% of the total solution; the second and the third largest are V_zz and V_yy, with average values of 28.04% and 26.08%, respectively; the component V_xz contributes 11.81%. Validation with external data shows that each component has its characteristic value and that the information content of the component V_xz is not negligible and should be included for gravity field recovery. The orbit part as derived from high-low satellite-to-satellite tracking (SST-hl) to the GPS contributes mostly to the coefficients below degree and order (d/o) 20, and to non-zonal coefficients from d/o 20 to 80. The mean value of the contribution of the polar stabilization is the smallest with a value of 0.22%, nevertheless it is important. In addition to the contribution analysis in terms of the normal matrices, each individual component of the gradiometer has been combined with SST and polar stabilization, to give a set of single component gravity field models. These partially combined solutions are compared to the fully combined solution in terms of geoid differences. They show that the partially combined solution with V_zz is closest to the complete solution. Even closer is a combination with V_xx and V_yy. In addition to the GOCE-only solution, a GOCE-GRACE (Gravity Recovery And Climate Experiment) combined gravity field model is derived and the information content of GOCE and an available set of normal equations of GRACE are investigated. Results show that, as expected, GRACE dominates the solution below degree 90 and GOCE above degree 140.     

基于GRACE数据的尼泊尔MS8.1地震北向重力梯度变化

利用德克萨斯大学空间研究中心发布的GRACE RL05月重力场模型数据,采用300 km扇形滤波,得到2015年尼泊尔M_S8.1地震北向重力梯度的时空分布,之后利用最小二乘拟合方法分析尼泊尔及其邻近区域北向重力梯度的长期变化趋势及研究区内6个特征点的北向重力梯度月变化时间序列,并结合黏弹性分层位错理论对GRACE检测尼泊尔M_S8.1地震北向重力梯度同震变化的可能性进行分析. 研究结果表明:在尼泊尔M_S8.1地震发生前北向重力梯度表现出比较明显的正负异常变化,而该地震即发生在北向重力梯度正负变化的零值线附近;研究区北向重力梯度年变率在印度板块边界及其垂直方向所形成的四象限呈正负相间分布,6个特征点的北向重力梯度在2012年4月出现比较明显的跳变. 由此推断,北向重力梯度的动态变化过程反映了震前区域物质迁移和震后壳幔物质黏滞性调整等问题.      

Comparative study of models for predicting permeability from nuclear magnetic resonance (NMR) logs in two Chinese tight sandstone reservoirs

Based on the analysis of mercury injection capillary pressure (MICP) and nuclear magnetic resonance (NMR) experimental data for core plugs, which were drilled from two Chinese tight sandstone reservoirs, permeability prediction models, such as the classical SDR, Timur-Coates, the Swanson parameter, the Capillary Parachor, the R10 and R35 models, are calibrated to estimating permeabilities from field NMR logs, and the applicabilities of these permeability prediction models are compared. The processing results of several field examples show that the SDR model is unavailable in tight sandstone reservoirs. The Timur-Coates model is effective once the optimal T _2cutoff can be acquired to accurately calculate FFI and BVI from field NMR logs. The Swanson parameter model and the Capillary Parachor model are not always available in tight sandstone reservoirs. The R35 based model cannot effectively work in tight sandstone reservoirs, while the R10 based model is optimal in permeability prediction.     

The intensity of the earth's magnetic field in the past derived from the archaeomagnetic data

a n¶rt;u¶rt; amuu amua uuu anmuaum n a n¶rt;u 8500 m nm um au aam m uuu m naama. ¶rt;u u ¶rt;a ¶rt; u ¶rt; m¶rt; mumu auu nau nuu aau mama [1–3] (Puc. 1.B). nma aau auau (Puc.2) mnaua m m m¶rt;a, ¶rt; au n¶rt; amu u, m n¶rt;u 2- ma m, mmu naa auu nu¶rt; n¶rt;a 1000 m u 350 m, au a aumu, mmm 80%- mmu. am au ¶rt;a u anmu ¶rt; m¶rt; mumu (mum) naa auu u¶rt; ma u ¶rt; m naama[5] (Puc. 3). uma a mam nma aaua naaa auu nu¶rt; n¶rt;a 750 m, 200 m u a — 300 m. (Puc. 4.) nmu u ¶rt;a nma amua an¶rt;u n¶rt;u n u ma m (au) auumu m ¶rt;m m (Puc. 5). aa a m¶rt;m nu naamaum n aam ¶rt;uu mu au, anum¶rt;a m uua nm m¶rt;a.

---

---

Dedicated to Academician Alois Zátopek on His 65th Birthday