“岩岩”并非“积石貌”

李健吾先生的著名散文《雨中登泰山》(见高中语文课本第一册)中写道:“古诗人形容泰山,说‘泰山岩岩’,注解人告诉你:‘岩岩,积石貌。’”笔者认为“岩岩”并非“积石貌”而应为“高峻貌”。“泰山岩岩”出自《诗经》中的《鲁颂·(门心)宫》,原诗是赞美鲁僖公功绩的,文多溢美,以泰山之高大来喻鲁僖公的所谓伟业。金启华先生在其所著的《诗经全译》中,将“岩岩”译为“高大”。

Prediction and critical transition mechanism for granite fracture:Insights from critical slowing down theory

Rock fracture warning is one of the significant challenges in rock mechanics.Many true triaxial and synchronous acoustic emission(AE)tests were conducted on granite samples.The investigation focused on the characteristics of AE signals preceding granite fracture,based on the critical slowing down(CSD)theory.The granite undergoes a transition from the stable phase to the fracture phase and exhibits a clear CSD phenomenon,characterized by a pronounced increase in variance and autocorrelation coefficient.The variance mutation points were found to be more identifiable and suitable as the primary criterion for predicting precursor information related to granite fracture,compared to the autocorrelation coefficient.It is noteworthy to emphasize that the CSD factor holds greater potential in elucidating the underlying mechanisms responsible for the critical transition of granite fracture,in comparison to the AE timing parameters.Furthermore,a novel multi-parameter collaborative prediction method for rock fracture was developed by comprehensively analyzing predictive information,including abnormal variation modes and the CSD factor of AE characteristic parameters.This method enhances the understanding and prediction of rock fracture-related geohazards.

Research on anisotropy of shale oil reservoir based on rock physics model

Rock physics modeling is implemented for shales in the Luojia area of the Zhanhua topographic depression. In the rock physics model, the clay lamination parameter is introduced into the Backus averaging theory for the description of anisotropy related to the preferred alignment of clay particles, and the Chapman multi-scale fracture theory is used to calculate anisotropy relating to the fracture system. In accordance with geological features of shales in the study area, horizontal fractures are regarded as the dominant factor in the prediction of fracture density and anisotropy parameters for the inversion scheme. Results indicate that the horizontal fracture density obtained has good agreement with horizontal permeability measured from cores, and thus confirms the applicability of the proposed rock physics model and inversion method. Fracture density can thus be regarded as an indicator of reservoir permeability. In addition, the anisotropy parameter of the P-wave is higher than that of the S-wave due to the presence of horizontal fractures. Fracture density has an obvious positive correlation with P-wave anisotropy, and the clay content shows a positive correlation with S-wave anisotropy, which fully shows that fracture density has a negative correlation with clay and quartz contents and a positive relation with carbonate contents.

Effect of pore structure on seismic rock-physics characteristics of dense carbonates

The Ordovician carbonate rocks of the Yingshan formation in the Tarim Basin have a complex pore structure owing to diagenetic and secondary structures. Seismic elastic parameters（e.g., wave velocity） depend on porosity and pore structure. We estimated the average specific surface, average pore-throat radius, pore roundness, and average aspect ratio of carbonate rocks from the Tazhong area. High P-wave velocity samples have small average specific surface, small average pore-throat radius, and large average aspect ratio. Differences in the pore structure of dense carbonate samples lead to fluid-related velocity variability. However, the relation between velocity dispersion and average specific surface, or the average aspect ratio, is not linear. For large or small average specific surface, the pore structure of the rock samples becomes uniform, which weakens squirt fl ow and minimizes the residuals of ultrasonic data and predictions with the Gassmann equation. When rigid dissolved（casting mold） pores coexist with less rigid microcracks, there are significant P-wave velocity differences between measurements and predictions.

Aeromagnetic data and geological structure of continental China:A review

We review the latest aeromagnetic geological data of continental China. We discuss the latest achievements in geological mapping and the newly detected features based on aeromagnetic data. Using aeromagnetic data collected for more than 50 years, a series of 1：5000000 and 1：1000000 aeromagnetic maps of continental China were compiled using state-of-the-art digital technology, and data processing and transformation. Guided by plate tectonics and continental dynamics, rock physical properties, and magnetic anomalies, we compiled maps of the depth of the magnetic basement of continental China and the major geotectonic units, and presented newly detected geological structures based on the aeromagnefic data.

Analyzing the mid-low porosity sandstone dry frame in central Sichuan based on effective medium theory

