关于矿产资源勘查开采过程中关键节点的探讨

随着采矿等矿业活动在全球范围内的大面积进行,人们对矿山及其周边环境问题的关注度持续增加。电感耦合多接收等离子体质谱仪(MC-ICP-MS)的出现推动了同位素的地球化学研究,也使得同位素示踪技术被广泛应用于探究矿山环境中的各类问题。为强调同位素示踪技术在复杂矿山环境中应用的重要性及其能解决科学问题的多样性,本文调研和分析了截至2022年7月国内外学者公开发表的借助同位素示踪技术测试、分析矿山水文环境中的地球化学过程及污染物来源/影响等方面的论文及其数据,研究区涵盖二十多个国家、四十多个地区。通过总结发现：水体氢、氧同位素示踪技术是矿山水源解析、水力联系研究及酸性矿山废水(AMD)源识别的有效工具;硫酸盐硫、氧同位素示踪技术为研究矿山环境中的硫酸盐来源、AMD酸化过程及污染、细菌硫酸盐的还原作用与元素迁移转化等提供重要支持;重金属(铅、镉、锌、汞等)同位素示踪技术是探究矿山及附近环境中的金属污染来源及不同来源贡献率的有效手段。大量研究表明,虽然同位素技术在解析矿山环境污染物来源和特征污染物迁移转化机制以及揭示矿山水文地球化学过程等方面起到重要作用,但目前的大部分研究局限于应用单一/少数同位素对矿山环境介质进行短时间示踪研究。因此,未来需进一步发展多同位素示踪技术,并对矿山环境中存在的各类问题进行长期、持续地监测调查,提出有效的污染防治新方法。

     

Chemical signature of intermediate water masses along western Portuguese margin

Although the circulation of intermediate water masses in the eastern North Atlantic remains poorly defined, the presence of fresher intermediate waters, the Sub-Artic (SAIW) and the Antarctic Intermediate Water (AAIW), as well the saline intermediate Mediterranean Water (MW), has been tracked using biogeochemical properties. Here we assess the hydrographic and chemical structures of intermediate waters along the western Portuguese margin by examining the vertical distributions and property-property plots of chemical tracers (oxygen and nutrients). AAIW was traced by low oxygen and high nutrients, while SAIW was recognized by low nutrients. The Mediterranean Water (MW) undercurrent is shown to spread towards the eastern flank of Gorringe bank. Concurrently, the fresher waters gained salt by direct incorporation of MW, while this water was enriched in nutrients on its way northward and westward owing, to a great extent, to the entrainment of an AAIW branch. The distributions of nutrients and apparent oxygen utilization are discussed in terms of regional ocean circulation. Our analysis suggests a circulation pattern of the various intermediate waters along the western Portuguese margin: MW extends all over the area, but its presence is more pronounced around cape St. Vincent; SAIW apparently moves southward, reaching the Gorringe bank region, and AAIW flows northward along the coast and around the bank.     

Models of travel demand with endogenous preference change and heterogeneous agents

In the literatures of regional science, urban economics, and urban development planning, a working assumption is that individuals respond to incentives and regulations, given their preferences. Models for planning and policy analyses are used to consider what might occur if the incentives or regulations were different. In these models, however, preferences are usually assumed to be given and stable, and agents are usually assumed to be homogeneous. This paper focuses on the implications of making preferences in models of policy implementation endogenously determined and time varying heterogeneous agents. We consider first the recent literature on intertemporal choice and preference change, which cuts across many disciplines, and more briefly the literature on norm-regarding behavior. We then elaborate a simple model of transportation demand—from a static to a dynamic orientation, from fixed and exogenously given preferences of strictly self-regarding agents to endogenously determined and policy-induced preferences of heterogeneous agents—and illustrate its characteristics with simple numerical examples.     

基于便携式X射线荧光光谱法的深海沉积物现场成分快速检测及适用性评估

