Consolidation computation of aquitard considering non-Darcy flow

In the deep aquitard, the clay soil was stiff and compact for the diagenesis. The flow of pore water in these clay soil layers may not obey Darcy's flow. But the consolidation computation of deep aquitard still utilizes these consolidation equations that based on Darcy's flow. In order to improve the precision of consolidation computation of deep aquitard, Non-Darcy flow (described by the threshold gradient i(l) and exponent m) was incorporated in the one-dimensional consolidation equation of Terzaghi. And semi-analysis was used to solve non-Darcy flow one-dimensional consolidation equation of Terzaghi. A consolidation case was computed by Darcy and non-Darcy flow one-dimensional consolidation equation of Terzaghi respectively; and the authors found that (1) for the special degree of consolidation, the consolidation time would increase when non-Darcy flow was considered; (2) The calculation results of non-Darcy flow consolidation equation were consistent with the measured values.