人工边界高频振荡失稳机理的一点注记

进一步完善了人工边界高频振荡失稳机理的解释. 解析证明了对于一维波动模型的离散模型, 当人工波速大于1.5倍的离散空间步距与时间步距的比值时, 则在某一高频段内, 其稳态波动解在人工边界上反射系数的模大于1的命题. 在此基础上对失稳的高频段作了进一步的讨论, 揭示了高频振荡失稳发生于对波动数值模拟无意义的高频段.      

Magnetic anisotropy of polycrystalline haematite induced by A D. C. magnetic field

Summary The changes in the anisotropy of initial susceptibility of haematite ores, induced by a D.C. magnetic field are studied. After applying stronger magnetic fields, the changes in the anisotropy parameters are considerable and depend on the manner in which the basal planes of the haematite grains are arranged in the samples. The experimental results are interpreted on the basis of the single-domain theory. The proposed simple model of induced anisotropy is based on the relation between the directional susceptibility of the sample and the preferred orientation of the magnetic moments of the grains in polycrystalline haematite.

---

am uu aumnuu aa nuuumuamum ¶rt;, a nm aum n. nuu u aum n uu aumnuu aam u aum m na n¶rt;u nmamum aa. numa. ¶rt;a umnmua unau ¶rt;¶rt; muu. ¶rt;a nma ¶rt; a aumnuu aa a auumu u nuuumu m umauu aum m nuumauamum.

---

---

Presented at the 5th Scientific Assembly of IAGA in Prague 1985.     

Assessment of hydro power potential of Swat,Kohistan Himalayas: A solution for energy shortfall in region

Pakistan is facing acute power dearth these days due to which people are suffering from privation and their economic condition is deflating. Expensive imported fossil fuel for thermal power generation is pressing them hard and exacerbating the situation. The total power produced in the country by all means is less by 4000 MW than the requirement of Pakistan. So this is the prime time to think about exploring aboriginal and renewable power resources. This article converses about hydro power resources of Swat, North Pakistan, one of the most promising areas in hydro power production. The area drained by the Swat River and its tributary net work is a part of Kohistan Himalayas and has a potential of more than 1000 MW that can contribute 25% to the total energy shortfall in the country and the same is 150% more than the shortfall experienced in the province of Khyber Pakhtunkhwa that hosts Swat Valley. Due to great relief difference acquired within a small stretch, the perennial Swat River and its tributaries become pertinent for power production as high head projects and also gain importance in irrigation as well.     

Ductile creep and compaction: A mechanism for transiently increasing fluid pressure in mostly sealed fault zones

A simple cyclic process is proposed to explain why major strike-slip fault zones, including the San Andreas, are weak. Field and laboratory studies suggest that the fluid within fault zones is often mostly sealed from that in the surrounding country rock. Ductile creep driven by the difference between fluid pressure and lithostatic pressure within a fault zone leads to compaction that increases fluid pressure. The increased fluid pressure allows frictional failure in earthquakes at shear tractions far below those required when fluid pressure is hydrostatic. The frictional slip associated with earthquakes creates porosity in the fault zone. The cycle adjusts so that no net porosity is created (if the fault zone remains constant width). The fluid pressure within the fault zone reaches long-term dynamic equilibrium with the (hydrostatic) pressure in the country rock. One-dimensional models of this process lead to repeatable and predictable earthquake cycles. However, even modest complexity, such as two parallel fault splays with different pressure histories, will lead to complicated earthquake cycles. Two-dimensional calculations allowed computation of stress and fluid pressure as a function of depth but had complicated behavior with the unacceptable feature that numerical nodes failed one at a time rather than in large earthquakes. A possible way to remove this unphysical feature from the models would be to include a failure law in which the coefficient of friction increases at first with frictional slip, stabilizing the fault, and then decreases with further slip, destabilizing it.     

大型引水隧道在平面地震波入射下动力响应的解析解

本文以波函数展开法为基础,给出了SV波与P波入射时大型引水隧道平面地震响应的解析解.并根据解析解,利用数值方法对建立的场地模型进行了计算分析,研究了隧道结构-水体体系的地震响应.重点讨论了入射波频率、入射角、衬砌力学性能对衬砌结构动力响应的影响.分析表明,大型引水隧道动力响应受入射波频率和入射角的影响较为明显,而衬砌力...     

Full moment tensor inversion constrained by double-couple focal mechanism for induced seismicity*

