Uncertainty quantification of overpressure buildup through inverse modeling of compaction processes in sedimentary basins

This study illustrates a procedure conducive to a preliminary risk analysis of overpressure development in sedimentary basins characterized by alternating depositional events of sandstone and shale layers. The approach rests on two key elements: (1) forward modeling of fluid flow and compaction, and (2) application of a model-complexity reduction technique based on a generalized polynomial chaos expansion (gPCE). The forward model considers a one-dimensional vertical compaction processes. The gPCE model is then used in an inverse modeling context to obtain efficient model parameter estimation and uncertainty quantification. The methodology is applied to two field settings considered in previous literature works, i.e. the Venture Field (Scotian Shelf, Canada) and the Navarin Basin (Bering Sea, Alaska, USA), relying on available porosity and pressure information for model calibration. It is found that the best result is obtained when porosity and pressure data are considered jointly in the model calibration procedure. Uncertainty propagation from unknown input parameters to model outputs, such as pore pressure vertical distribution, is investigated and quantified. This modeling strategy enables one to quantify the relative importance of key phenomena governing the feedback between sediment compaction and fluid flow processes and driving the buildup of fluid overpressure in stratified sedimentary basins characterized by the presence of low-permeability layers. The results here illustrated (1) allow for diagnosis of the critical role played by the parameters of quantitative formulations linking porosity and permeability in compacted shales and (2) provide an explicit and detailed quantification of the effects of their uncertainty in field settings.RésuméCette étude illustre une procédure conduisant à une analyse préliminaire du risque de développement de surpression dans un bassin sédimentaire caractérisé par l’alternance de niveaux de grès et de schiste argileux. L’approche repose sur deux éléments clef : (1) une modélisation préalable des écoulements et de la compaction et (2) une application d’une technique de réduction de la complexité du modèle basée sur un développement généralisé du polynôme du chaos (gPCE). Le modèle développé considère des processus de compaction verticale unidimensionnelle. Le modèle gPCE est ensuite utilisé en situation de modélisation inverse pour obtenir une estimation fiable des paramètres du modèle et une quantification de l’incertitude. La méthodologie est appliquée à deux cas de terrain décrits dans des publications antérieures, le Venture Field (Scotian Shelf, Canada), et le Basin de Navarin (Mer de Bering, Alaska, USA), reposant sur des données de porosité et de pression disponibles pour le calibrage du modèle. Le meilleur résultat est obtenu lorsque les données de porosité et de pression sont considérées simultanément dans la procédure de calibrage du modèle. La propagation de l’incertitude des paramètres d’entrée non connus sur les données de sortie du modèle, telles que la distribution verticale de la pression dans les pores, est étudiée et quantifiée. Cette stratégie de modélisation permet de quantifier l’importance relative du phénomène clef régissant l’interaction entre la compaction du sédiment et les processus d’écoulements, et conduisant à un excès de surpression hydraulique dans un bassin sédimentaire stratifié caractérisé par la présence de niveaux de faible perméabilité. Les résultats illustrés ici (1) permettent le diagnostic du rôle critique joué par les paramètres des équations reliant la porosité et la perméabilité dans des schistes argileux compactés et (2) fournissent une quantification explicite et détaillée des effets de leur incertitude pour des conditions naturelles de terrain.RESUMENEste estudio ilustra un procedimiento apropiada para el análisis preliminar del riesgo de desarrollo de sobrepresión en cuencas sedimentarias que se caracterizan por la alternancia de eventos de deposición de capas de areniscas y lutitas. El enfoque se basa en dos elementos fundamentales: (1) el modelado directo del flujo del fluido y la compactación, y (2) la aplicación de una técnica de reducción de la complejidad del modelo basado en una expansión polinómica generalizada del caos (gPCE). El modelo directo considera procesos unidimensionales de compactación vertical. El modelo gPCE se utiliza en un contexto de modelado inverso para obtener la estimación de parámetros de un modelo eficiente y la cuantificación de la incertidumbre. La metodología se aplica a dos escenarios de campo considerados en la literatura de trabajos anteriores, es decir, Venture Field (Scotian Shelf, Canada) y la Cuenca del Navarín (Mar de Bering, Alaska, EEUU) basándose en la información de la porosidad y la presión disponible para la calibración del modelo. Se encontró que se obtiene un mejor resultado cuando los datos de porosidad y de presión se consideran conjuntamente en el procedimiento de calibración del modelo. Se investigaron y cuantificaron la incertidumbre de propagación de los parámetros desconocidos de entrada en los resultados del modelo, y la distribución vertical de la presión poral. Esta estrategia de modelado permite cuantificar la importancia relativa de los fenómenos claves que regulan la retroalimentación entre los procesos de flujo de fluido en la compactación de los sedimentos y la conducción y acumulación del exceso de presión de fluido en cuencas sedimentarias estratificadas que se caracterizan por la presencia de capas de baja permeabilidad. Los resultados ilustran aquí por (1) permitir el diagnóstico del papel fundamental desempeñado por los parámetros en las formulaciones cuantitativas que enlazan la porosidad y la permeabilidad en las lutitas compactadas y (2) proporcionar una cuantificación explícita y detallada de los efectos de la incertidumbre en la configuración de campo.摘要本研究详细阐明了砂岩和页岩层交替沉积的沉积盆地中超压发展的整个过程,这个过程有利于对超压发展进行初步风险分析。该方法依赖于两个关键要素:(1)流体流动和压实的正演模拟,及(2)基于概化多项式混沌膨胀(gPCE)的模型-复杂度还原技术。正演模型考虑了一维垂直压实过程。然后概化多项式混沌膨胀(gPCE)模型用于反演模拟环境中,以获取有效的模型参数估算结果和不确定性量化结果。该方法应用于过去文献中考虑到的两个野外情况下,即(加拿大Scotian大陆架)的Venture牧场及(美国阿拉斯加州白令海)的Navarin盆地,依靠可现有的孔隙率和压力信息进行模型校准。发现在模型校准程序中,综合考虑孔隙率和压力就能获得最好的结果。调查和量化了模型输出中未知输入参数诸如孔隙压力垂直分布导致的不确定性传播。这个模拟策略能使人们量化关键现象的相关重要性,这些关键现象控制着沉积压实和液体流动过程之间的反馈及推进存在低渗透性地层的分层沉积盆地中液体超压的逐渐积累。这里描述的的结果:(1)能够诊断连接压实页岩中孔隙率和渗透性的定量公式参数所化发挥的关键作用以及 (2)提供野外背景下不确定性影响的明确和详细的量化结果。ResumoEste estudo ilustra um procedimento útil para uma análise preliminar do risco de sobrepressão desenvolvido em bacias sedimentares caracterizadas por eventos alternados de deposição de camadas de arenitos e xisto. A abordagem apoia-se em dois elementos chave: (1) modelagem direta do fluxo de fluido e compactação; e (2) aplicação de uma técnica de redução de complexidade do modelo baseado em uma expansão polinomial generalizada do caos (qPCE). A modelagem direta considera um processo de compactação vertical unidimensional. O modelo qPCE é então usado em um contexto de modelagem inversa para obter um modelo eficiente de estimativa de parâmetros e quantificação de incertezas. A metodologia é aplicada em duas configurações de campo consideradas em obras de literatura anteriores, por exemplo, o Campo Venture (plataforma Escocesa, Canadá) e a Bacia Navarin (Mar de Bering, Alasca, EUA), contando com informações de porosidade e pressão disponíveis para calibração do modelo. Verificou-se que o melhor resultado é obtido quando os dados de porosidade e pressão são considerados em conjunto no processo de calibração do modelo. Propagação de incerteza para entrada de parâmetros desconhecidos para modelar os resultados, tais como a distribuição vertical da pressão dos poros, é investigada e quantificada. Esta estratégia possibilita uma modelagem para quantificar a importância relativa dos fenômenos chave que regulam a retroalimentação entre compactação do sedimento e o processo de fluxos de fluido e levam ao acúmulo de fluido em sobrepressão em bacias sedimentares estratificadas caracterizadas pela presença de camadas pouco permeáveis. Os resultados aqui ilustrados (1) permitem diagnósticos do papel crítico desempenhado pelos parâmetros de formulação quantitativos que ligam a porosidade e permeabilidade em xistos compactados e (2) fornecem uma quantificação explicita e detalhada dos efeitos da sua incerteza em configurações de campo.