Coupled Simulation of Seismic Wave Propagation and Failure Phenomena by Use of an MPS Method     

Junichi Takekawa  Hitoshi Mikada  Tada-nori Goto  Yoshinori Sanada  Yuzuru Ashida 《Pure and Applied Geophysics》2013,170(4):561-570

The failure of brittle materials, for example glasses and rock masses, is commonly observed to be discontinuous. It is, however, difficult to simulate these phenomena by use of conventional numerical simulation methods, for example the finite difference method or the finite element method, because of the presence of computational grids or elements artificially introduced before the simulation. It is, therefore, important for research on such discontinuous failures in science and engineering to analyze the phenomena seamlessly. This study deals with the coupled simulation of elastic wave propagation and failure phenomena by use of a moving particle semi-implicit (MPS) method. It is simple to model the objects of analysis because no grid or lattice structure is necessary. In addition, lack of a grid or lattice structure makes it simple to simulate large deformations and failure phenomena at the same time. We first compare analytical and MPS solutions by use of Lamb’s problem with different offset distances, material properties, and source frequencies. Our results show that analytical and numerical seismograms are in good agreement with each other for 20 particles in a minimum wavelength. Finally, we focus our attention on the Hopkinson effect as an example of failure induced by elastic wave propagation. In the application of the MPS, the algorithm is basically the same as in the previous calculation except for the introduction of a failure criterion. The failure criterion applied in this study is that particle connectivity must be disconnected when the distance between the particles exceeds a failure threshold. We applied the developed algorithm to a suspended specimen that was modeled as a long bar consisting of thousands of particles. A compressional wave in the bar is generated by an abrupt pressure change on one edge. The compressional wave propagates along the interior of the specimen and is visualized clearly. At the other end of the bar, the spalling of the bar is reproduced numerically, and a broken piece of the bar is formed and falls away from the main body of the bar. Consequently, these results show that the MPS method effectively reproduces wave propagation and failure phenomena at the same time.  相似文献   

A new significance on the vertical component ratio of the power spectra between two sites in the application of array methods     

Xinrui Zhang  Hitoshi Morikawa 《Journal of Seismology》2016,20(1):185-196

  相似文献   

Possible variation in methane flux caused by gas hydrate formation on the northeastern continental slope off Sakhalin Island, Russia     

Hirotsugu Minami  Kazuya Tatsumi  Akihiro Hachikubo  Satoshi Yamashita  Hirotoshi Sakagami  Nobuo Takahashi  Hitoshi Shoji  Young Keun Jin  Anatoly Obzhirov  Nataliya Nikolaeva  Alexander Derkachev 《Geo-Marine Letters》2012,32(5-6):525-534

The Sakhalin Slope Gas Hydrate Project (SSGH) is an international collaborative effort by scientists from Japan, Korea, and Russia to investigate natural gas hydrates (GHs) that have accumulated on the continental slope off Sakhalin Island, Okhotsk Sea. From 2009 to 2011, field operations of the SSGH-09, -10, and -11 projects were conducted. GH-bearing and -free sediment cores were retrieved using steel hydro- and gravity corers. The concentrations of sulfate ions in sediment pore waters were measured to investigate sulfate concentration–depth profiles. Seventeen cores showed linear depth profiles of sulfate concentrations. In contrast, eight cores and two cores showed concave-up and -down profiles plausibly explained by sudden increase and decrease in methane flux from below, respectively, presumably caused by the formation of gas hydrate adjacent to the core sampling sites.  相似文献   

Seismic response of a cut‐and‐cover tunnel isolated by polymer material     

Satoshi Kiryu  Yoshitaka Murono  Hitoshi Morikawa 《地震工程与结构动力学》2012,41(14):2043-2057

