Modeling the growth and interaction of stylolite networks, using the discrete element method for pressure solution

Pressure solution (also termed chemical compaction) is considered the most important ductile deformation mechanism operating in the Earth's upper crust. This mechanism is a major player in a variety of geological processes, including evolution of sedimentary basins, hydrocarbon reservoirs, aquifers, earthquake recurrence cycles, and fault healing. Pressure solution in massive rocks often localizes into solution seams or stylolites. Field observations of stylolites often show elastic/brittle interactions in regions between pressure solution features, including and shear fractures, veins and pull-apart features. To understand these interactions, we use a grain-scale model based on the Discrete Element Method that allows granular dissolution at stressed contacts between grains. The new model captures both the slow chemical compaction process and the more abrupt brittle fracturing and sliding between grains. We simulate a sample of rock as a collection of particles, each representing either a grain or a unit of rock, bonded to each other with breakable cement. We apply external stresses to this sample, and calculate elastic and frictional interactions between the grains. Dissolution is modeled by an irreversible penetration of contacting grains into each other at a rate that depends on the contact stress and an adjustable rate constant. Experiments have shown that dissolution rates at grain contacts are greatly enhanced when there is a mineralogical contrast. Therefore, we dissolution rate constant can be increased to account for an amount of impurities (e.g. clay in a quartz or calcite sandstone) that can accumulate on dissolving contacts. This approach allows large compaction and shear strains within the rock, while allowing examination of local grain-scale heterogeneity. For example, we will describe the effect of pressure solution on the distribution of contact forces magnitudes and orientations. Contact forces in elastic granular packings are inherently heteregeneous, but stress-dependent dissolution tends to equalize them. We apply our model to the simulation of stylolite networks, particularly the interaction of stylolite tips. The stress concentrations from these tips are transmitted through the intervening rock, which can cause elastic strain, brittle damage and frictional sliding. Our model shows that grain rearrangement and compaction rate depend on the surface friction coefficient of grains. Simulation results show the development of shear zones between stylolites, and a high porosity process zone at the tips of stylolites. These features, which have been observed in field studies, are modeled and predicted for the first time. This modeling tool holds a promise to provide many new insights regarding the coupling between pressure solution and brittle deformation, i.e. between mechanical and chemical compaction.