Éditeur(s) : HAL CCSD Seismological Society of America
Résumé : International audience
We present an improvement of the maximum intersection (MAXI) method in which absolute earthquake location is defined by the maximum intersection number of hyperbolic equal differential time (EDT) volumes (one EDT volume being described as all grid nodes satisfying the arrival‐time differences between two stations, ± a tolerance value known as TERR). This 3D technique is well adapted to a strongly heterogeneous environment, avoids the depth versus origin‐time trade‐off (even using P arrivals alone), and objectively filters possible erroneous arrival times. Improvements consist of introducing an iterative multiscale approach on the TERR parameter and conducting a hierarchy of grid discretization to refine the preliminary solution (multigrid algorithm). The outlier cleaning process is strengthened by conducting statistical analyses based on the number of EDT volume intersections per station, rather than on travel‐time residues. Processing two series of synthetic arrival times computed in the 3D seismological configuration of the southernmost Ryukyu subduction zone, offshore Taiwan, allows us to (1) confirm MAXI robustness, even when arrival times are highly disturbed or azimuthal coverage largely reduced, and (2) show that MAXI, by limiting trade‐offs between origin time and depth and between epicentral position and depth, is less affected by a wrong velocity model than a technique that minimizes travel‐time residues. To conclude, we apply the MAXI method to locate earthquakes recorded by the combined Central Weather Bureau of Taiwan (CWB) and the Japanese Meteorological Agency (JMA) networks (1991-2008) within a georealistic a priori 3D P‐velocity model. Hypocenter determinations are selected based on MAXI confidence factors, and earthquake clustering along well‐known tectonic features witnesses hypocenter determination accuracy.
ISSN: 0037-1106

hal-00725227
https://hal.archives-ouvertes.fr/hal-00725227
DOI : 10.1785/0120100311

Éditeur(s) : HAL CCSD Institute for Scientific Computing and Information
Résumé : International audience
In this paper, we solve a problem describing the mechanical deformations of a porous medium in the presence of a monophasic linear flow or a two phase nonlinear flow with the purpose of modelizing subsidence of hydrocarbon reservoirs. An essential characteristics of this problem is that the mechanical deformation and the flow have different time scales. In petroleum industry, one uses different very efficient simulators for the flow problem and the mechanical deformations, which enables to handle complex models. Therefore it is necessary to be able to combine as efficiently as possible the exploitation of these simulators. We propose two alternative splitting approaches. The first one is the staggered algorithm used by engineers, which amounts to a Gauss-Seidel method in the one phase linear case. The second approach is based upon the preconditioned conjugate gradient method. We use a numerical multi-scale method in both of these algorithms. We compare these two approaches and we show that the preconditioned conjugate gradient algorithm is faster and more robust than the staggered algorithm. Applying the preconditioned conjugate gradient algorithm therefore seems to compensate for the fact that the inf-sup condition for the mixed discretization method is not satisfied when combining the simulators for the mechanical deformations and for the flow computations.
ISSN: 1705-5105

Droits : http://creativecommons.org/licenses/by/
hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039
https://hal.archives-ouvertes.fr/hal-01449039/document
https://hal.archives-ouvertes.fr/hal-01449039/file/DHLE.pdf

Éditeur(s) : HAL CCSD Institute for Scientific Computing and Information
Résumé : International audience
The production of oil in soft highly compacted reservoirs produces a significant reduction in the pore space which in turn causes the surface to subside. The authors couple two simulators to model the mechanical deformations in the presence of a Darcy flow, both one- and two-phase. The simulators solve for flow in a porous medium using a finite volume method, and solve for the mechanical deformation of the porous media using a finite element method. The flow equations and the deformation equations are coupled using Biot's law. The authors use two different methods for the coupling: the first uses a staggered algorithm; the second is based on a preconditioned conjugate method. The authors use a multiscale method in both approaches. The authors show that the inf-sup condition is not satisfied, because of the mixed discretization which is done by using two separate simulators. The authors give convergence comparisons for the two methods for several examples. Two of the examples are for the linear problem (one-phase) where they know the exact solution and one is for the nonlinear problem (two-phase). In all cases, the preconditioned conjugate algorithm is faster and more robust, in particular when the inf-sup condition is not satisfied. In the nonlinear problem the authors state that the results agree with physical predictions but do not give actual empirical data. Finally, they show the effects of grid orientation on the nonlinear problem.
ISSN: 1705-5105

hal-00730673
https://hal.archives-ouvertes.fr/hal-00730673

Éditeur(s) : HAL CCSD
Résumé : International audience
Half away between edge detection and object recognition this study aim to characterize singularities in order to guide this two phases. The caracterization is obtained with a generalization of the Lipschitz exponents .for complex structures. This new characterization is called " maxima chain " and integrates the notion of value and spatial repartition of the maxima. The maxima chains are obtained with a multiscale decomposition using a Diadic Discrete Wavelet Transform (DDWT). We use the maxima of the gradient image computed with the wavelet details. The local maxima detect the position of irregular structures (edges) in the image (Canny algorithm). Then we characterize each maxima using their evolution across the different scales. A classification is made using the maxima chains. We integrate the classification to eliminate unwanted maxima (from noise or other objects) to drive the search of closed contours. In fact, there is a totally different signature from one maxima to another depending on the parameters of the objects : nature, shape, size, gray scale, texture, ...
RFIA Proceedings
Toulouse, France

hal-00602265
https://hal.archives-ouvertes.fr/hal-00602265