Indirect unstructured hex-dominant mesh generation using tetrahedra recombination

Corner-point gridding is widely used in reservoir and basin modeling but generally yields approximations in the representation of geological interfaces. This paper introduces an indirect method to generate a hex-dominant mesh conformal to 3D geological surfaces and well paths suitable for finite-element and control-volume finite-element simulations. By indirect, we mean that the method first generates an unstructured tetrahedral mesh whose tetrahedra are then merged into primitives (hexahedra, prisms, and pyramids). More specifically, we focus on determining the optimal set of primitives that can be recombined from a given tetrahedral mesh. First, we detect in the tetrahedral mesh all the feasible volumetric primitives using a pattern-matching algorithm (Meshkat and Talmor Int. J. Numer. Meth. Eng. 49(1-2), 17–30 2000) that we re-visit and extend with configurations that account for degenerated tetrahedra (slivers). Then, we observe that selecting the optimal set of primitives among the feasible ones can be formalized as a maximum weighted independent set problem (Bomze et al. 1999), known to be \(\mathcal {N}\mathcal {P}\)-Complete. We propose several heuristic optimizations to find a reasonable set of primitives in a practical time. All the tetrahedra of each selected primitive are then merged to build the final unstructured hex-dominant mesh. This method is demonstrated on 3D geological models including a faulted and folded model and a discrete fracture network.     

Enrichment characteristics and risk assessment of Hg in bird feathers from Caohai wetland in Guizhou Province,China

Based on the analysis of the enrichment characteristics of Hg and MeHg in bird feathers from Caohai National Nature Reserve in Guizhou, the risks of Hg pollution to the birds from Caohai wetland have been evaluated. The total Hg content of bird feathers ranges from 40 to 5058 ng/g with an average of 924 ng/g. The content of MeHg is significantly correlated with total Hg (r = 0.68, p < 0.01), and the content are among 0.75 and 113 ng/g. The total Hg content in the birds feathers is significantly dependent on their feeding habits, which is mainly in accordance with the following rule: carnivorous birds > omnivorous birds that are mainly carnivorous > omnivorous birds that are mainly herbivorous. There are also differences in the Hg enrichment ability in the different parts of bird feathers, and the total Hg and MeHg content in the wing feathers are significantly higher than that in the other parts of feathers. The bioaccumulation coefficients of aqueous Hg and MeHg by bird feathers are 0.9 × 10⁴–112.13 × 10⁴ (mean value is 20.47 × 10⁴) and 0.47 × 10⁴–70.4 × 10⁴ (mean value is 9.52 × 10⁴), respectively. Although the whole Hg level in Caohai bird feathers is not too high, the Hg content in some carnivorous birds exceeds over or approaches the abnormal threshold when birds are breeding (5 μg/g), which indicates that the birds in Caohai wetland are faced with some risks of ecological Hg pollution.     

论卡以头组的时代

中国科学院南京地质古生物研究所的有关人员 (1980 )曾做过大量地层古生物工作 ,证实卡以头砂页岩层或卡以头组的时代为三叠纪最早期 ,是飞仙关组底部地层的相变。“论卡以头组”一文认为该组下部的时代为二叠纪 ,上部则属三叠纪。对该文提供的地层古生物资料进行分析后 ,本文作者认为该论点缺乏证据 ,不能成立 ;将陆相的威宁哲觉剖面代替标准地点宣威羊场的过渡相剖面作为卡以头组的标准剖面 ,也是不合适的     

钻孔摩擦灌注桩岩土设计参数取值商榷

桩周岩土极限摩阻力及其桩尖处岩土的容许承载力是公路工程钻(挖)孔摩擦灌注桩单桩竖向承载力计算的重要岩土工程设计参数,其取值是否适当事关工程安危和工程费用的经济合理.对此,工程勘察设计人员不仅应遵照<公路桥涵地基与基础设计规范>有关规定外,还应考虑场地地基岩土的地质时代、场地工程地质条件、岩土工程性质及有关因素合理取值,不可机械地依规范查表确定,这是岩土工程师应具有的最基本技能.     

