The invention discloses a fracture-following cavern-finding
acid fracturing numerical
simulation method, which comprises the following steps: establishing a stress-flow-acid rock reaction
coupling mathematical model: establishing a stress-flow-acid rock reaction
coupling control equation, and setting initial conditions and boundary conditions; adopting a maximum
strain energy release rate criterion to judge crack propagation; combining the
control equation and the boundary condition by adopting a weighted margin method to obtain an equivalent integral weak form of the
control equation; discretizing the control equation by adopting an
extended finite element method to obtain an extended finite element discrete equation; establishing a
physical model containing natural fractures and
karst caves; and calculating the
physical model based on the extended finite element discrete equation. According to the method, a set of crack hole type
carbonate reservoir
acid etching crack propagation model comprising a stress module, a flow module and an acid rock reaction module is established, a control equation is discretized by adopting an expansion
finite element method, crack propagation is judged by utilizing a maximum
strain energy release rate criterion, and an
acid etching crack propagation path is researched.