Tight gas sandstone reservoirs in Guang＇an are characterized by wide distribution and low abundance. Sandstone samples from this area usually have low porosity and poor connectivity. We analyze the observed velocity data of tight sandstone samples with the Mori- Tanaka model, and give the sandstone framework physical model in this area based on theory and experiment analysis. The matrix modulus was obtained by an empirical relationship and then the experiment data were compared with the values predicted by the Mori-Tanaka model with different pore shapes. The results revealed that the experiment data were close to the model with low pore aspect ratio. Considering the matrix modulus and pore shape variation, we find that, under the condition of small mineral composition change, the effective pore aspect ratio of these samples increased with porosity evidently.

Application of Prediction Techniques in Carbonate Karst Reservoir in Tarim Basin

Carbonate karst reservoir is the emphases of Tarim＇s carbonate exploration. However, it is buried at a large depth, which results in Weak seismic reflection signal and low S/N ratio. In addition, the karst reservoir contains great heterogeneity, so reservoir prediction is very difficult. Through many years of research and exploration, we have established a suite of comprehensive evaluation technology for carbonate karst reservoir using geophysical characteristics and a geological concept model, including a technique for reconstructing the paleogeomorphology of buried hills based on a sequence framework, seismic description of the karst reservoir, and strain variant analysis for fracture estimation. The evaluation technology has been successfully applied in the Tabei and Tazhong areas, and commercial production of oil and gas has been achieved. We show the application of this technology in the Lunguxi area in North Tarim in this paper.

A rock-physical modeling method for carbonate reservoirs at seismic scale

Strong heterogeneity and complex pore systems of carbonate reservoir rock make its rock physics model building and fluid substitution difficult and complex. However, rock physics models connect reservoir parameters with seismic parameters and fluid substitution is the most effective tool for reservoir prediction and quantitative characterization. On the basis of analyzing complex carbonate reservoir pore structures and heterogeneity at seismic scale, we use the gridding method to divide carbonate rock into homogeneous blocks with independent rock parameters and calculate the elastic moduli of dry rock units step by step using different rock physics models based on pore origin and structural feature. Then, the elastic moduli of rocks saturated with different fluids are obtained using fluid substitution based on different pore connectivity. Based on the calculated elastic moduli of rock units, the Hashin-Shtrikman-Walpole elastic boundary theory is adopted to calculate the carbonate elastic parameters at seismic scale. The calculation and analysis of carbonate models with different combinations of pore types demonstrate the effects of pore type on rock elastic parameters. The simulated result is consistent with our knowledge of real data.

Sandstone Diagenesis and Porosity Evolution of Paleogene in Huimin Depression

The sandstones of the paleogene in the Huimin Depression have undergone numerous diagenetic processes, such as compaction, cementation, dissolution, metasomasis and recrystallization. The first three — compaction, cementation and dissolution — affect reservoir properties most intensively. An average porosity loss due to compaction is 0.78% per 100 meters, slightly higher in the central and southern parts and lower in the northern. Cementation has resulted in the porosity loss from 8% to 20% at a depth below 1,500 meters. Dissolution, which in most cases is the dissolution of feldspar and of carbonate cement, primarily occurs in two depth intervals, from 1,400 to 2,500 meters and 2,700 to 4,000 meters respectively. New porosity is created through dissolution from 4% to 14%. The sandstones experienced stage A and stage B of early diagenesis and stage A of late diagenesis. At present, most of them are experiencing Phase A2 of late diagenesis. The types of pores in the sand reservoir can be grouped into primary and secondary ones.

桂中坳陷下石炭统鹿寨组页岩气研究

桂中坳陷石炭系发现多处油气苗及沥青，具有良好的页岩气储集条件。对桂中坳陷下石炭统鹿寨组(C J )烃源岩进行了研究，认为烃源岩主要为台前盆地相暗色炭质页岩，鹿寨组下部页岩具有良好的生烃条件。烃源岩有机质丰度高，平均TOC含量为3.08°% ，干酪根类型为IIi 型和II2,来源于水生浮游生物和菌藻类母质，烃源岩最高热解峰温度460~560°C ，处于高成熟阶段，现今以生成天然气为主。页岩中脆性矿物含量为35.6% ~75.2%，粘土矿物含量为3.1%~62.1%，页岩硅质来源中生物成因硅质具有较大贡献，鹿寨组底部富有机质页岩段发育的有机质微孔隙和纳米级一微米级微裂缝是页岩气富集的主要空间。桂中坳陷经历多期叠加构造改造，抬升剥蚀对页岩气保存有重要影响，断层和岩浆热液活动对页岩气保存也不利。断层影响相对较小、盖层厚度超过1 8 0 m 的桂中坳陷西北部地区是下石炭统页岩气勘探的有利方向。

Failure characteristics of three-body model composed of rock and coal with different strength and stiffness