In this study, we propose a new method to determine full moment tensor solution for induced seismicity. This method generalizes the full waveform matching algorithm we have developed to determine the double-couple (DC) focal mechanism based on the neighbourhood algorithm. One major difference between the new method and the former one is that we adopt a new misfit function to constrain the candidate moment tensor solutions with respect to a reference DC solution in addition to other misfit terms characterizing the waveform matching. Through synthetic tests using a real passive seismic survey geometry, the results show the new constraint can help better recover the DC components of inverted moment tensors. We further investigate how errors in the velocity model and source location affect the moment tensor solution. The synthetic test results indicate that the constrained inversion is robust in recovering both the DC and non-DC components. We also test the proposed method on several real induced events in an oil/gas field in Oman using the same observation system as synthetic tests. While it is found that the full moment tensor solutions without using the DC constraints have much larger non-DC components than solutions with the DC constraints, both solutions are able to fit the observed waveforms at similar levels. The synthetic and real test results suggest the proposed DC constrained inversion method can reliably retrieve full moment tensor solutions for the induced seismicity.     

Research Note: A simple method of image solution for a sphere of constant electrical potential in a conducting half‐space: implications for the applied potential method

The applied potential, or mise‐à‐la‐masse, method is used in mineral exploration and environmental applications to constrain the shape and extent of conductive anomalies. However, few simple calculations exist to help gain understanding and intuition regarding the pattern of measured electrical potential at the ground surface. While it makes intuitive sense that the conductor must come close to the ground surface in order for the lateral extent of the potential anomaly to be affected by the dimensions of the conductor rather than simply by the depth, no simple calculation exists to quantify this effect. In this contribution, a simple method of images solution for the case of a sphere of constant electrical potential in a conducting half‐space is presented. The solution consists of an infinite series where the first term is the same as the method of images solution for a point current source in an infinite half‐space. The higher order terms result from the interaction of the constant potential sphere with the no‐flux boundary condition representing the ground surface and cause the change in the shape of the potential anomaly that is of interest in the applied potential method. The calculation is relevant to applied potentials when the conductive anomaly is limited in all three space dimensions and is highly conductive. Using the derived formula, it is shown that, while the electrical potential at the ground surface caused by the sphere is affected even when the sphere is quite deep, the ratio of the potential to the current, a quantity that is more relevant to the applied potential method, is not affected until the centre of the sphere is within two radii of the ground surface. An expression for the contact resistance of the sphere as a function of depth is also given, and the contact resistance is shown to increase by roughly 45% as the sphere is moved from great depth to the ground surface.     

Development of a simplified model for pore water pressure build-up induced by cyclic loading

In this paper the formulation of a simplified model for predicting pore water pressure build-up under seismic loading is updated and applied to different soils. The model is directly based on the results of cyclic laboratory tests and it is based on the damage parameter concept, avoiding any arbitrary equivalence criterion necessary to compare the seismic demand to the cyclic strength of liquefiable soils. The model is suitable to be implemented into non-linear coupled seismic response analyses since it operates in the time domain. The analytical formulation is fully described and the calibration and the physical meaning of the model parameters are analysed in detail. Simple applications show the practical usefulness of the model with respect to other literature approaches.     

Scattering of SH waves induced by a symmetrical V-shaped canyon: a unified analytical solution

This paper reports a series solution of wave functions for two-dimensional scattering and diffraction of plane SH waves induced by a symmetrical V-shaped canyon with different shape ratios. A half-space with a symmetrical V-shaped canyon is divided into two sub-regions by using a circular-arc auxiliary boundary. The two sub-regions are represented by global and local cylindrical coordinate systems, respectively. In each coordinate system, the wave field satisfying the Helmholtz equation is represented by the separation of variables method, in terms of the series of both Bessel functions and Hankel functions with unknown complex coefficients. Then, the two wave fields are described in the local coordinate system using the Graf addition theorem. Finally, the unknown coefficients are sought by satisfying the continuity conditions of the auxiliary boundary. To consider the phase characteristics of the wave scattering, a parametric analysis is carried out in the time domain by assuming an incident signal of the Ricker type. Surface and subsurface transient responses demonstrate the characteristics and mechanisms of wave propagating and scattering.     

水轮机流道压力脉动诱发厂房振动分析

水轮机压力脉动是引起厂房内源振动的主要诱因,如何在振动系统中合理、准确地施加水力荷载的振源特性,是目前机组振动分析的主要难点.文中将压力脉动假定为简谐荷载,通过跟踪楼板的振动变化和共振校核确定机组的振源特性,采用有限元法精细模拟了水轮机流道压力脉动诱发的厂房振动问题,并结合电站特点提出振动控制标准的建议值,由此对振动进...     

Magnetic smooth zones and slope anomalies: a mechanism to explain both