An effective method for the seismic retrofitting of the sidewalls of cut‐and‐cover tunnels has not yet been established. Thus, a new method of seismic retrofitting for cut‐and‐cover tunnels, called the ‘polymer isolation method’, is developed here. In this method, thin walls made of polymer materials, called ‘isolation walls’, are inserted between the ground and the sidewalls of a tunnel. We demonstrate the effectiveness of the proposed method in reducing the seismic response of tunnels by using some numerical simulations. It is found that the proposed method depends on the thickness of the soil cover over the tunnel and the ratio of the stiffness between the soil and the structure. Furthermore, a simple chart is proposed for convenience to represent the applicability of the polymer isolation method to the design of seismic retrofits for cut‐and‐cover tunnels. Although this particular chart has been obtained through a limited case, some other typical cases in which the chart can also be implemented for special conditions are addressed, and the applicability and the limitations of this chart are suggested for possible conditions of the ground and of tunnel structures. Copyright © 2012 John Wiley & Sons, Ltd.  相似文献   

An efficient method to compute the dynamic response of a fluid-filled crack     

Mare Yamamoto  Hitoshi Kawakatsu 《Geophysical Journal International》2008,174(3):1174-1186

The vibration of a fluid-filled crack is considered to be one of the most plausible source mechanisms for the long-period events and volcanic tremors occurring around volcanoes. As a tool for the quantitative interpretation of source process of such volcanic seismic signals, we propose a method to numerically simulate the dynamic response of a fluid-filled crack. In this method, we formulate the motions of the fluid inside and the elastic solid outside of the crack, using boundary integrals in the frequency domain and solve the dynamic interactions between the fluid and the elastic solid using the point collocation method. The present method is more efficient compared with the time-domain finite difference method, which has been used in simulations of a fluid-filled crack and enables us to study the dynamics of a fluid-filled crack over a wide range of physical parameters. The method also allows us direct calculation of the attenuation quality factor of the crack resonance, which is an indispensable parameter for estimating the properties of the fluid inside the crack. The method is also designed to be flexible to many applications, which may be encountered in volcano seismology, and thus, extensions of the method to more complicated problems are promising.  相似文献   

Gas hydrate saturation at Site C0002, IODP Expeditions 314 and 315, in the Kumano Basin,Nankai trough     

Ayumu Miyakawa  Saneatsu Saito  Yasuhiro Yamada  Hitoshi Tomaru  Masataka Kinoshita  Takeshi Tsuji 《Island Arc》2014,23(2):142-156

The degree of gas hydrate saturation at Integrated Ocean Drilling Program (IODP) Site C0002 in the Kumano Basin, Nankai Trough, was estimated from logging‐while‐drilling logs and core samples obtained during IODP Expeditions 314 and 315. Sediment porosity data necessary for the calculation of saturation were obtained from both core samples and density logs. Two forms of the Archie equation (‘quick‐look’ and ‘standard’) were used to calculate gas hydrate saturation from two types of electrical resistivity log data (ring resistivity and bit resistivity), and a three‐phase Biot‐type equation was used to calculate gas hydrate saturation from P‐wave velocity log data. The gas hydrate saturation baseline calculated from both resistivity logs ranges from 0% to 35%, and that calculated from the P‐wave velocity log ranges from 0% to 30%. High levels of gas hydrate saturation (>60%) are present as spikes in the ring resistivity log and correspond to the presence of gas hydrate concentrations within sandy layers. At several depths, saturation values obtained from P‐wave velocity data are lower than those obtained from bit resistivity data; this discrepancy is related to the presence of free gas at these depths. Previous research has suggested that gas from deep levels in the Kumano Basin has migrated up‐dip towards the southern and seaward edge of the basin near Site C0002. The high saturation values and presence of free gas at site C0002 suggest that a large gas flux is flowing to the southern and seaward edge of the basin from a deeper and/or more landward part of the Kumano Basin, with the southern edge of the Kumano Basin (the location of site C0002) being the main area of fluid accumulation.  相似文献   

Preferred orientation of the trace fossil Entradichnus meniscus in eolian dune strata (Djadokhta Formation) at Tugrikiin Shiree,southern Mongolia and its paleoecological implications     

Koji Seike  Hitoshi Hasegawa  Niiden Ichinnorov 《Island Arc》2010,19(4):622-627

