Abstract We report a new, open-source, MATLAB-based 3D code for numerically simulating the self-potential (SP) in subsurface reservoirs. The code works as a post-processor, using outputs from existing reservoir flow and transport simulators at a selected timestep to calculate the SP throughout the reservoir model. The material properties required to calculate the SP are user defined and may be constant or vary in each cell. The code solves the equations governing flow and transport of electrical charge and global charge conservation using a control-volume-finite-difference scheme. Electrical currents associated with the SP may spread beyond the reservoir model domain, and the code allows for the domain to be extended vertically and laterally to account for this. Here, we present the governing equations and the numerical method used and demonstrate application of the code using an example in which we predict the SP signals associated with oil production from a subsurface reservoir supported by water injection.

中文翻译：

---

地下储层自电位数值模拟

摘要 我们报告了一种新的、开源的、基于 MATLAB 的 3D 代码，用于数值模拟地下储层中的自电势 (SP)。该代码用作后处理器，在选定的时间步使用现有油藏流量和输运模拟器的输出来计算整个油藏模型的 SP。计算 SP 所需的材料属性是用户定义的，并且在每个单元格中可能是恒定的或变化的。该代码使用控制体积有限差分方案求解控制电荷流动和传输以及全局电荷守恒的方程。与 SP 相关的电流可能会扩散到储层模型域之外，代码允许将域垂直和横向扩展以解决这一问题。这里，