Marginal magnetic smooth zones occur only along continental margins where these are oblique (not normal) to the regional fracture-zone trend.This relationship suggests a simple geometric mechanism which can generate marginal smooth zones and abyssal smooth zones during times of typical reversal history.The mechanism provides an explanation of “slope anomalies” and other smooth-zone characteristics, and casts doubt on the validity of using smooth zones as evidence for long fixed polarity intervals.     

A new theoretical solution of the effect of atmospheric pressure on water level

Ａｎｅｗｔｈｅｏｒｅｔｉｃａｌｓｏｌｕｔｉｏｎｏｆｔｈｅｅｆｆｅｃｔｏｆａｔｍｏｓｐｈｅｒｉｃｐｒｅｓｓｕｒｅｏｎｗａｔｅｒｌｅｖｅｌ赵鹏君Ｐｅｎｇ－ＪｕｎＺＨＡＯ（ＳｅｉｓｍｏｌｏｇｉｃａｌＢｕｒｅａｕｏｆＡｎｈｕｉＰｒｏｖｉｎｃｅ，Ｈｅｆｅｉ２３...     

Advances in theory of plastic mechanism control: closed form solution for MR‐Frames

In this paper new advances in the application of ‘Theory of Plastic Mechanism Control’ (TPMC) are presented. TPMC is aimed at the design of structures assuring a collapse mechanism of global type. The theory has been developed in the nineties with reference to moment‐resisting frames (MRFs) and progressively extended to all the main structural typologies commonly adopted as seismic‐resistant structural systems. In particular, the outcome of the theory is the sum of the plastic moments of the columns required, at each storey, to prevent undesired failure modes, i.e. partial mechanisms and soft‐storey mechanisms. The theory is used to provide the design conditions to be satisfied, in the form of a set of inequalities where the unknowns are constituted by the column plastic moments. This set of inequalities was originally solved by means of an algorithm requiring an iterative procedure. The advances presented in this paper are constituted by the identification of a ‘closed form solution’ and by the use of TPMC in a more systematic design approach. This result is very important, because the practical application of TPMC can now be carried out even with very simple hand calculations. The practical application of TPMC is herein presented with reference to the design of a multi‐storey frame whose pattern of yielding is validated by means of both push‐over analysis and incremental dynamic analyses. Copyright © 2014 John Wiley & Sons, Ltd.     

利用全波形匹配方法确定水力压裂诱发地震震源机制

准确确定水力压裂诱发地震的震源机制对于描述裂缝破裂类型以及工区地应力状态十分重要.本文采用全波形匹配的方法确定诱发地震震源机制解,在拟合波形的同时对实际数据P波初动极性和纵横波振幅比也进行匹配,并在此基础上发展了一种基于邻域算法分级优化确定震源机制解的新方法.新方法的优势在于它对解具有更强的约束,并且能够解决现有方法中存在的权值系数选择问题.合成数据测试结果表明新方法能够得到可靠的震源机制解.本文将该方法应用于国内某页岩气田水力压裂诱发地震监测数据,反演结果表明诱发地震的震源类型以走滑破裂为主,表明该地区的最大水平主应力大于垂直主应力.

     

Perturbation solutions for vertical infiltration problems: A uniformly valid solution

The problem of vertical infiltration is solved by means of perturbation methods. The structure of the solution exhibits uniform validity for all times and all depths. For comparison between prediction and an exact result, the Burgers equation has been selected. This equation is mildly nonlinear and hence constitutes the most unfavourable known example.     

Back front of seismicity induced by non‐linear pore pressure diffusion

Borehole fluid injections are accompanied by microseismic activity not only during but also after termination of the fluid injection. Previously, this phenomenon has been analysed, assuming that the main triggering mechanism is governed by a linear pressure diffusion in a hydraulically isotropic medium. In this context the so‐called back front of seismicity has been introduced, which allows to characterize the hydraulic transport from the spatiotemporal distribution of post‐injection induced events. However, rocks are generally anisotropic, and in addition, fluid injections can strongly enhance permeability. In this case, permeability becomes a function of pressure. For such situations, we carry out a comprehensive study about the behaviour and parametrization of the back front. Based on a model of a factorized anisotropic pressure dependence of permeability, we present an approach to reconstruct the principal components of the diffusivity tensor. We apply this approach to real microseismic data and show that the back front characterizes the least hydraulic transport. To investigate the back front of non‐linear pore‐fluid pressure diffusion, we numerically consider a power‐law and an exponential‐dependent diffusivity. To account for a post‐injection enhanced hydraulic state of the rock, we introduce a model of a frozen (i.e., nearly unchanged after the stimulation) medium diffusivity and generate synthetic seismicity. We find that, for a weak non‐linearity and 3D exponential diffusion, the linear diffusion back front is still applicable. This finding is in agreement with microseismic data from Ogachi and Fenton Hill. However, for a strong non‐linear fluid–rock interaction such as hydraulic fracturing, the back front can significantly deviate from a time dependence of a linear diffusion back front. This is demonstrated for a data set from the Horn River Basin. Hence, the behaviour of the back front is a strong indicator of a non‐linear fluid–rock interaction.     