Eolian sand dune deposits of the Upper Cretaceous Djadokhta Formation at Tugrikiin Shiree, southern Mongolia, yield not only dinosaur skeletal remains but also numerous trace fossils produced by invertebrates. This paper describes the trace fossil Entradichnus meniscus, a long unlined and unbranched trail that is filled with meniscate laminae and occurs characteristically in positive epirelief. The trail is straight to gently meandering, parallel to the foreset laminae of the eolian dunes, and exhibits a significant preferred orientation parallel to the depositional dip of the cross‐stratification laminae. In addition, almost all the crescentic internal laminae of the trail show concave down‐dips. These features indicate that the trails were produced beneath the slipface of eolian dunes by the downward burrowing of the trace‐makers. This occurrence mode of E. meniscus of the Mongolian Cretaceous is very similar to that described from the Jurassic eolian dune deposits in North America. Hence, the downward burrowing of the E. meniscus animal might be a common feature in arid eolian dune deposits at least during the Jurassic and Cretaceous, and possibly reflecting a behavioral response to the morphology of large sand dunes under an arid climate.  相似文献   

Tsunami Observations in Rivers from a Perspective of Tsunami Interaction with Tide and Riverine Flow     

Elena Tolkova  Hitoshi Tanaka  Min Roh 《Pure and Applied Geophysics》2015,172(3-4):953-968

  相似文献   

Detailed Distribution of Accelerating Foreshocks before a M 5.1 Earthquake in Japan     

K. Yamaoka  T. Ooida  Y. Ueda 《Pure and Applied Geophysics》1999,155(2-4):335-353

—The M 5.1 event (May 23, 1993) which occurred in one of the most active swarm areas of Japan was preceded by foreshock activity. We obtained precise hypocenters of the foreshock-mainshock-aftershock sequences with a temporary seismic network installed just above the source region twenty days before the mainshock. The foreshocks are very unique in their accelerating activity; the acceleration in the number of foreshocks enabled us to estimate the time of the mainshock with time-to-failure analysis proposed by . Although substantial snow remained in the swarm area, we quickly installed the network because the time-to-failure analysis disclosed that the mainshock was impending. The temporary network provided detailed information on both the temporal and spatial distribution of the foreshock-mainshock-aftershock sequences. Foreshocks started fifty days before the mainshock and were distributed linearly at the base of the seismogenic layer with a length of 5 km and horizontal and vertical widths of about 1 km. The temporal change of the number of foreshocks is approximated by a power law, and the time of the mainshock can be estimated by extrapolating plots of the inverse of the daily number of events. An area of seismic quiescence appeared 40 hours before the mainshock and propagated with a rate of 20 m/hour. The mainshock occurred 2 km westward from the primary foreshock area. It was located at the base of the aftershock region. This process can be interpreted as source nucleation; preslip on the fault prior to the mainshock.  相似文献   

Isotopic composition of dissolved inorganic carbon in subsurface sediments of gas hydrate-bearing mud volcanoes, Lake Baikal: implications for methane and carbonate origin   总被引：1，自引：1，他引：0  

Alexey A. Krylov  Oleg M. Khlystov  Akihiro Hachikubo  Hirotsugu Minami  Yutaka Nunokawa  Hitoshi Shoji  Tamara I. Zemskaya  Lieven Naudts  Tatyana V. Pogodaeva  Masato Kida  Gennady V. Kalmychkov  Jeffrey Poort 《Geo-Marine Letters》2010,30(3-4):427-437

We report on the isotopic composition of dissolved inorganic carbon (DIC) in pore-water samples recovered by gravity coring from near-bottom sediments at gas hydrate-bearing mud volcanoes/gas flares (Malenky, Peschanka, Peschanka 2, Goloustnoe, and Irkutsk) in the Southern Basin of Lake Baikal. The δ¹³C values of DIC become heavier with increasing subbottom depth, and vary between ?9.5 and +21.4‰ PDB. Enrichment of DIC in ¹³C indicates active methane generation in anaerobic environments near the lake bottom. These data confirm our previous assumption that crystallization of carbonates (siderites) in subsurface sediments is a result of methane generation. Types of methanogenesis (microbial methyl-type fermentation versus CO₂-reduction) were revealed by determining the offset of δ¹³C between dissolved CH₄ and CO₂, and also by using δ¹³C and δD values of dissolved methane present in the pore waters. Results show that both mechanisms are most likely responsible for methane generation at the investigated locations.  相似文献