东南极格罗夫山地区夏季的天气特征

利用中国南极科学考察队 1 998年 1 2月至 1 999年 1月在东南极格罗夫山地区地质考察期间所获得的天气资料 ,对该地区的温度、雪面温度、风速、风向等天气要素进行了分析。结果表明 :夏季格罗夫地区温度日变化和盛行风向与中山站相类似 ,但温度日较差和强风频率比中山站大。当受来自北方暖湿气流影响时 ,常出现降雪天气 ;在东风气流控制下 ,天气以晴为主     

Evaluation of CBM reservoir fracture characteristics by azimuthal AVA inversion

Exploitation of coalbed methane (CBM) is directly affected by the strike and density of fractures, particularly because fractures can act as conduits for water and gas (after desorption). To realize the azimuthal amplitude versus angle (AVA) inversion for a CBM reservoir, we treated reservoirs whose fractures are aligned vertically as being equivalent to a horizontal transverse isotropic medium and simplified the P-P wave reflection coefficient formula to an approximately linear equation. Then, we adopted a modified damping least squares algorithm to solve the azimuthal AVA inversion problem. With edited macro-bins, corrected residual moveout, and normalized amplitude, the azimuthal gathers exhibited unambiguous azimuthal AVA characteristics. The #3 CBM reservoir in the Qinshui Basin was selected for a case study. The density and strike of the reservoir fractures were obtained and analyzed, and the reservoir was divided into three classes according to fracture density. Finally, in comparison with the gas production data, we confirmed that the inversion results and reservoir classification were both reasonable and credible.     

Solar-, monsoon- and Kuroshio-influenced thermocline depth and sea surface salinity in the southern Okinawa Trough during the past 17,300 years

Factors influencing millennial-scale variability in the thermocline depth (vertical mixing) and sea surface salinity (SSS) of the southern Okinawa Trough (OT) during the past 17,300 years were investigated based on foraminifer oxygen isotope records of the surface dweller Globigerinoides ruber sensu stricto and the thermocline dweller Pulleniatina obliquiloculata in the AMS ¹⁴C dated OKT-3 core. The thermocline depth is influenced by surface thermal buoyancy (heat) flux, in turn controlled by the annual mean insolation at 30°N and the strength of the East Asian winter monsoon (EAWM). Strong insolation and weak EAWM tend to increase buoyancy gain (decrease buoyancy loss), corresponding to shallow thermocline depths, and vice versa. Regional SSS is influenced by the global ice volume, the Kuroshio Current (KC), and vertical mixing. A deep thermocline coincides with a high SSS because strong vertical mixing brings more, saltier subsurface KC water to the surface, and vice versa. Local SSS (excluding the global ice volume effect) became lower in the northern OT than in the southern OT after ~9.2 ka, implying that Changjiang diluted water had stronger influence in the northern sector. SSS show no major changes during the Bølling/Allerød and Younger Dryas events, probably because the KC disturbed the North Atlantic signals. This argues against earlier interpretations of sea surface temperature records of this core. Wavelet and spectral analyses of the Δδ¹⁸O_P-G (δ¹⁸O of P. obliquiloculata minus G. ruber s.s.) and δ¹⁸O_local records display 1,540-, 1,480-, 1,050-, 860-, 640-, and 630-year periods. These are consistent with published evidence of a pervasive periodicity of 1,500 years in global climate as well as EAWM and KC signatures, and a fundamental solar periodicity of 1,000 years and intermediary derived periodicity of 700 years.     

BVOC emissions from mechanical wounding of leaves and branches of <Emphasis Type="Italic">Eucalyptus sideroxylon</Emphasis> (red ironbark)