Assessment of potential strength of an induced earthquake by using fuzzy multifactorial evaluation

The potential strength of an induced earthquake depends on a series of factors at different levels. It can be effectively assessed and predicted by using the method of fuzzy multifactorial evaluation from a fuzzy set theory. As an illustration of the above-mentioned method, this paper has applied the method to assessing the potential strength of induced earthquakes due to water reservoir.     

Detection of potential leakage pathways from geological carbon storage by fluid pressure data assimilation

One of the main concerns of geological carbon storage (GCS) systems is the risk of leakage through “weak” permeable areas of the sealing formation or caprock. Since the fluid pressure pulse travels faster than the carbon dioxide (CO₂) plume across the storage reservoir, the fluid overpressure transmitted into overlying permeable formations through caprock discontinuities is potentially detectable sooner than actual CO₂ leakage occurs. In this work, an inverse modeling method based on fluid pressure measurements collected in strata above the target CO₂ storage formation is proposed, which aims at identifying the presence, the location, and the extent of possible leakage pathways through the caprock. We combine a three-dimensional subsurface multiphase flow model with ensemble-based data assimilation algorithms to recognize potential caprock discontinuities that could undermine the long-term safety of GCS. The goal of this work is to examine and compare the capabilities of data assimilation algorithms such as the ensemble smoother (ES) and the restart ensemble Kalman filter (REnKF) to detect the presence of brine and/or CO₂ leakage pathways, potentially in real-time during GCS operations. For the purpose of this study, changes in fluid pressure in the brine aquifer overlying to CO₂ storage formation aquifer are hypothetically observed in monitoring boreholes, or provided by time-lapse seismic surveys. Caprock discontinuities are typically characterized locally by higher values of permeability, so that the permeability distribution tends to fit to a non-Gaussian bimodal process, which hardly complies with the requirements of the ES and REnKF algorithms. Here, issues related to the non-Gaussianity of the caprock permeability field are investigated by developing and applying a normal score transform procedure. Results suggest that the REnKF is more effective than the ES in characterizing caprock discontinuities.     

Analytical solution for the tidally induced Lagrangian residual current in a narrow bay

In a weakly nonlinear tidal system, the depth-averaged equations for the first-order Lagrangian residual velocity (LRV) are deduced systematically. For the case of a narrow bay, the equations are solved analytically and the results for a specific bottom profile are discussed in detail. According to the pattern of the first-order LRV, the bay can be divided into three parts, namely an inner part, a transitional zone, and an outer part. For the given depth profile, the streamline of the first-order LRV for a shorter bay is a part of that for a longer bay. The first-order LRV depends on a nondimensional parameter that combines the influences of the bottom friction coefficient, the tidal period and the averaged water depth. The form of the bottom friction also has a significant influence on the first-order LRV. The second-order LRV, i.e., the Lagrangian drift, is analytically solved and shows dependence on the initial tidal phase. The LRV differs from the Eulerian residual transport velocity both quantitatively and qualitatively. It is demonstrated that the residual currents obtained according to other definitions may cause misunderstanding of the mass transport in water exchange applications.     

Analysis of induced seismicity for stress field determination and pore pressure mapping

The focal mechanisms of some one hundred microseismic events induced by various water injections have been determined. Within the same depth interval, numerous stress measurements have been conducted with the HTPF method. When inverted simultaneously, the HTPF data and the focal plane solutions help determine the complete stress field in a fairly large volume of rock (about 15×10⁶ m³). These results demonstrate that hydraulically conductive fault zones are associated with local stress heterogeneities. Some of these stress heterogeneities correspond to local stress concentrations with principal stress magnitudes much larger than those of the regional stress field. They preclude the determination of the regional stress field from the sole inversion of focal mechanisms. In addition to determining the regional stress field, the integrated inversion of focal mechanisms and HTPF data help identify the fault plane for each for each of the focal mechanisms. These slip motions have been demonstrated to be consistent with Terzaghi's effective stress principle and a Coulomb friction law with a friction coefficient ranging from 0.65 to 0.9. This has been used for mapping the pore pressure in the rock mass. This mapping shows that induced seismicity does not outline zones of high flow rate but only zones of high pore pressure. For one fault zone where no significant flow has been observed, the local pore pressure has been found to be larger than the regional minimum principal stress but no hydraulic fracturing has been detected there.