基于数值模拟的再生水对地下水质影响实例

为了探讨在再生水地表回灌条件下,再生水对地下水中氨素和主要盐分浓度变化的影响,选取现实存在的清河再生水湿地公园为例,以氨氮、Cl-离子为主要研究对象,利用Hydrus-1D溶质运移模型和Modflow、MT3DMSD耦合地下水模型,模拟了再生水长期利用后的地下水氨氮和Cl-指标变化情况。研究结果显示,再生水到达包气带底部时,水中氮氨、Cl-离子含量高于研究区潜水地下水中的相应指标含量,氨氮超出地下水质量Ⅲ类标准。预测5~20年后,潜水中氨氮超标面积不断扩展,20年后将达14km2,氨氮超标垂向延伸至第一层弱透水层;Cl-离子在潜水中不超标,但在湿地公园周边浓度有所增加。湿地公园再生水的长期地表回灌将给当地地下水质量带来负面影响。     

Crust structure of the northeastern Barents Sea basin area

Processing of data from regional geophysical surveys completed in the northern Barents Sea has provided updates to gravity and magnetic databases, structural maps of seismic interfaces, and positions of anomaly sources, which made a basis for 3D density and magnetic models of the crust. The new geological and geophysical results placed constraints on the boundaries between basement blocks formed in different settings and on the contours of deposition zones of different ages in the northeastern Barents Sea. The estimated thicknesses of sedimentary sequences that formed within certain time spans record the deposition history of the region. There is a 20-50 km wide deep suture between two basins of Mesozoic and Paleozoic ages in the eastern part of the region, where pre-Late Triassic reflectors have no clear correlation. The suture slopes eastward at a low angle and corresponds to a paleothrust according to seismic and modeling data. In the basement model, the suture is approximated by a zone of low magnetization and density, which is common to active fault systems. The discovery of the suture has important geological and exploration implications.     

Role of heterogeneities in jointing (fracturing) of geological media: Numerical analysis of fracture mechanisms

Fracture process is investigated using finite-difference simulations with a new constitutive model. It is shown that both geometry and fracture mechanism itself depend on the preexisting heterogeneities that are stress concentrators. In the brittle regime (low pressure, P), Mode-I fractures propagate normal to the least stress σ₃ from the imposed weak zones. At high P, shear deformation bands are formed oblique to σ₃. At intermediate values of P, the fracture process involves both shear banding and tensile cracking and results in the initiation and propagation of pure dilation bands. The propagating band tip undulates, reacting on the failure mechanism changes, but its global orientation is normal to σ₃. The σ₃-normal fractures are joints. There are thus two types of joints resulting from Mode-I cracking and dilation banding, respectively. The obtained numerical results are in good agreement with and explain the results from previous similar experimental study.     

内蒙古巴林右旗晚二叠世埃达克质火山岩特征及其地质意义

内蒙古东南部巴林右旗地区发育晚二叠世埃达克质火山岩,岩石组合为安山岩、粗安岩、英安岩及辉石安山岩,其LA-ICP-MS锆石U-Pb测年结果为256.7±2.7Ma,指示其形成于晚二叠世。地球化学特征显示,该套火山岩属准铝质-弱铝质中钾钙碱性岩石系列,具富Si(SiO_256%)、高Al(Al_2O_315%)、富Na、贫K、高Sr、低Yb和Y等特征,Na_2O/K_2O值为2.33~3.90,Mg~#值为35.3~60.8;稀土元素总量为96.69×10~(-6)~192.4×10~(-6)、轻重稀土元素分馏较明显((La/Yb)_N值为6.27~13.82),具正的Eu(δEu=1~1.67)异常,在原始地幔标准化蛛网图中,富集大离子亲石元素Rb、Ba、U,亏损高场强元素Nb、Ta,为O型(大洋型)埃达克质火山岩地球化学特征。综合区域资料,巴林右旗埃达克质火山岩是残留在地幔中的古亚洲洋残余洋壳部分熔融并受到地幔橄榄岩混染形成的,暗示晚二叠世存在古亚洲洋向华北板块俯冲消亡事件。     

