Extended Coulomb failure criteria for the Zipingpu reservoir and Longmenshan slip

July 13, 2011
This figure is a presentation of viscous stress and Reynolds stress at multiple micro-spatial scales under ultra-high temperature and pressure conditions. Credit: ©Science China Press

The extended Coulomb failure stress (ECFS) criteria and anisotropic porosity and permeability tensor at micro/meso/macro scale under ultra‑high temperature and pressure (UTP) conditions were developed employing the flow driven pore‑network crack (FDPNC) model under multiple temporal–spatial scales and the hybrid hypersingular integral equation‑lattice Boltzmann method (HHIE‑LBM). The correlation of the Zipingpu reservoir and Longmenshan slip was then analyzed and the fluid–solid coupled three‑dimensional facture mechanism of the reservoir and earthquake fault was explored.

Describing the correlation of a reservoir and earthquake slip remains a great challenge. Professor SHI Yaolin and his group at the Key Laboratory of Computational Geodynamics of the Chinese Academy of Sciences and the College of Science of the Graduate University of the Chinese Academy of Sciences are attempting to solve this problem. On the basis of the laboratory's parallel central processing unit (CPU) and graphics processing unit (GPU) hardware platform, the micro FDPNC model for UTP conditions was established and the ECFS criteria and anisotropic porosity and permeability , which includes pore pressure, pressure gradient, viscous and Reynolds stress, have been given. In the case of the Zipingpu reservoir and Longmenshan fault, the simulation results show that the extended viscous stress and Reynolds stress cannot be neglected. Their work, entitled "Application of flow driven pore-network crack model to Zipingpu reservoir and Longmenshan slip", was published in SCIENCE CHINA Physics, Mechanics & Astronomy, 2011, Vol. 54(8).

The interaction mechanism of a reservoir and fault slip has importance in terms of the economy, social security and other national issues as it relates to the effects of hydro-projects, landslides and reservoir earthquakes. Several approaches have been taken to study the mechanism owing to its complexity and social importance. The use of Coulomb failure theory is a classical method for addressing this problem and has provided many landmark achievements. However, it is not clear yet how water affects a fault slip, especially at a depth of about 10 km or more, which is an important consideration. It is suggested that Coulomb criteria should be improved by consideration of the viscous stress and Reynolds stress in addition to pore pressure and its gradient.

This figure shows a presentation of the multi-spatial-scale FDPNC model. Credit: ©Science China Press

At this time, ultra-large high-performance computing hardware, parallel CPU & GPU simulation technology, fluid–solid coupled crack initiation and propagation theory and the multi-spatial-scale fluid FDPNC method provide research opportunities.

First, using micro-tomography data, the multi-spatial-scale flow-driven pore-network crack model was established, and the effects of viscous stress and Reynolds stress under UTP were determined for the first time (Figure 1).

Second, the hybrid hypersingular integral equation (HHIE) method and lattice Boltzmann (LBM) were combined to establish a link between the flow pore pressure and solid skeleton structure deformation/strength. The HHIE method deals with the crack initiation, propagation and transmission problem for the solid skeleton on multiple temporal–spatial scales (yellow in Figure 2 indicates the solid skeleton core), and the LBM deals with the flow-driven fluid pore problem on multiple temporal–spatial scales (black in Figure 2 indicates flow in pores). Combining the two methods, it is possible to study the coupled extended stress of the micro/meso/ structure at the interface of the pore and core by establishing the relationship between the extended pore fluid pressure and extended core stress, and analyze porosity and permeability under different initial and boundary conditions.

Finally, the relationship among the water-drainage sluice processes of a reservoir, stress triggers, earthquake shadows and porosity variability in a fault slip zone has been analyzed. The vadose energy (which results from pore pressure and can flow to the fault slip tip) differs for undrained and drained zones, with more energy being released in a drained zone than in an undrained zone. If the fault slip is a stable creep rupture process, the energy criteria (strain energy function factors) must increase with the rate at which the fault spreads. When penetration reaches a stable stage, the role of the initiation and propagation of the fluid FDPNC becomes dominant. As the time scale increases, the solid–fluid interface weakens and blurs microscopically, while macroscopically, the porosity increases. Strain energy is released during the slip process, and decreases as spreading increases in the drained zone.

The reservoir loading and earthquake trigger relationship depends on the fault-slip geometry and character, the porosity variability of the surrounding geological structure and the temporal and spatial scales. In the case of the Zipingpu reservoir and 2008 Wenchuan earthquake, the and temporal scale are important factors, and they need to be investigated in further field and laboratory experiments.

Explore further: Why do earthquakes stop?

More information: Zhu B J, Liu C, Shi Y L* et al. Application of flow driven pore-network crack model to Zipingpu reservoir and Longmenshan slip. SCI CHINA Physics, Mechanics & Astronomy, 2011, 54(8): 1-9

Related Stories

Why do earthquakes stop?

February 6, 2008

The underlying structure of a fault determines whether an earthquake rupture will jump from one fault to another, magnifying its size and potential devastation. Understanding why some earthquakes terminate along a fault, ...

Sun and moon trigger deep tremors on San Andreas Fault

December 23, 2009

The faint tug of the sun and moon on the San Andreas Fault stimulates tremors deep underground, suggesting that the rock 15 miles below is lubricated with highly pressurized water that allows the rock to slip with little ...

New way to monitor faults may help predict earthquakes

October 1, 2009

Scientists at the Carnegie Institution have found a way to monitor the strength of geologic faults deep in the Earth. This finding could prove to be a boon for earthquake prediction by pinpointing those faults that are likely ...

Sumatra earthquake three times larger than originally thought

February 8, 2005

Northwestern University seismologists have determined that the Dec. 26 Sumatra earthquake that set off a deadly tsunami throughout the Indian Ocean was three times larger than originally thought, making it the second largest ...

Recommended for you

Carbon coating gives biochar its garden-greening power

October 20, 2017

For more than 100 years, biochar, a carbon-rich, charcoal-like substance made from oxygen-deprived plant or other organic matter, has both delighted and puzzled scientists. As a soil additive, biochar can store carbon and ...

Cool roofs have water saving benefits too

October 20, 2017

The energy and climate benefits of cool roofs have been well established: By reflecting rather than absorbing the sun's energy, light-colored roofs keep buildings, cities, and even the entire planet cooler. Now a new study ...

1 comment

Adjust slider to filter visible comments by rank

Display comments: newest first

not rated yet Jul 13, 2011
That doesn't seem like a statement intended for release to the public and reprint in a place like this.

That sounds like a straight up quotation of a study abstract to me. I wouldn't post that here.

Please sign in to add a comment. Registration is free, and takes less than a minute. Read more

Click here to reset your password.
Sign in to get notified via email when new comments are made.