远洋粘土、硅藻软泥、铁锰结核及结壳是常见的深海沉积成因物质类型,在开展深海地质调查工作时及时获得这些物质的成分信息,对于提升海上地质调查效率和认知能力十分重要。通过便携式X射线荧光光谱法（pXRF）对菲律宾海深水区域60个样品进行分析,并结合实验室测试结果,评估pXRF能测出的24种元素含量数据的稳定性、准确性及相关性,探讨pXRF在海上地质调查工作中对于深海沉积物现场成分快速检测的适用性。通过综合对比研究发现,Ca、Cu、Fe、K、P、Pb、Sr、Zn、Zr九种元素稳定性、准确性和相关性较好,可以直接用于定性乃至定量研究工作。Al、Ba、Mn、Mo、Ni、Rb、Si、Ti、Th、V十种元素的三项指标参数等级略低,可用于含量高低判定和趋势分析等定性研究工作。Bi、Cs、Mg、Sb、Sc五种元素的测试效果较差,本法不建议使用。取样量充足且分布均匀、湿样烘干至恒量后研磨过筛、封装时充分压实平整、测试时间增长、选择合适的标准物质校正和检验、重点和异常样品增加测试次数等措施的执行有利于获得高精度的pXRF海上测试数据。未来对更多类型和数量深海样品开展pXRF现场分析,将有利于建立更精准的测试方法流程,达到现场快速揭示样品成分特征、确定底质类型及圈定矿化异常等目的,为海上关键决策的部署和实施提供参考和借鉴。     

Application of saddlepoint approximation in reliability analysis of dynamic systems

The application of the saddlepoint approximation to reliability analysis of dynamic systems is investigated. The failure event in reliability problems is formulated as the exceedance of a single performance variable over a prescribed threshold level. The saddlepoint approximation technique provides a choice to estimate the cumulative distribution function (CDF) of the performance variable. The failure probability is obtained as the value of the complement CDF at a specif ied threshold. The method requires computing the saddlepoint from a simple algebraic equation that depends on the cumulant generating function (CGF) of the performance variable. A method for calculating the saddlepoint using random samples of the performance variable is presented. The applicable region of the saddlepoint approximation is discussed in detail. A 10-story shear building model with white noise excitation illustrates the accuracy and effi ciency of the proposed methodology.     

Lateral continuity of basement seismic reflections in 15 Ma ultrafast-spreading crust at ODP Site 1256

The Ocean Drilling Program (ODP) initiated drilling at Site 1256D in the Guatemala Basin, about 1,000 km off the East Pacific Rise to penetrate plutonic rocks, anticipated to be relatively shallow in this region, formed at an ultra-fast spreading rate. IODP Expedition E312 successfully drilled into gabbros at ~1,150 m in basement. Multi-channel seismic traces show weak laterally coherent sub-basement reflections at borehole depths. Synthetic reflectivity seismograms were computed using a Ricker wavelet and impedance profiles from borehole sonic logs. These seismograms show significant sub-basement amplitude peaks. A zero-offset vertical seismic profile, shot on E312, was processed to investigate the authenticity of these reflections and their relationship to borehole geology. A dual scheme of the median filtering and F–K dip filtering was used. Tests with synthetic seismograms indicate the approach is effective at reasonable SNR levels. Downgoing energy is clearly identified but negligible upgoing energy is visible over random noise. These results indicate that lava flows and igneous contacts in upper ocean crust have significant topography on lateral scales less than the Fresnel Zone (~300 m) due to igneous and tectonic processes.     

An Analytical Model For Threshold Crossing Rates Of Concentration Fluctuations In Dispersing Plumes

A simple analytical model is developed for the meanupcrossing rate of plume concentration fluctuations assuming that thisprocess can be well approximated by a lognormal process. The resultingexpression requires only the specification of the in-plume fluctuationintensity and in-plume Taylor micro-time scale and, hence, does notexplicitly involve the joint probability density function of theconcentration and its derivative. The analytical model provides agood fit to some field measurements of the mean upcrossing rate ina dispersing plume.     

Study on uncertainty of inversion and resolution limit in the presence of noise

This study analyzed the uncertainty of inversion and the resolution limit in the presence of noise by means of statistical experiments. The exhaustive method is adopted to obtain the global optimal solution in each experiment. We found that even with small level of noise, solutions fluctuate in a large range for the thin bed. The distribution of solutions in the presence of noise is closely related to the spread of the cost function in the absence of noise. As a result, the area of a certain neighborhood around the true solution on the spread of the cost function in the absence of noise is used to evaluate the uncertainty of inversion and the resolution limit in the presence of noise. In the case that the SNR (signal-to-noise ratio) is 5 in this study, solutions focus around the true solution with a very small uncertainty only when the bed thickness is greater than the reciprocal of the double predominant frequency of the convoluting wavelet.     

An inter-comparison of tidal solutions computed with a range of unstructured grid models of the Irish and Celtic Sea Regions

