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.     

Periodic orbits and their bifurcations in a 3-D system

We study some simple periodic orbits and their bifurcations in the Hamiltonian . We give the forms of the orbits, the characteristics of the main families, and some existence diagrams and stability diagrams. The existence diagram of the family 1a contains regions that are stable (S), simply unstable (U), doubly unstable (DU) and complex unstable (). In the regionsS andU there are lines of equal rotation numberm/n. Along these lines we have bifurcations of families of periodic orbits of multiplicityn. When these lines reach the boundary of the complex unstable region, they are tangent to it. Inside the region there are linesm/n, along which the orbits 1a, describedn-times, are doubly unstable; however, along these lines there are no bifurcations ofn-ple periodic orbits. The families bifurcating from 1a exist only in certain regions of the parameter space (, ). The limiting lines of these regions join at particular points representing collisions of bifurcations. These collisions of bifurcations produce a nonuniqueness of the various families of periodic orbits. The complicated structure of the various bifurcations can be understood by constructing appropriate stability diagrams.     

Effect of perturbed potentials on the non-linear stability of libration pointL 4 in the restricted problem

The non-linear stability of the libration pointL ₄ in the restricted problem has been studied when there are perturbations in the potentials between the bodies. It is seen that the pointL ₄ is stable for all mass ratios in the range of linear stability except for three mass ratios depending upon the perturbing functions. The theory is applied to the following four cases:

Nonlinear stability of the libration points in the photogravitational restricted three body problem

The stability of the triangular libration points in the case when the first and the second order resonances appear was investigated. It was proved that the first order resonances do not cause instability. The second order resonances may lead to instability. Domains of the instability in the two-dimensional parameter space were determined.     

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.     

Stability of outer planetary orbits around binary stars: A comparison of Hill's and Laplace's stability criteria

A comparison is made between the stability criteria of Hill and that of Laplace to determine the stability of outer planetary orbits encircling binary stars. The restricted, analytically determined results of Hill's method by Szebehely and co-workers and the general, numerically integrated results of Laplace's method by Graziani and Black are compared for varying values of the mass parameter =m ₂/(m ₁+m ₂). For 00.15, the closest orbit (lower limit of radius) an outer planet in a binary system can have and still remain stable is determined by Hill's stability criterion. For >0.15, the critical radius is determined by Laplace's stability criterion. It appears that the Graziani-Black stability criterion describes the critical orbit within a few percent for all values of .     

涠西南地区岩石特性及钻井液技术研究

北部湾盆地涠西南油田群是南海西部重要的原油产区，同时油田群紧邻众多国家级自然保护区，属于环境敏感区域。针对北部湾盆地涠西南地区环保要求，以及涠洲组和流沙港组钻井过程经常出现井壁失稳导致的憋卡、起下钻阻卡等问题，开展了涠西南地区地层泥页岩特性及环保钻井液技术研究。通过地层岩石黏土矿物分析、孔喉结构分析以及理化性能评价，明确了涠西南地区泥页岩井壁失稳的机理；提出以类油基的水溶性复合基基液为核心，构建了一套具有油基钻井液工程特性和水基钻井液环保特性的新型环保防塌钻井液技术，并进行了相关的现场应用。现场应用显示：应用井与邻井相比，12-1/4″井段阻卡划眼时间减少，井径扩大率大幅度降低，并且在128 h的长时间浸泡过程中没有复杂情况产生，有效地解决了涠洲组、流沙港组易失稳地层的井壁稳定问题；同时，钻井液环保性能达到一级海域环保要求，可以替代目前应用的油基钻井液体系，解决了海上使用油基钻井液存在配制成本高、含油钻屑需全回收及环境危害隐患大等技术难题。该研究对于涠西南油田环保、井壁失稳和钻井液技术发展具有较高的理论价值和借鉴意义。     

施工顺序对软土地坪复合地基的影响探讨

场地回填土与搅拌桩施工顺序将直接影响软土地坪复合地基的处理效果,其中沉降控制是影响地坪安全运行的关键。以近海软土地区工业厂房水泥土搅拌桩复合地基工程为例,探讨了复合地基水泥土搅拌桩、回填土两种不同施工顺序的影响因素及相应处理效果;通过数值模拟试验,计算分析了不同填土厚度下两种工况的地基固结总沉降、施工沉降、工后沉降。结果表明,采用“先土后桩”的施工顺序的施工期沉降大于“先桩后土”,但工后沉降大大减小,提高了软土地坪地基处理效果,有利于处理后地坪的安全运行,验证了前述分析结论,为合理设计方案的确定提供了重要的参考依据。