Calibration of cyclic constitutive models by oscillating functions

Most structures are subjected to more cyclic loads during their life time than static loads. These cyclic action could be a result of either natural or man-made activities and may lead to soil failure. In order to understand the response of the foundation and its interaction with these complex cyclic loadings, various researchers have over the years developed different constitutive models. Although a lot of research is being carried out on these relatively new models, little or no details exist in literature about the model-based identification of the cyclic constitutive parameters which to a large extent govern the quality of the model output. This could be attributed to the difficulties and complexities of the inverse modeling of such complex phenomena. A variety of optimisation strategies are available for the solution of the sum of least-squares problems as usually done in the field of model calibration. However, for the back analysis (calibration) of the soil response to oscillatory load functions, this article gives insight into the model calibration challenges and also puts forward a method for the inverse modeling of cyclic loaded foundation response such that high-quality solutions are obtained with minimum computational effort.     

Middle Jurassic–Early Cretaceous tectonic evolution of the Bayanhushuo area,southern Great Xing’an Range,NE China: constraints from zircon U–Pb geochronological and geochemical data of volcanic and subvolcanic rocks

This study presents new zircon U–Pb geochronology, geochemistry, and zircon Hf isotopic data of volcanic and subvolcanic rocks that crop out in the Bayanhushuo area of the southern Great Xing’an Range (GXR) of NE China. These data provide insights into the tectonic evolution of this area during the late Mesozoic and constrain the evolution of the Mongol–Okhotsk Ocean. Combining these new ages with previously published data suggests that the late Mesozoic volcanism occurred in two distinct episodes: Early–Middle Jurassic (176–173 Ma) and Late Jurassic–Early Cretaceous (151–138 Ma). The Early–Middle Jurassic dacite porphyry belongs to high-K calc-alkaline series, showing the features of I-type igneous rock. This unit has zircon ε_Hf(t) values from +4.06 to +11.62 that yield two-stage model ages (T_DM2) from 959 to 481 Ma. The geochemistry of the dacite porphyry is indicative of formation in a volcanic arc tectonic setting, and it is derived from a primary magma generated by the partial melting of juvenile mafic crustal material. The Late Jurassic–Early Cretaceous volcanic rocks belong to high-K calc-alkaline or shoshonite series and have A₂-type affinities. These volcanics have ε_Hf(t) and T_DM2 values from +5.00 to +8.93 and from 879 to 627 Ma, respectively. The geochemistry of these Late Jurassic–Early Cretaceous volcanic rocks is indicative of formation in a post-collisional extensional environment, and they formed from primary magmas generated by the partial melting of juvenile mafic lower crust. The discovery of late Mesozoic volcanic and subvolcanic rocks within the southern GXR indicates that this region was in volcanic arc and extensional tectonic settings during the Early–Middle Jurassic and the Late Jurassic–Early Cretaceous, respectively. This indicates that the Mongol–Okhotsk oceanic plate was undergoing subduction during the Early–Middle Jurassic, and this ocean adjacent to the GXR may have closed by the Late Middle Jurassic–Early Late Jurassic.     

Early Carboniferous ophiolite in central Qiangtang,northern Tibet: record of an oceanic back-arc system in the Palaeo-Tethys Ocean

