有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

Numerical Simulation of Residual Circulation due to Bottom Rouglmess Variability Under Tidal Flows in A Semi-Enclosed Bay

Nowadays there are some chronic serious environmental problems, such as eutrophication, blue tide and so on, in a complicated coastal zone or a semi-enclosed bay, because the water exchanges between an inner bay and an outer sea is weak compared with the supply of contaminant. Under this situation, a method to improve the water quality by 3-dimensional small unsymmetrical structures has been proposed by Komatsu et al. In this paper, several numerical simulations of the tidal current and concentration for various arrangements of bottom roughness in a semi-enclosed model bay are carfled out with a depth-averaged 2-D numerical model. The model is solved by the hybrid finite analytic method with nonstaggered grid. And the SIMPLES algorithm with Rhie and Chow＇ s momentum interpolation technique is used for the simulation. The effect of Komatsu＇ s method for water purification is examined by numerical simulation. The result of numerical experiment indicates that it is possible to generate a new tidal residual current and to activate a tidal exchange by bottom roughness arrangement only.     

Numerical Analysis of Wave Forces on Arbitrarily Shaped Multi-Pillars Over Varying Topography

The modified hybrid element method (MHEM) is utilized to predict and analyze wave forces on arbitrarily shaped multiple bodies. This method can be applied to waves of all water depths, i. e. shallow, intermediate, and deep waters, on slowly varying seabed. The MHEM employs the ICCG method to save CPU and storage, thus the computation of wave forces for large multi-body systems can be carried out on microcomputers. Numerical results of the present method are compared with experimental data and other solutions. It is shown that the MHEM provides more accurate solutions of the wave forces than other numerical methods do. Therefore, the methodology presented herein can be used in the design of coastal and ocean structures.     

Influence of large-diameter pipe pile driving on surrounding marine soils and adjacent submarine pipelines

Abstract

With the large-scale development and utilization of ocean resources and space, it is inevitable to encounter existing submarine facilities in pile driving areas, which necessitates a safety assessment. In this article, by referring to a wharf renovation project as a reference, the surrounding soil response and buried pipe deformation during pile driving in a near-shore submarine environment are investigated by three-dimensional (3D) numerical models that consider the pore water effect. Numerical studies are carried out in two different series: one is a case of a single pile focusing on the effect of the minimum plane distance of the pile–pipe, and the other is a case of double piles focusing on the effect of the pile spacing.     

利用臭氧氧化实现复合生物反应器污泥减量

采用复合生物反应器,对应用臭氧氧化实现污泥减量进行了研究。2个相同的复合生物反应器平行运行,1个作为对照系统,1个作为氧化系统。反应器内装有半软性填料,投加量为10%。曝气池中悬浮污泥浓度为1 500 mg/L左右,生物膜浓度为2 000 mg/L左右。试验结果表明,随着臭氧氧化污泥比例的增加,污泥表观产率系数也随之降低,当臭氧投量为0.05 g O3/gSS,每天氧化的污泥分别为反应器内污泥的10%,20%,30%时,污泥表观产率系数分别减少28.2%,44.9%,75.8%。虽然随着污泥氧化比例的增大,氧化系统出水CODcr略有增加,但氧化系统仍能保持其生物处理能力,CODcr去除率在92%以上;2个系统之间氨氮的去除率相差不大,氧化系统的硝化能力基本没有受到臭氧氧化的影响。     

城区水文循环机制与改善策略分析

在城区化过程中,由于人口快速成长且集中,随之而来的建筑物、道路等不透水面密度增加,导致地下水补给减少及破坏水文循环。为了解土地利用方式改变对长期水文量影响,建立一个适宜的城区水循环模式是非常重要的。选定台湾台北的六馆抽水站集水区进行分析,模式参数经过敏感度分析、率定与验证。采用改变不同透水铺面设置比例的方案,评估其对径流量、入渗量、蒸发散量的改变。并分析设置不同比例屋顶雨水贮集系统对地表径流量及供水量的改变。所建立的城区水循环模式是评估土地利用改变造成水文量的改变的一个有效工具,并可提供决策者做为城区土地利用决策的参考。     

临涣矿区界沟煤矿工程地质条件研究

根据主采煤层顶底板岩体岩性类型及结构面发育程度对其稳定性进行了评价：72煤层顶板以中等稳定一不稳定类型为主,82煤层顶板以中等稳定类型为主,10煤层顶板以中等稳定一稳定类型为主。在大量统计见煤点硬质岩百分含量（k值）的基础上,结合顶板岩石的单轴抗压强度,对各主采煤层顶底板工程地质类型进行了划分：72煤层以一、二类顶板为主,82、10煤层均以五类顶板为主。综合评价该矿工程地质条件为中等即Ⅲ类二型。根据研究成果,指出该矿在采掘过程中可能出现的工程地质问题,并依此提出了治理建议。     

三维地震技术在煤层顶板稳定性研究中的应用

龙固井田为全隐蔽的华北型煤田,位于巨野煤田中部,其首采扩大区主采煤层为3号煤层。考虑3煤层顶板稳定性主要受其顶板的构造信息和岩性信息影响,因此首先依据三维地震勘探综合解释成果及波阻抗反演解释成果对二者进行定量化,然后对波阻抗数据进行归一化处理,使得波阻抗数据和量化后的构造数据具有相同的变化范围及等量贡献。在此距离范围内构造和岩性的权值各为0．5,依此生成综合因素煤层顶板稳定性隶属度。分析3煤层顶板以上10m、20m处的综合因素煤层顶板稳定性隶属度图可以发现,该区3煤层顶板稳定性比较好。且其稳定性主要是受构造因素控制,岩性因素相对影响较小。     

Hybrid seismic modeling based on discrete-wave number and finite-difference methods

Any calculation of seismic wave propagation comprising the seismic source, the travel path, and the receiver site in a single finite-difference (FD) model requires a considerable amount of computer time and memory. Moreover, the methods currently available for including point sources in the 2D FD calculations are far-field approximations only. Therefore we have developed a new hybrid method for treating the seismic wave fields at localized 2D near-surface structures embedded in a 1D background medium, and excited by a point source. The source radiation and propagation in the background model is solved by the discrete-wave number (DW) method, while the propagation in the local 2D structure is calculated by the FD method. The coupling between the two sets of calculations is performed on a rectangular excitation box surrounding the local structure. We show the usefulness of the method in ground-motion studies where both near-field source effects and local site effects are important. Technical problems connected with the inconsistency between the 3D source radiation and the 2D FD calculation are minor for the relatively distant in-plane point explosive sources, but are more serious for the in-plane dislocation sources.