应用综合勘查方法分析平朔安太堡三号井工矿突水原因

平朔公司安太堡三号井工矿2007年3月在首采工作面辅运顺槽掘进时,迎头冒顶突水,初期突水量达6000～8000m^3/d,将1289．16m标高以下淹没。在分析研究区地层、构造及水文地质条件的基础上,对沙沟河道进行了渗漏观测,对地表排污水、沟水和井下突水点的水质进行了分析化验,有针对性的布置了6条物探线,最后进行了钻孔验证,综合分析认为突水水源主要为大沙沟河道的侧向补给,突水通道为NW向的强风化带。提出了加大井下排水,地表对大沙沟河进行土工膜防渗处理的综合治理措施。     

努力提高三维地震勘探精度，为西部煤炭工业做出新贡献

以国家重大产业技术开发专项“西部煤炭资源高精度三维地震勘探技术”项目的由来、意义和总体研究目标为引,概括的介绍了项目依托工程中各个专项技术研究完成情况,并对非均匀介质成像技术、高精度三维地震静校正技术、高密度采集技术、特观技术、岩性反演技术、属性体解释技术等六项重大关键技术取得的突破性进展进行了重点说明。指出随着我国煤炭生产重点的逐步西移,应加强诸如叠前、叠后深度偏移技术的研究,以解决复杂山区三维地震面元内地震反射波散射问题,提高其三维地震勘探精度,为西部煤炭工业做出新贡献!     

Stability Analysis and Numerical Simulation of Differential Frost Heave

Differential frost heave is often implicated in the formation of patterned ground in regions subject to recurrent freezing and thawing. A linear stability analysis (LSA) indicates that a continuum model of frost heave is linearly unstable under typical natural freezing conditions of silty-clay soils. A two-dimensional non-linear numerical analysis corroborates the frozen time LSA results, and also indicates the importance of non-linear and time-dependent terms that ultimately lead to a preferred mode, which the LSA fails to predict. Instability of the one-dimensional solution occurs at shallow freezing depths and near-zero surface loads when positive perturbations in the ice content at the freezing front lead to a concomitant increase in thermomolecular pressure and upward ice velocity. Differential frost heave can then occur because of the increased heat flux from the perturbed surfaces. A three-dimensional model using random initial surface perturbations indicates that regular surface patterns will evolve with a length scale in the order of 2–4 meters, which corresponds quite closely with naturally-occurring non-sorted patterned ground.     

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.     

Accelerating ground warming in the Canadian Prairie provinces: Is it a result of global warming?

The results of precision temperature logs made to depths of several hundred meters in some 80 wells in Western Canada, most of which are located in the Prairie Provinces, show evidence of warming at the ground surface in the 0.5 K to 3.5 K range (average=2.2±0.7 K, for 80 unevenly distributed sites). Modeling shows that this warming mostly pertains to this century and it has been most substantal in the last four decades if the ramp function of the linear increase of surface temperature is assumed. Using the step function model's increase of surface temperature (land clearing, forest fires, etc.) the calculated onset of warming would pertain mostly to the last two decades. Contour maps of ground temperatures currently and previously and a contour map of the ground warming magnitude dilineate a large regional character of the ground temperature change at the southern marigin of permafrost for the large area of the Prairie Provinces. In many cases however, the magnitude of ground warming is much larger than the magnitude of air warming. This is especially evident for the northern areas of Alberta in the boreal forest ecoprovince. The magnitude of ground warming is equal to the magnitude of surface air warming in southern Alberta in the grassland and aspen parkland ecoprovinces. The analysis of the temperature depth response to the surface warming from well data shows the integrated effect of surface air warming together with the increases in ground temperature due to natural terrain effects and other anthropogenical changes to the surface of the earth.     

Blind prediction of the site effects at Ashigara Valley,Japan, and its comparison with reality

Weak and strong ground motions were numerically predicted for three stations of the Ashigara Valley test site. The prediction was based on the records from a rock-outcrop station, one weak-motion record from a surface-sediments station, and the standard geotechnical model. The data were provided by the Japanese Working Group on the Effects of Surface Geology as a part of an international experiment. The finite-difference method for SH waves in a 2-D linear viscoelastic medium (a causalQ model) was employed.Comparison with the real records shows that at two stations the predictions fit better than at the third one. Strangely, the two better predictions were for stations situated at larger distances from the reference rock station (one station was on the surface, the other in a borehole). The strong ground motion (the peak acceleration of about 200 cm s^–2) was not predicted qualitatively worse than the weak motion (8 cm s^–2). A less sophisticated second prediction (not submitted during the experiment), in which we did not attempt to fit the available weak-motion record at the sedimentary station, agrees with the reality significantly better.     

The emplacement of intermediate volume ignimbrites: A case study from Roccamonfina Volcano,Southern Italy

A model is presented for the emplacement of intermediate volume ignimbrites based on a study of two 6 km³ volume ignimbrites on Roccamonfina Volcano, Italy. The model considers that the flows were slow moving, and quickly deflated from turbulent to non-turbulent conditions. Yield strength and density increased whereas fluidisation decreased with time and runout of the pyroclastic flows. In proximal locations, on the caldera rim, heterogeneous exposures including discontinuous lithic breccias, stratified and cross-stratified units interbedded with massive ignimbrite suggest deposition from turbulent flows. In medial locations thick, massive ignimbrite occurs associated with three types of co-ignimbrite lithic breccia which we interpret as being emplaced by non-turbulent flows. Multiple grading of different breccia/lithic concentration types within single flow units indicates that internal shear occurred producing overriding or overlapping of the rear of the flow onto the slower-moving front part. This overriding of different parts of non-turbulent pyroclastic flows could be caused by at least two different mechanisms: (1) changes in flow regime, such as hydraulic jumps that may occur at breaks in slope; and (2) periods of increased discharge rate, possibly associated with caldera collapse, producing fresh pulses of lithic-rich material that sheared onto the slower-moving part of the flow in front.We propose that ground surge deposits enriched in pumice compared with their associated ignimbrite probably formed by a flow separation mechanism from the top and front of the pyroclastic flow. These turbulent clouds moved ahead of the non-turbulent lower part of the flow to form stratified pumice-rich deposits. In distal regions well-developed coarse, often clast-supported, pumice concentrations zones and coarse intra-flow-unit lithic concentrations occur within the massive ignimbrite. We suggest that the flows were non-turbulent, possessed a relatively high yield strength and may have moved by plug flow prior to emplacement.     

The arrangement in mean elements space of the periodic orbits close to that of the Moon

In a simplified model of the Earth-Moon-Sun system based on the restricted circular 3-dimensional 3-body problem, it is possible to find numerically a set of 8 periodic orbits whose time evolutions closely resemble that of the Moon's orbit. These orbits have a period of 223 synodic months (i.e. the period of the Saros cycle known for more than two millennia as a means of predicting eclipses), and are characterized by a secular rotation of the argument of perigee . Periodic orbits of longer durations exhibiting this last feature are very abundant in Earth-Moon-Sun dynamical models. Their arrangement in the space of the mean orbital elements- for various values of the lunar mean motion is presented.     

Generalized Lyapunov exponents indicators in Hamiltonian dynamics: An application to a double star system

The Lyapunov characteristic numbers (LCNs) which are defined as the mean value of the distribution of the local variations of the tangent vectors to the flow (=ln _k ⁱ ) (see Froeschlé, 1984) have been found to be sensitive indicators of stochasticity. So we computed the distribution of these local variations and determined the moments of higher order for the integrable and stochastic regions in a binary star system with =0.5.