Four different types of three-body model composed of rock and coal with different strength and stiffness were established in order to study the failure characteristics of compound model such as roof-coal-floor. Through stress analysis of the element with variable strength and stiffness extracted from the strong-weak interface, the tri-axial compressive strength of the weak body and strong body near the interface as well as the areas away from the contact surface was found. Then, on the basis of three-dimensional fast Lagrangian method of continua and strain softening constitutive model composed of Coulomb-Mohr shear failure with tensile cut-off, stress and strain relationship of the four three-body combined models were analyzed under different confining pressures by numerical simulation. Finally, the different features of local shear zones and plastic failure areas of the four different models and their development trend with increasing confining pressure were discussed. The results show that additional stresses are derived due to the lateral deformation constraints near the strong-weak interface area, which results in the strength increasing in weak body and strength decreasing in strong body. The weakly consolidated soft rock and coal cementation exhibit significant strain softening behavior and bear compound tension-shear failure under uni-axial compression. With the increase of confining pressure, the tensile failure disappears from the model, and the failure type of composed model changes to local shear failure with different number of shearing bands and plastic failure zones. This work shows important guiding significance for the mechanism study of seismic, rock burst, and coal bump.

Anisotropy rock physics model for the Longmaxi shale gas reservoir,Sichuan Basin,China

The preferred orientation of clay minerals dominates the intrinsic anisotropy of shale. We introduce the clay lamination （CL） parameter to the Backus averaging method to describe the intrinsic shale anisotropy induced by the alignment of clay minerals. Then, we perform the inversion of CL and the Thomsen anisotropy parameters. The direct measurement of anisotropy is difficult because of the inability to measure the acoustic velocity in the vertical direction in boreholes and instrument limitations. By introducing the parameter CL, the inversion method provides reasonable estimates of the elastic anisotropy in the Longmaxi shale. The clay content is weakly correlated with the CL parameter. Moreover, the parameter CL is abnormally high at the bottom of the Longmaxi and Wufeng Formations, which are the target reservoirs. Finally, we construct rock physics templates to interpret well logging and reservoir properties.

也谈“坳”“拗”,“心”“芯”的使用——兼与刘全稳先生商榷

针对《西南石油大学学报》刘全稳副主编发表在该刊2007年第6期上的一篇关于"坳""拗""心""芯"使用方法说明的文章,对石油科技文献中经常出现的"坳陷""拗陷""取芯""取心""岩芯""岩心"的合理使用提出了自己的观点和看法。进而指出:"坳陷""拗陷"两者略有差别,不能一概将所有"坳陷"均改为"拗陷";同时,刘先生认为的"取心""岩心"都必须改为"取芯""岩芯"也是不尽合理的,"取芯"与"取心","岩芯"与"岩心"都可使用,且意思完全一样。

Petrogenesis of skarn in Shizhuyuan W-polymetallic deposit, southern Hunan,China:Constraints from petrology,mineralogy and geochemistry

Skarn is the main altered rock type and is of great importance to mineralization and ore-prospecting in the Shizhuyuan area of Hunan province, China. Its features of petrography, mineralogy and geochemistry were studied systematically. The results show that the skarn mainly consists of garnet skarn, secondary wollastonite-garnet skarn, tremolite-clinozoisite skarn, and few wolframine garnet skarn, idocrase-garnet skarn and wollastonite skarn with granoblastic texture, granular sheet crystalloblastic texture, and massive structure, disseminated structure, mesh-vein structure, comb structure, and banded structure. And, it is mainly composed of garnet, fluorite, chlorite, hornblende, epidote, tremolite, plagioclase, biotite, muscovite, plagioclase, quartz, idocrase, and calcite and so on. The chemical components mainly include SiO2, Al2O3, Fe2O3, MgO and CaO, and the trace elements and REEs consist of Li, Be, V, Co, Zn, Ga, Rb, Sr, Y, Ce, Nd, Pb and Bi, etc. And, the obvious fractionation exists between LREE and HREE, and it shows typical features of Nanling ore-forming granite for W？Sn polymetallic deposit. Skarn is derived from the sedimentary rock, such as limestone, mudstone, argillaceous rock, and few pelitic strips. It is affected by both Shetianqiao formation strata and Qianlishan granite during the diagenesis, indicating a strong reduction environment. The occurrence of skarn, whose mutation site is favorable to the mineralization enrichment, is closely related to the mineralization and prospecting.

Strength characteristics of dry and saturated rock at different strain rates

The strength of rock materials is largely affected by water and loading conditions, but there are few studies on mechanical properties of saturated rocks at high strain rates. Through compressive tests on dry and saturated sandstone specimens, it was found that the dynamic compressive strength of both dry and saturated sandstone specimens increased with the increase of strain rates. The saturated rock specimens showed stronger rate dependence than the dry ones. The water affecting factor （WAF）, as the ratio of the strength under dry state to that under saturated state, was introduced to describe the influence of water on the compressive strength at different strain rates. The WAF under static load was close to 1.38, and decreased with the increase of strain rate. When the strain rate reached 190 s^-1, the WAF reduced to 0.98. It indicates that the compressive strength of saturated specimens can be higher than that of dry ones when the strain rate is high enough. Furthermore, the dual effects of water and strain rate on the strength of rock were discussed based on sliding crack model, which provided a good explanation for the experimental results.

