Theoretical analyses of acidization dissolution front instability in fluid‐saturated carbonate rocks

This paper presents an instability theory that can be used to understand the fundamental behavior of an acidization dissolution front when it propagates in fluid‐saturated carbonate rocks. The proposed theory includes two fundamental concepts, namely the intrinsic time and length of an acidization dissolution system, and a theoretical criterion that involves the comparison of the Zhao number and its critical value of the acidization dissolution system. The intrinsic time is used to determine the time scale at which the acidization dissolution front is formed, while the intrinsic length is used to determine the length scale at which the instability of the acidization dissolution front can be initiated. Under the assumption that the acidization dissolution reaction is a fast process, the critical Zhao number, which is used to assess the instability likelihood of an acidization dissolution front propagating in fluid‐saturated carbonate rocks, has been derived in a strictly mathematical manner. Based on the proposed instability theory of a propagating acidization dissolution front, it has been theoretically recognized that: (i) the increase of the mineral dissolution ratio can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (ii) the increase of the final porosity of the carbonate rock can destabilize the acidization dissolution front, while the increase of the initial porosity can stabilize the acidization dissolution front in fluid‐saturated carbonate rocks; (iii) the increase of the mineral dissolution ratio can cause an increase in the dimensionless propagation speed of the acidization dissolution front; (iv) the increase of the initial porosity can enable the acidization dissolution front to propagate faster, while the increase of the final porosity can enable the acidization dissolution front to propagate slower in the acidization dissolution system. Copyright © 2012 John Wiley & Sons, Ltd.     

A numerical simulation of gas flow during coal/gas outbursts

Summary A model of the gas flow in airways during an instantaneous outburst of coal and gas is formulated and solved numerically using MacCormack's explicit finite-difference scheme. This model is based on the assumption that geological structures, in-situ stresses and high-gas-pressure gradients play important roles in initiating an outburst, with the gas content and gas-pressure gradients being the most dominant factors. The fluid-dynamic processes that occur after an outburst are computed by the numerical integration of the complete time-dependent Navier-Stokes equations. The mixture velocity, the density and the gas-concentration profiles in both time and space domains (immediately after an outburst) are presented. The global results are useful in gaining an improved understanding of gas-flow patterns during coal/gas outbursts and in determining the range of the disturbance so that effective methods of control can be developed.Nomenclature a speed of sound - C gas concentration (mass) - C _v specific heat at constant volume - C _p specific heat at constant pressure - CFL Courant-Friedrichs-Lewy number - E total energy of the mixture - F the vector defined in Equation 27 - G the vector defined in Equation 28 - h entropy - h _i entropy at inlet - J _xi diffusion flux in the i-direction - k thermal conductivity of a gas - M Mach number - P pressure of a mixture (the partial pressure of gas) - P gas pressure - P _a atmospheric pressure     

Run-up from impact tsunami

We report on calculations of the on-shore run-up of waves that might be generated by the impact of subkilometre asteroids into the deep ocean. The calculations were done with the COULWAVE code, which models the propagation and shore-interaction of non-linear moderate- to long-wavelength waves  ( kh < π)  using the extended Boussinesq approximation. We carried out run-up calculations for several different situations: (1) laboratory-scale monochromatic wave trains onto simple slopes; (2) 10–100 m monochromatic wave trains onto simple slopes; (3) 10–100 m monochromatic wave trains onto a compound slope representing a typical bathymetric profile of the Pacific coast of North America; (4) time-variable scaled trains generated by the collapse of an impact cavity in deep water onto simple slopes and (5) full-amplitude trains onto the Pacific coast profile. For the last case, we also investigated the effects of bottom friction on the run-up. For all cases, we compare our results with the so-called 'Irribaren scaling': The relative run-up   R / H ₀=ξ= s ( H ₀/ L ₀)^−1/2  , where the run-up is   R , H ₀  is the deep-water waveheight, L ₀ is the deep-water wavelength, s is the slope and ξ is a dimensionless quantity known as the Irribaren number. Our results suggest that Irribaren scaling breaks down for shallow slopes   s ≤ 0.01  when  ξ < 0.1 − 0.2  , below which   R / H ₀  is approximately constant. This regime corresponds to steep waves and very shallow slopes, which are the most relevant for impact tsunami, but also the most difficult to access experimentally.     

High-resolution structures of the Landers fault zone inferred from aftershock waveform data

High-frequency body waves recorded by a temporary seismic array across the surface rupture trace of the 1992 Landers, California, earthquake were used to determine fault-zone structures down to the seismogenic depth. We first developed a technique to use generalized ray theory to compute synthetic seismograms for arbitrarily oriented tabular low-velocity fault-zone models. We then generated synthetic waveform record sections of a linear array across a vertical fault zone. They show that both arrival times and waveforms of P and S waves vary systematically across the fault due to transmissions and reflections from boundaries of the low-velocity fault zone. The waveform characteristics and arrival-time patterns in the record sections allow us to locate the boundaries of the fault zone and to determine its P - and S -wave velocities independently as well as its depth extent. Therefore, the trade-off between the fault-zone width and velocities can be avoided. Applying the method to the Landers waveform data reveals a low-velocity zone with a width of 270–360 m and a 35–60 per cent reduction in P and S velocities relative to the host rock. The analysis suggests that the low-velocity zone extends to a depth of ∼7 km. The western boundary of the low-velocity zone coincides with the observed main surface rupture trace.     

Gaussian beams in inhomogeneous media: A review

The paper outlines the most important results of the paraxial complex geometrical optics (CGO) in respect to Gaussian beams diffraction in the smooth inhomogeneous media and discusses interrelations between CGO and other asymptotic methods, which reduce the problem of Gaussian beam diffraction to the solution of ordinary differential equations, namely: (i) Babich’s method, which deals with the abridged parabolic equation and describes diffraction of the Gaussian beams; (ii) complex form of the dynamic ray tracing method, which generalizes paraxial ray approximation on Gaussian beams and (iii) paraxial WKB approximation by Pereverzev, which gives the results, quite close to those of Babich’s method. For Gaussian beams all the methods under consideration lead to the similar ordinary differential equations, which are complex-valued nonlinear Riccati equation and related system of complex-valued linear equations of paraxial ray approximation. It is pointed out that Babich’s method provides diffraction substantiation both for the paraxial CGO and for complex-valued dynamic ray tracing method. It is emphasized also that the latter two methods are conceptually equivalent to each other, operate with the equivalent equations and in fact are twins, though they differ by names. The paper illustrates abilities of the paraxial CGO method by two available analytical solutions: Gaussian beam diffraction in the homogeneous and in the lens-like media, and by the numerical example: Gaussian beam reflection from a plane-layered medium.