Well stimulation

1. A method of forming a wormhole in a porous medium, comprising:

running ( 20 ) a stimulation simulator comprising:

gridding ( 22 ) a treatment region of the porous medium into a plurality of cells comprising a first portion designated as matrix cells and a second portion designated as wormhole cells;

populating ( 24 ) the simulator with static properties of the porous medium and reaction kinetic properties for reaction of the porous medium with a reactant in a treatment fluid;

modeling ( 26 ) the matrix cells wherein a medium of the matrix cells comprises matrix material behaving as a single permeability, single porosity system;

modeling ( 28 ) the wormhole cells in a wormhole initiation stage wherein a medium of the respective wormhole initiation stage cells has a solid saturation above a respective critical solid saturation and is comprised of the matrix material behaving as a single permeability, single porosity system;

modeling ( 30 ) at least a portion of the wormhole cells in a wormhole growth stage wherein the respective wormhole cells have a solid saturation equal to or less than the respective critical solid saturation, and wherein the wormhole growth stage cells comprise two different interconnected media comprised respectively of the matrix material and a wormhole material having a fluid mobility as a function of solid saturation; and

obtaining ( 32 ) optimized treatment fluid injection parameters; and

injecting ( 34 ) the treatment fluid into the treatment region of the porous medium according to the optimized treatment fluid injection parameters to form the wormhole.

2. The method of claim 1 , wherein the stimulation simulator uses a finite difference numerical method.

3. The method of claim 1 , wherein the stimulation simulator accounts for the presence in the treatment region of a multicomponent fluid selected from the group consisting of gas, aqueous and oil phases.

4. The method of claim 1 , wherein the stimulation simulator accounts for the presence in the treatment region of a plurality of solid phases.

5. The method of claim 1 , wherein the treatment region comprises a subterranean formation comprising calcium carbonate rock and the treatment fluid comprises acid delivered to the treatment region through a wellbore penetrating the subterranean formation.

6. The method of claim 1 , wherein the fluid mobility as a function of solid saturation is specified independently for each cell to characterize different behaviors of different rock types in the respective cells.

7. The method of claim 1 , wherein the wormhole initiation stage modeling accounts for dissolution of the matrix material to increase permeability and pore volume in the respective cells.

8. The method of claim 1 , wherein the media of the wormhole cells in the wormhole initiation stage modeling comprise the matrix material and the wormhole material, and the wormhole initiation stage modeling further comprises assigning values to a matrix-fracture coupling transmissibility multiplier (Sigma) such that the reactant in the treatment fluid does not interact with the wormhole material.

9. The method of claim 1 , wherein the media of the cells in the wormhole initiation stage modeling comprise the matrix material and the wormhole material, and wherein the wormhole initiation stage further comprises assigning initial values to a matrix-fracture coupling transmissibility multiplier (Sigma) such that reactant in the treatment fluid does not interact with the material of the wormhole material, and further comprising transitioning to the wormhole growth stage modeling by increasing the matrix-fracture coupling transmissibility multiplier above the respective initial values.

10. The method of claim 1 , wherein reaction of the treatment fluid with the matrix material and, where the solid saturation is equal to or less than the respective critical solid saturation, with the wormhole material, is independently parameterized to account for dissolution of the respective material(s) in the respective cells.

11. The method of claim 1 , wherein a reaction rate R r between the treatment fluid and a solid material in the cells is given by:

R r V b ·A r ·Πc ri n ri ·ΠD mijk

wherein V b is bulk volume of the respective cell, A r is a reaction rate constant, c ri is the product of reactant and solid concentrations, n ri is the order of each concentration term, and D mijk is an equilibrium deviation reaction term given by:

D mijk =θ·( F k ( a i )− C a )

wherein θ is porosity of the respective cell, F k (a i ) is a function of the reactant concentration and C a is the solid concentration.

12. The method of claim 1 , wherein the stimulation simulator comprises a table of mobility function versus solid saturation.

13. The method of claim 1 , comprising calibrating the stimulation simulator using experimental data derived from a specimen representing rock from the treatment region.

14. The method of claim 1 , comprising calibrating the stimulation simulator using experimental data derived from a specimen representing rock from the treatment region to determine a reaction rate function for reaction between the treatment fluid and a solid material in the cells and to populate a table of the fluid mobility versus the solid saturation.

15. The method of claim 1 , comprising running the stimulation simulator a plurality of times to obtain datapoints comprising pore volume to breakthrough as a function of treatment fluid injection rate.

16. The method of claim 1 , comprising running the stimulation simulator a plurality of times to obtain datapoints comprising pore volume to breakthrough as a function of treatment fluid injection rate to determine the treatment fluid injection rate corresponding to a minimum pore volume to breakthrough.

17. The method of claim 1 , wherein the treatment region comprises a near-wellbore region of a subterranean formation, and further comprising running the stimulation simulator to determine an optimum treatment fluid injection rate to treat the near wellbore region.

18. The method of claim 1 , wherein the treatment region comprises a sector of a subterranean formation, and wherein the optimized treatment fluid injection parameters comprise an optimum treatment fluid injection rate to treat the sector.

19. The method of claim 1 , wherein the treatment region comprises a field of a subterranean formation, and wherein the optimized treatment fluid injection parameters comprise an optimum treatment fluid injection rate to treat the field.

20. A method, comprising:

modeling ( 40 ) a stimulation treatment involving a chemical reaction between a treatment fluid and a porous medium in a subterranean formation using a computerized model, comprising:

gridding a treatment region of the subterranean formation into a plurality of cells;

modeling the cells in a wormhole initiation stage wherein the medium of the cells having a solid saturation above a respective critical solid saturation is comprised of matrix material behaving as a single permeability, single porosity system;

modeling the cells having a solid saturation equal to or less than the respective critical solid saturation in a wormhole growth stage wherein the cells comprise two different interconnected media comprised of the matrix material and a wormhole material having a fluid mobility as a function of solid saturation; and

performing a stimulation treatment in a wellbore based on the modeling.