Three finite element codes, namely TELEMAC, ADCIRC and QUODDY, are used to compute the spatial distributions of the M₂, M₄ and M₆ components of the tide in the sea region off the west coast of Britain. This region is chosen because there is an accurate topographic dataset in the area and detailed open boundary M₂ tidal forcing for driving the model. In addition, accurate solutions (based upon comparisons with extensive observations) using uniform grid finite difference models forced with these open boundary data exist for comparison purposes. By using boundary forcing, bottom topography and bottom drag coefficients identical to those used in an earlier finite difference model, there is no danger of comparing finite element solutions for “untuned unoptimised solutions” with those from a “tuned optimised solution”. In addition, by placing the open boundary in all finite element calculations at the same location as that used in a previous finite difference model and using the same M₂ tidal boundary forcing and water depths, a like with like comparison of solutions derived with the various finite element models was possible. In addition, this open boundary was well removed from the shallow water region, namely the eastern Irish Sea where the higher harmonics were generated. Since these are not included in the open boundary, forcing their generation was determined by physical processes within the models. Consequently, an inter-comparison of these higher harmonics generated by the various finite element codes gives some indication of the degree of variability in the solution particularly in coastal regions from one finite element model to another. Initial calculations using high-resolution near-shore topography in the eastern Irish Sea and including “wetting and drying” showed that M₂ tidal amplitudes and phases in the region computed with TELEMAC were in good agreement with observations. The ADCIRC code gave amplitudes about 30 cm lower and phases about 8° higher. For the M₄ tide, in the eastern Irish Sea amplitudes computed with TELEMAC were about 4 cm higher than ADCIRC on average, with phase differences of order 5°. For the M₆ component, amplitudes and phases showed significant small-scale variability in the eastern Irish Sea, and no clear bias between the models could be found. Although setting a minimum water depth of 5 m in the near-shore region, hence removing wetting and drying, reduced the small-scale variability in the models, the differences in M₂ and M₄ tide between models remained. For M₆, a significant reduction in variability occurred in the eastern Irish Sea when a minimum 5-m water depth was specified. In this case, TELEMAC gave amplitudes that were 1 cm higher and phases 30° lower than ADCIRC on average. For QUODDY in the eastern Irish Sea, average M₂ tidal amplitudes were about 10 cm higher and phase 8° higher than those computed with TELEMAC. For M₄, amplitudes were approximately 2 cm higher with phases of order 15° higher in the northern part of the region and 15° lower in the southern part. For M₆ in the north of the region, amplitudes were 2 cm higher and about 2 cm lower in the south. Very rapid M₆ tidal-phase changes occurred in the near-shore regions. The lessons learned from this model inter-comparison study are summarised in the final section of the paper. In addition, the problems of performing a detailed model–model inter-comparison are discussed, as are the enormous difficulties of conducting a true model skill assessment that would require detailed measurements of tidal boundary forcing, near-shore topography and precise knowledge of bed types and bed forms. Such data are at present not available.     

Tectonics of southwestern Mexico,isotopic evidence,nuclear Central America,Late Cretaceous break up

In the paleogeographic reconstruction of Mexico and northern Central America, an ever-increasing amount of evidence shows that the entire region is a collage of suspect terranes transported from abroad, whose timing and sense of motion are now beginning to be understood. Among these, the Chortis block (nuclear Central America) and the Baja California Peninsula have been proposed as pieces of continent separated from the Pacific coast of southwestern Mexico, that have moved either southeastward by the Farallon plate or northwestward by the Kula plate. Previous studies mainly confined to the northern margin of the Chortis block, confirmed a left-lateral displacement of 130 km in Neogene time. Further studies made northwestward along the Mexican coast provided a better understanding of magmatic and metamorphic processes in the area, and suggested times of detachment increased to 30 Ma, 40 Ma, and 66 Ma. The pre-detachment westernmost position of the block has changed, depending on the model chosen, from Puerto Vallarta and beyond, to the current position. Here we show that the isotopic mineral ages from coastal granites along the coast from Puerto Vallarta, Jalisco (80 Ma) to Puerto Angel, Oaxaca (11 Ma) record systematic decrease of cooling ages from NW to SE. This pattern is interpreted to result from the progressive uplift of rocks exposed at the present-day coast in that direction, such uplift occurred in response to the development of the Middle America Trench at the newly formed continental margin when the Chortis block was sliding at an average rate of 1.5 cm/year in a sinistral sense to its present position. Our results also constrain the position of the Kula-Farallon spreading axis north of Puerto Vallarta. These observations led us to conclude that several indicators point to this time and region for the onset of strike-slip drifting of the Chortis block toward its current position. Here, we also present several view points in terms of other possilble interpretations to different tectonic, geologic and isotopic data sets published recently by different authors.