Rock skeleton models and seismic porosity inversion

By substituting rock skeleton modulus expressions into Gassmann approximate fluid equation, we obtain a seismic porosity inversion equation. However, conventional rock skeleton models and their expressions are quite different from each other, resuling in different seismic porosity inversion equations, potentially leading to difficulties in correctly applying them and evaluating their results. In response to this, a uniform relation with two adjusting parameters suitable for all rock skeleton models is established from an analysis and comparison of various conventional rock skeleton models and their expressions including the Eshelby-Walsh, Pride, Geertsma, Nur, Keys-Xu, and Krief models. By giving the two adjusting parameters specific values, different rock skeleton models with specific physical characteristics can be generated. This allows us to select the most appropriate rock skeleton model based on geological and geophysical conditions, and to develop more wise seismic porosity inversion. As an example of using this method for hydrocarbon prediction and fluid identification, we apply this improved porosity inversion, associated with rock physical data and well log data, to the ZJ basin. Research shows that the existence of an abundant hydrocarbon reservoir is dependent on a moderate porosity range, which means we can use the results of seismic porosity inversion to identify oil reservoirs and dry or water-saturated reservoirs. The seismic inversion results are closely correspond to well log porosity curves in the ZJ area, indicating that the uniform relations and inversion methods proposed in this paper are reliable and effective.

Feasibility of stochastic gradient boosting approach for predicting rockburst damage in burst-prone mines

The database of 254 rockburst events was examined for rockburst damage classification using stochastic gradient boosting （SGB） methods. Five potentially relevant indicators including the stress condition factor, the ground support system capacity, the excavation span, the geological structure and the peak particle velocity of rockburst sites were analyzed. The performance of the model was evaluated using a 10 folds cross-validation （CV） procedure with 80%of original data during modeling, and an external testing set （20%） was employed to validate the prediction performance of the SGB model. Two accuracy measures for multi-class problems were employed： classification accuracy rate and Cohen’s Kappa. The accuracy analysis together with Kappa for the rockburst damage dataset reveals that the SGB model for the prediction of rockburst damage is acceptable.

Computation of elastic properties of 3D digital cores from the Longmaxi shale

The dependence of elastic moduli of shales on the mineralogy and microstructure of shales is important for the prediction of sweet spots and shale gas production. Based on 3D digital images of the microstructure of Longmaxi black shale samples using X-ray CT, we built detailed 3D digital images of cores with porosity properties and mineral contents. Next, we used finite-element （FE） methods to derive the elastic properties of the samples. The FE method can accurately model the shale mineralogy. Particular attention is paid to the derived elastic properties and their dependence on porosity and kerogen. The elastic moduli generally decrease with increasing porosity and kerogen, and there is a critical porosity （0.75） and kerogen content （ca. ≤3%） over which the elastic moduli decrease rapidly and slowly, respectively. The derived elastic moduli of gas- and oil-saturated digital cores differ little probably because of the low porosity （4.5%） of the Longmaxi black shale. Clearly, the numerical experiments demonstrated the feasibility of combining microstructure images of shale samples with elastic moduli calculations to predict shale properties.

DESCRIPTION OF A NEW SPECIES OF THE GENUS Laudakia FROM XIZANG (TIBET)

This paper describes a new species of the rock agamid genus Laudakia . The new species is closely similar to L. tuberculata (Hardwicke et Gray),but the new species differs from the latter in that:①nostril is situated in center of elliptic nasal,directing outwards and rearwards;②one supranasal;③superciliary ridge poorly developed,blunt and not everting upwards;④dorsum and flanks with many small,light colored spots,the large,conic scales scattered on flanks are not situated in those spots.

Calculations of rock matrix modulus based on a linear regression relation

The rock matrix bulk modulus or its inverse, the compressive coefficient, is an important input parameter for fluid substitution by the Biot-Gassmann equation in reservoir prediction. However, it is not easy to accurately estimate the bulk modulus by using conventional methods. In this paper, we present a new linear regression equation for calculating the parameter. In order to get this equation, we first derive a simplified Gassmann equation by using a reasonable assumption in which the compressive coefficient of the saturated pore fluid is much greater than the rock matrix, and, second, we use the Eshelby- Walsh relation to replace the equivalent modulus of a dry rock in the Gassmann equation. Results from the rock physics analysis of rock sample from a carbonate area show that rock matrix compressive coefficients calculated with water-saturated and dry rock samples using the linear regression method are very close （their error is less than 1%）. This means the new method is accurate and reliable.