Emissions of biogenic volatile organic compounds (BVOCs) from mechanical wounding of leaves and branches of plants can contribute to the atmospheric burden of volatile organic compounds (VOCs) in both (a) urban airsheds (from urban garden maintenance) and (b) the global atmosphere (from large scale forest harvesting). These emissions of BVOCs are poorly understood and quantified, and their role in urban and global emissions inventories neglected. This paper presents measurements of the magnitude, duration and composition of emissions of BVOCs, carbon dioxide (CO₂) and methane (CH₄) from freshly cut leaf mulch and wood chips derived from a common eucalypt tree, Eucalyptus sideroxylon (red ironbark), found in southeastern Australian forests and gardens. The emissions of BVOCs from freshly cut and shredded leaves and wood of E. sideroxylon were found to be 2.3 ± 0.6 and 0.05 ± 0.04 mg g^-1 DM (Dry Mass) from leaf mulch and wood chips respectively and to last typically for 1 day following cutting. Three sampling techniques were used for VOC speciation and the 12 most abundant BVOCs released from the mulch materials were identified. The specific BVOCs emitted in order of decreasing abundance from leaf mulch are: (a) stored plant oils, 1,8-cineole, α–pinene and o-cymene which make up the major part of the emissions, (b) a minor contribution from chemicals associated with environmental stress and wound defence, (Z)–3–hexenyl acetate, (E)-2-hexenal and (Z)-3-hexen-1-ol, and (c) a second minor contribution from metabolic products, acetaldehyde and acetone. The observed integrated emissions of BVOCs from leaves following mulching are equivalent to more than half and perhaps all of the likely stored plant oils in the leaves. For the two comparable studies available, one of a plant with stored oils (this study) and one of a plant without stored plant oils, the emissions of leaf wound defence BVOCs are in the same range for both plants. In the plant with stored plant oils, the plant oil emissions are about a factor of 11 larger in emission rate than the plant wound defence BVOCs. A compilation of available leaf wounding BVOC emission studies indicates that for plants with stored plant oils, plant oil emissions dominate, whereas with other plants, leaf wound defence BVOCs dominate the emissions.     

流域地形自动分割研究

本文提供了基于栅格数字高程模型（ＤＥＭ）提取流域空间结构,进行流域自动分割的技术方法。该方法包括地貌形态结构定义、流域结构定义和算法设计３个部分。基于本文设计的算法,可以从ＤＥＭ提取完整的汇流网络和分水线网络,两者组成的多边形集把流域分割成具有拓扑结构的多个子汇流区,这个子汇流区集作为存贮流域信息的基本地貌单元,是流域信息进入地学数据库的框架,也是进行空间分布式机理──过程流域分析模型的基础。并以山西离石黄土丘陵沟壑区的王家沟流域为试验区,获得了较为满意的结果。     

Three-dimensional forward modeling of DC resistivity using the aggregation-based algebraic multigrid method

To speed up three-dimensional (3D) DC resistivity modeling, we present a new multigrid method, the aggregation-based algebraic multigrid method (AGMG). We first discretize the differential equation of the secondary potential field with mixed boundary conditions by using a seven-point finite-difference method to obtain a large sparse system of linear equations. Then, we introduce the theory behind the pairwise aggregation algorithms for AGMG and use the conjugate-gradient method with the V-cycle AGMG preconditioner (AGMG-CG) to solve the linear equations. We use typical geoelectrical models to test the proposed AGMG-CG method and compare the results with analytical solutions and the 3DDCXH algorithm for 3D DC modeling (3DDCXH). In addition, we apply the AGMG-CG method to different grid sizes and geoelectrical models and compare it to different iterative methods, such as ILU-BICGSTAB, ILU-GCR, and SSOR-CG. The AGMG-CG method yields nearly linearly decreasing errors, whereas the number of iterations increases slowly with increasing grid size. The AGMG-CG method is precise and converges fast, and thus can improve the computational efficiency in forward modeling of three-dimensional DC resistivity.