二维潮流数学模型的四边形等参单元法

本文在有限差分和有限元法的基础上，采用四边形等参单元，对在河海域进行了潮流数学模型计算，计算表明：该方法既能较好地拟合边界，又可减少网格划分等前期工作，所得结果与实测资料基本一致。     

文县天池环境调查

根据文县天池首次实地科学考察资料，本文从自然地理，物理化学、生物及旅游诸方面论述了天池的基本特。     

The combined effect of F and H2O on interdiffusion between peralkaline dacitic and rhyolitic melts

The effective binary diffusion coefficient (EBDC) of silicon has been measured during the interdiffusion of peralkaline, fluorine-bearing (1.3 wt% F), hydrous (3.3 and 6 wt% H₂O), dacitic and rhyolitic melts at 1.0 GPa and temperatures between 1100°C and 1400°C. From Boltzmann-Matano analysis of diffusion profiles the diffusivity of silicon at 68 wt% SiO₂ can be described by the following Arrhenius equations (with standard errors): $$\begin{gathered} {\text{with 1}}{\text{.3 wt\% F and 3}}{\text{.3\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.66}} \times {\text{10}}^{ - {\text{9}}} } \\ { - {\text{1}}{\text{.86}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{86}}{\text{.1}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ {\text{with 1}}{\text{.3 wt\% F and 6}}{\text{.0\% H}}_{\text{2}} {\text{O:}} \hfill \\ {\text{D}}_{{\text{Si}}} = \begin{array}{*{20}c} { + {\text{3}}{\text{.59}}} \\ {{\text{3}}{\text{.51}} \times {\text{10}}^{ - {\text{8}}} } \\ { - {\text{1}}{\text{.77}}} \\ \end{array} {\text{exp}}\left( {{{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} \mathord{\left/ {\vphantom {{ - {\text{109}}{\text{.5}} \pm {\text{8}}{\text{.9}}} {{\text{RT}}}}} \right. \kern-\nulldelimiterspace} {{\text{RT}}}}} \right) \hfill \\ \end{gathered} $$ where D is in m²s^?1 and activation energies are in kJ/mol. Diffusivities measured at 64 and 72 wt% SiO₂ are only slightly different from those at 68 wt% SiO₂ and frequently all measurements are within error of each other. Silicon, aluminum, iron, magnesium, and calcium EBDCs were also calculated from diffusion profiles by error function inversion techniques assuming constant diffusivity. With one exception, silicon EBDCs calculated by error function techniques are within error of Boltzmann-Matano EBDCs. Average diffusivities of Fe, Mg, and Ca were within a factor of 2.5 of silicon diffusivities whereas Al diffusivities were approximately half those of silicon. Alkalies diffused much more rapidly than silicon and non-alkalies, however their diffusivities were not quantitatively determined. Low activation energies for silicon EBDCs result in rapid diffusion at magmatic temperatures. Assuming that water and fluorine exert similar effects on melt viscosity at high temperatures, the viscosity can be calculated and used in the Eyring equation used to determine diffusivities, typically to within a factor of three of those measured in this study. This correlation between viscosity and diffusivity can be inverted to calculate viscosities of fluorine- and water-bearing granitic melts at magmatic temperatures; these viscosities are orders of magnitude below those of hydrous granitic melts and result in more rapid and effective separation of granitic magmas from partially molten source rocks. Comparison of Arrhenius parameters for diffusion measured in this study with Arrhenius parameters determined for diffusion in similar compositions at the same pressure demonstrates simple relationships between Arrhenius parameters, activation energy-E_a, kJ/mol, pre-exponential factor-D_o, m²s^?1, and the volatile, X=F or OH^?, to oxygen, O, ratio of the melt {(X/X+O)}: $$\begin{gathered} {\text{E}}a = - {\text{1533\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }} + {\text{213}}{\text{.3}} \hfill \\ {\text{D}}_{\text{O}} = {\text{2}}{\text{.13}} \times {\text{10}}^{ - {\text{6}}} {\text{exp}}\left[ { - {\text{6}}{\text{.5\{ }}{{\text{X}} \mathord{\left/ {\vphantom {{\text{X}} {\left( {{\text{X}} + {\text{O}}} \right)}}} \right. \kern-\nulldelimiterspace} {\left( {{\text{X}} + {\text{O}}} \right)}}{\text{\} }}} \right] \hfill \\ \end{gathered} $$ These relationships can be used to estimate diffusion in various melts of dacitic to rhyolitic composition containing both fluorine and water. Calculations for the contamination of rhyolitic melts by dacitic enclaves at 800°C and 700°C provide evidence for the virtual inevitability of diffusive contamination in hydrous and fluorine-bearing magmas if they undergo magma mixing of any form.     

Lithofacies modeling by multipoint statistics and economic evaluation by NPV volume for the early Cretaceous Wabiskaw Member in Athabasca oilsands area,Canada

The static modeling and dynamic simulation are essential and critical processes in petroleum exploration and development. In this study, lithofacies models for Wabiskaw Member in Athabasca, Canada are generated by multipoint statistics(MPS) and then compared with the models built by sequential indicator simulation(SIS). Three training images(Tls) are selected from modern depositional environments;the Orinoco River Delta estuary, Cobequid bay-Salmon River estuary, and Danube River delta environment. In order to validate lithofacies models, average and variance of similarity in lithofacies are calculated through random and zonal blind-well tests.In random six-blind-well test, similarity average of MPS models is higher than that of SIS model. The Salmon MPS model closely resembles facies pattern of Wabiskaw Member in subsurface. Zonal blind-well tests show that successful lithofacies modeling for transitional depositional setting requires additional or proper zonation information on horizontal variation, vertical proportion, and secondary data.As Wabiskaw Member is frontier oilsands lease, it is impossible to evaluate the economics from production data or dynamic simulation. In this study, a dynamic steam assisted gravity drainage(SAGD)performance indicator(SPIDER) on the basis of reservoir characteristics is calculated to build 3 D reservoir model for the evaluation of the SAGD feasibility in Wabiskaw Member. SPIDER depends on reservoir properties, economic limit of steam-oil ratio, and bitumen price. Reservoir properties like porosity,permeability, and water saturation are measured from 13 cores and calculated from 201 well-logs. Three dimensional volumes of reservoir properties are constructed mostly based on relationships among properties. Finally, net present value(NPV) volume can be built by equation relating NPV and SPIDER. The economic area exceeding criterion of US$ 10,000 is identified, and the ranges of reservoir properties are estimated. NPV-volume-generation workflow from reservoir parameter to static model provides costand time-effective method to evaluate the oilsands SAGD project.     

Temperature‐Activated Auto‐Decomposition Reactions: An Under‐Utilized In Situ Remediation Solution

The remediation industry has witnessed multiple innovations arising from a greater understanding of the physical, chemical, and biological processes that control the fate and transport of chemicals in the subsurface environment. In addition, increasing emphasis is being placed on remediation solutions that are greener, simpler, and more resource efficient. The positive impacts that can be derived from this emphasis include reduced energy consumption, reduced waste emissions, and lower costs. Temperature‐activated auto‐decomposition reactions represent a potentially underutilized option for the in situ remediation of certain organic contaminants, and an option that can be both highly effective and greener than other available technologies.     

Vesiculation of high fountaining Hawaiian eruptions: episodes 15 and 16 of 1959 Kīlauea Iki

The 1959 summit eruption of Kīlauea volcano produced the highest recorded Hawaiian fountain in Hawai‘i. Quantitative analysis of closely spaced samples from the final two high-fountaining episodes of the eruption result in a fine-scale textural study of pyroclasts and provide a record of postfragmentation processes. As clast vesicularity increases, the vesicle number density decreases and vesicle morphology shifts from small and round to larger and more irregular. The shift in microtexture corresponds to greater degrees of postfragmentation expansion of clasts with higher vesicularity. We suggest the range of clast morphologies in the deposit is related to thermal zonation within a Hawaiian fountain where the highest vesicularity clasts traveled in the center and lowest traveled along the margins. Vesicle number densities are greatest in the highest fountaining episode and therefore scale with intensity of activity. Major element chemical analyses and fasciculate crystal textures indicate microlite-rich zones within individual clasts are portions of recycled lava lake material that were incorporated into newly vesiculating primary melt.     

远震体波有限破裂过程反演敏感性测试及其在汶川地震上的应用

本文利用远震P和SH波反演得到2008年5月12日汶川大地震(M_W=7.9)的一系列有限破裂模型。使用一种基于小波变换的模拟退火非线性反演方法, 我们将主断层划分成若干个子断层, 在反演时同时确定每个子断层上的滑移量、 滑动角、 上升时间(rise time)以及平均破裂速度。我们首先根据一个假定的破裂模型生成理论地震图, 将该理论地震数据作为输入进行反演, 对该有限破裂反演方法进行了一系列测试, 以验证反演对断层倾角、 平均破裂速度、 最大破裂深度等参数的敏感性。然后我们采用4个不同倾角的断层面来对汶川地震远震体波记录进行反演。结果表明, 若对只在一个断层面上模拟该地震, 30°倾角是个较为合适的值。反演的结果还表明, 此次地震有两个主要的能量释放区域, 并且主断层面存在倾角变化的可能性。在将来的研究中, 可以结合GPS, InSAR测地学以及强震等数据, 来对强震的破裂过程做更细致的研究。