Matlab files used to produce the figures in the publication: On the Derivation of Closed-Form Expressions for Displacements, Strains and Stresses Inside Poroelastic Reservoirs

MATLAB code to recreate figures from the paper:<br><em>Cornelissen, P., Meulenbroek, B.J., Jansen, J.D., (2023). On the Derivation of Closed-Form Expressions for Displacements, Strains and Stresses Inside Poroelastic Reservoirs. Submitted to JGR Solid Earth.</em><br>The MATLAB scripts are used to compute the analytical solution for the poroelastic strain and stress field for rectangular and triangular reservoirs. The analytical solutions have been derived using inclusion theory (Equations 40-42 and 45-47 of the associated manuscript). The purpose of these scripts is to verify the analytical solutions by computing the strain discontinuity across the inclusion boundary and computing the volumetric strain in and outside the inclusion, and to show an example of an application of inclusion theory to reservoir depletion. These scripts produce Figures 5, 6, and 9 of the associated manuscript.<br><br><br>

用于复现以下论文图件的MATLAB代码：<em>Cornelissen, P.、Meulenbroek, B.J.、Jansen, J.D.（2023）。《孔隙弹性（poroelastic）储层内部位移、应变与应力闭合式表达式的推导》，已提交至《地球物理学研究杂志：固体地球》（JGR Solid Earth）</em>。本MATLAB脚本用于计算矩形与三角形储层的孔隙弹性应变场与应力场解析解。该解析解基于夹杂理论（inclusion theory）推导得到，对应关联手稿中的公式40-42与45-47。本脚本的设计目的为：通过计算夹杂边界处的应变不连续性（strain discontinuity），以及夹杂内部与外部的体积应变（volumetric strain），验证该解析解的正确性；同时展示将夹杂理论应用于储层耗竭（reservoir depletion）的实例。本脚本可复现关联手稿中的图5、图6与图9。