Zircon U–Pb dating of two samples of metagabbro from the Riwanchaka ophiolite yielded early Carboniferous ages of 354.4 ± 2.3 Ma and 356.7 ± 1.9 Ma. Their positive zircon εHf(t) values (+7.9 to +9.9) indicate that these rocks were derived from a relatively depleted mantle. The metagabbros can be considered as two types: R1 and R2. Both types are tholeiitic, with depletion of high-field-strength elements (HFSE) and enrichment of large-ion lithophile elements (LILE) similar to those of typical back-arc basin basalts (BABB), such as Mariana BABB and East Scotia Ridge BABB. Geochemical and isotopic characteristics indicate that the R1 metagabbro originated from a back-arc basin spreading ridge with addition of slab-derived fluids, whereas the R2 metagabbro was derived from a back-arc basin mantle source, with involvement of melts and fluids from subducted ocean crust. The Riwanchaka ophiolite exhibits both mid-ocean ridge basalts- and arc-like geochemical affinities, consistent with coeval ophiolites from central Qiangtang. Observations indicate that the Qiangtang ophiolites developed during the Late Devonian–early Carboniferous (D₃–C₁) in a back-arc spreading ridge above an intra-oceanic subduction zone. Based on our data and previous studies, we propose that an oceanic back-arc basin system existed in the Longmuco–Shuanghu–Lancang Palaeo-Tethys Ocean during the D₃–C₁ period.     

准噶尔盆地腹部地区生烃增压定量化评价——东道海子北凹陷为例

目前准噶尔盆地腹部地区超压演化过程仅提供了相对时间、概念、定性的认识,并未达到定量化。本文以东道海子北凹陷为例,基于流体包裹体、单井一维和剖面二维烃源岩成熟生烃史模拟技术,采用生烃增压模型定量化评价了侏罗纪八道湾组顶部和底部烃源岩的生烃增压演化过程。研究结果表明侏罗纪八道湾组为东道海子北凹陷主力烃源岩,具有幕式排烃的特点。八道湾组底部烃源岩经历了3 次超压增加和两次超压释放的过程,顶部烃源岩经历了两次超压增加和一次超压释放的过程。八道湾底部烃源岩3 次超压增加的时间分别为距今170～70 Ma、55～23 Ma、10～0 Ma,对应最大压力系数均为2.0,最大剩余流体压力50 MPa,两次超压释放的时间分别为70～55 Ma、23～10 Ma;八道湾组顶部烃源岩两次超压增加的时间分别为130～40 Ma、23～0 Ma,对应流体最大压力系数最大亦可达2.0 左右,最大剩余流体压力44 MPa,超压释放的时间为40～23 Ma。     

基于人工鱼群算法和SCA等效理论的井中横波速度预测方法

在建立页岩岩石物理模型的基础上,根据等效自相容近似(SCA)岩石物理模型,构建出岩石的纵波速度、横波速度与岩石密度、组分和孔隙度等的定量关系,得出使理论纵波速度和实际纵波速度最接近的孔隙纵横比,进而将该孔隙纵横比作为约束条件来实现横波速度预测。反演算法利用人工鱼群算法来计算最佳孔隙纵横比,并将预测的横波速度与实际测得的横波速度对比,证明了人工鱼群算法的有效性。     

波浪在珊瑚岸礁礁坪上传播变形的数值研究

波浪的传播变形对珊瑚岛礁附近的营养物质输送以及近岸珊瑚沙运动都起着关键的作用,同时也是礁盘顶部建筑物在设计和施工过程中需要重点考虑的因素,因此,对波浪在礁坪上的传播变形进行研究具有重要意义。利用FUNWAVE 2.0数值模型模拟了不规则波浪在不同形状岸礁上的传播变形过程,讨论了礁前斜坡坡度以及珊瑚底面粗糙程度对波浪在礁坪上传播变形的影响。结果显示,随着礁坡斜率增大,不对称度参数绝对值的增长趋势逐渐变缓,波浪在光滑礁坪上的不对称度绝对值要略大于粗糙底面上的值,而不同地形下计算得到的波浪偏度参数则基本相同。说明礁坡坡度对波浪的不对称度特征具有显著影响,礁坪粗糙度的影响则相对较小,而上述因素对波浪偏度的影响则完全可以忽略。给出了考虑这两个影响因素的礁坪区域内波浪偏度和不对称度经验公式。