时域瞬变场电磁场有限差分法

模拟电磁场在地下介质中的传播规律是理解地下介质响应的重要手段，直接应用有限差分法在时间域对二维、三维瞬变场进行数值分析是一种有效方法。应用反映电磁场基本规律的麦克斯韦方程组，导出时域场量的齐次扩散方程，对所研究的空间区域作差分离散，源作为初始条件加入，采用合适的边界条件，利用差分方程进行计算，从而得到场量的数值结果，展现瞬变场在地下随时间扩散的全过程，对电磁场的瞬态响应有更加直观的了解。     

南秦岭十里坪锑矿床成矿时代及成因的初步研究

十里坪锑矿床受赵川陆缘隆_滑构造的主滑脱拆离带的控制。矿体呈脉状赋存于韧_脆性主滑脱带上部的脆性次级断层_节理中,矿石类型主要为萤石石英辉锑矿型。围岩为太古宙_元古宙变质岩系,围岩蚀变弱。成矿流体属H2O_CO2_NaCl体系,流体包裹体盐度w(NaCleq)为3.6%～11.3%,均一温度为109～232℃,形成压力大致为800×105Pa。硫、铅同位素研究表明,矿质主要来源于变火山岩围岩;氢、氧同位素显示,成矿流体以大气降水为主,初步将该矿床定为变质岩源就地式大气降水热液矿床。矿石中萤石Sm_Nd等时线年龄为(392±24)Ma,与南秦岭北部晚古生代拗陷区热水喷流_沉积成矿时代相一致,它们都形成于秦岭微板块泥盆纪非造山裂解阶段。     

Web-cam’s potential for collaborative activities in the Earth sciences

We report here our experiences from using easily acquired web-cam components for collaborative ventures in the Earth sciences. We have used a variety of hardware and different software. We demonstrate from various locations in the U.S.A. the feasibility of using web-cam in multitudinous activities, relevant for conducting research and knowledge dissemination. We summarize the quality of the connections from the various combinations of communicating parties. Today, web-cams can be utilized as an economical and viable means of point-to-point communication in the Earth science community. Greater bandwidth is sorely needed for activities such as multi-party conferencing on the present internet network. We propose that web-cam can be readily deployed as a web-service for facilitating collaborative research over the GRID infrastructure, using the middleware Narada-Brokering. Web-cam will play an important role in the emerging field of geoinformatics.     

Indicator Simulation Accounting for Multiple-Point Statistics

Geostatistical simulation aims at reproducing the variability of the real underlying phenomena. When nonlinear features or large-range connectivity is present, the traditional variogram-based simulation approaches do not provide good reproduction of those features. Connectivity of high and low values is often critical for grades in a mineral deposit. Multiple-point statistics can help to characterize these features. The use of multiple-point statistics in geostatistical simulation was proposed more than 10 years ago, on the basis of the use of training images to extract the statistics. This paper proposes the use of multiple-point statistics extracted from actual data. A method is developed to simulate continuous variables. The indicator kriging probabilities used in sequential indicator simulation are modified by probabilities extracted from multiple-point configurations. The correction is done under the assumption of conditional independence. The practical implementation of the method is illustrated with data from a porphyry copper mine.     

Applicability of a Backprojection Algorithm to Reconstruct Images of Subsurface Horizontal Planes for Laboratory Experiments in Electrical Resistance Tomography

This paper describes the use of a backprojection algorithm to reconstruct subsurface images of the electrical resistivity in horizontal planes parallel to the surface. The algorithm can be applied to detect buried objects such as tanks or pipes and possible leakages from them. Two imaging strategies are compared: juxtaposition of vertical planes, and 3D reconstruction from the sensitivity matrix corresponding to the entire volume whose surface is explored. The electrode arrays used for voltage measurement are the dipole–dipole array and a modified Schlumberger array. A personal computer controls current injection, electrode switching, and voltage detection. The system injects 1 kHz, 20 V peak-to-peak square waveforms, thus avoiding electrode polarization effects. Experimental laboratory measurements show that the algorithm detects localized objects such as an insulating sphere and a conductive cylinder immersed in water. Furthermore, covering half of the cylinder by a rubber sleeve to simulate a nonconductive leak, yields a distinct image for the leak. The backprojection algorithm does not need any regularization parameter and it is very fast in inverting the sensitivity matrix because it approximates the inverse matrix by its transposed. The dipole–dipole array usually yields a lower overall pixel error than the modified Schlumberger array but both allow the detection of simulated underground leaks.     

Signatures in flowing fluid electric conductivity logs

Flowing fluid electric conductivity logging provides a means to determine hydrologic properties of fractures, fracture zones, or other permeable layers intersecting a borehole in saturated rock. The method involves analyzing the time-evolution of fluid electric conductivity (FEC) logs obtained while the well is being pumped and yields information on the location, hydraulic transmissivity, and salinity of permeable layers. The original analysis method was restricted to the case in which flows from the permeable layers or fractures were directed into the borehole (inflow). Recently, the method was adapted to permit treatment of both inflow and outflow, including analysis of natural regional flow in the permeable layer. A numerical model simulates flow and transport in the wellbore during flowing FEC logging, and fracture properties are determined by optimizing the match between simulation results and observed FEC logs. This can be a laborious trial-and-error procedure, especially when both inflow and outflow points are present. Improved analyses methods are needed. One possible tactic would be to develop an automated inverse method, but this paper takes a more elementary approach and focuses on identifying the signatures that various inflow and outflow features create in flowing FEC logs. The physical insight obtained provides a basis for more efficient analysis of these logs, both for the present trial and error approach and for a potential future automated inverse approach. Inflow points produce distinctive signatures in the FEC logs themselves, enabling the determination of location, inflow rate, and ion concentration. Identifying outflow locations and flow rates typically requires a more complicated integral method, which is also presented in this paper.     

沙坡头没蒿和柠条叶片δ13C 和 N 含量的季节变化特征

通过对腾格里沙漠东南缘沙坡头人工固沙区油蒿和柠条叶片稳定碳同位素分辨率（△）和N 含量的测定，研究了不同种植方式下油蒿和柠条叶片△ 和 N 含量的季节变化及其关系。结果表明：在不同种植方式下，两种植物的△ 与 N 含量在不同月份之间均存在显著差异。油蒿叶片 △ 极显著高于柠条的：在不同种植方式下，油蒿单种下 △ 显著高于混种，柠条则反之。 油蒿叶片 N 含量极显著低于柠条的：不同种植方式下柠条叶片 N 含量无显著差异，而单种油蒿叶片 N 含量则显著低于混种。在 △ 与 N 含量的关系中，柠条叶片 △ 与 N 在两种生境、单种和混种下均极显著正相关，雨明柠条叶片 N 含量可以作为其 △ 与 N 含量的季节变化及其关系在不同微生境有所差异。     

Continental rift systems and anorogenic magmatism

Precambrian Laurentia and Mesozoic Gondwana both rifted along geometric patterns that closely approximate truncated-icosahedral tessellations of the lithosphere. These large-scale, quasi-hexagonal rift patterns manifest a least-work configuration. For both Laurentia and Gondwana, continental rifting coincided with drift stagnation, and may have been driven by lithospheric extension above an insulated and thermally expanded mantle. Anorogenic magmatism, including flood basalts, dike swarms, anorthosite massifs and granite-rhyolite provinces, originated along the Laurentian and Gondwanan rift tessellations. Long-lived volcanic regions of the Atlantic and Indian Oceans, sometimes called hotspots, originated near triple junctions of the Gondwanan tessellation as the supercontinent broke apart. We suggest that some anorogenic magmatism results from decompression melting of asthenosphere beneath opening fractures, rather than from random impingement of hypothetical deep-mantle plumes.     

北祁连山白山子花岗闪长岩成岩时代

寒山大型金矿是近年来在北祁连山西段发现的与侵人岩有关的构造蚀变岩型金矿。笔者首次利用锆石U-Pb方法测得白山子花岗闪长岩的形成年龄为370±25 Ma，属于华力西期。寒山金矿的成矿主要在213.95~339 Ma间。在多期热液叠加，多期成矿作用中，早期成矿的热液很可能是白山子花岗闪长岩（370±25 Ma) ,寒山辉长岩(347.1±6. 4 Ma）共同提供的。由于该区有较多的中酸性岩体存在，它们可能为金的主要来源，因而，确定这些侵入体的形成年龄，对于在该区寻找蚀变岩型金矿不仅有重要的理论意义，而且有重要的现实意义。