Modelling geological faults
To model hydrocarbon reservoirs more accurately, models of a geological volume using a three dimensional grid of cells populated with data from seismic studies are augmented by providing fine-grained stochastic models of the distribution of relays in fault zones and modelling the flow of hydrocarbons taking account the modelled relays.
1. A method of simulating fluid flow in a geological volume comprising a geological fault, the method comprising:
generating a three dimensional grid of cells representative of the structure of the geological volume, the cells comprising base cells;
classifying a plurality of base cells as fault cells corresponding to the fault in the geological volume;
classifying one or more sets of fault cells as relay cell sets corresponding to assumed locations of relays within the fault in the geological volume;
generating a three dimensional relay grid of relay sub-cells in the place of each relay cell set, the relay sub-cells of each relay grid being representative of a relay structure within the geological fault and having a finer grid-cell dimension than that of the cells of the relay cell set;
determining transmissibilities between the relay sub-cells;
determining transmissibilities between the relay cell sets based on the determined transmissibilities between the relay sub-cells, and
determining transmissibilities between the base cells based on the determined transmissibilities between the relay sub-cells, wherein the classification of the one or more relay cell sets provides a spatially even distribution of relay cell sets along the direction of the length of the fault.
2. The method according to claim 1 , wherein the classification of the one or more relay cell sets corresponds to a random spatial distribution of relay cell sets along the direction of the length of the fault.
3. The method according to claim 1 , wherein a poisson point process is employed in order to classify the one or more relay cell sets.
4. The method according to claim 1 , wherein the number of classified relay cell sets is determined based on the length of the fault represented by the fault cells.
5. The method according to claim 1 , wherein the fluid flow is simulated over multiple iterations, each new iteration comprising reclassifying one or more sets of cells as relay cell sets.
6. The method according to claim 1 , wherein one or more relay shape parameters are assigned to the relay cell sets and the calculation of the transmissibilities between the relay sub-cells takes into account the one or more relay shape parameters.
7. The method according to claim 6 , wherein the one or more relay shape parameters include width, length, and vertical extent of the relay.
8. The method according to claim 1 , wherein the classifying of the relay cell sets, generation of the relay grid of relay sub-cells and calculation of transmissibilities between relay sub-cells is repeated for each fault identified in the geological volume.
9. The method according to claim 1 , wherein the one or more sets of cells classified as relay cell sets are base cells.
10. The method according to claim 1 , wherein:
the vertical faces of the fault cells are vertically offset from the fault; and
the fault cells comprise cutted fault cells and non-cutted fault cells, the cutted fault cells having a reduced height as compared to the height of the non-cutted fault cells.
11. A processing device arranged to implement the method as claimed in claim 1 .
12. A computer readable medium comprising computer program instructions that when executed on a processor perform the method according to claim 1 .
13. A method of simulating fluid flow in a geological volume comprising a geological fault, the method comprising:
generating a three dimensional grid of cells representative of the structure of the geological volume, the cells comprising base cells;
classifying a plurality of base cells as fault cells corresponding to the fault in the geological volume;
classifying one or more sets of fault cells as relay cell sets corresponding to assumed locations of relays within the fault in the geological volume;
generating a three dimensional relay grid of relay sub-cells in the place of each relay cell set, the relay sub-cells of each relay grid being representative of a relay structure within the geological fault and having a finer grid-cell dimension than that of the cells of the relay cell set;
determining transmissibilities between the relay sub-cells;
determining transmissibilities between the relay cell sets based on the determined transmissibilities between the relay sub-cells, and
determining transmissibilities between the base cells based on the determined transmissibilities between the relay sub, wherein the number of classified relay cell sets is determined based on a random distribution.
14. A method of simulating fluid flow in a geological volume comprising a geological fault, the method comprising:
generating a three dimensional grid of cells representative of the structure of the geological volume, the cells comprising base cells;
classifying a plurality of base cells as fault cells corresponding to the fault in the geological volume;
classifying one or more sets of the fault cells as relay cell sets corresponding to assumed locations of relays within the fault in the geological volume;
generating a three dimensional relay grid of relay sub-cells in the place of each relay cell set, the relay sub-cells of each relay grid being representative of a relay structure within the geological fault and having a finer grid-cell dimension than that of the cells of the relay cell set;
determining transmissibilities between the relay sub-cells;
determining transmissibilities between the relay cell sets based on the determined transmissibilities between the relay sub-cells, and
determining transmissibilities between the base cells based on the determined transmissibilities between the relay sub-cells; and wherein:
the vertical faces of the fault cells are vertically offset from the fault; and
the fault cells comprise cutted fault cells and non-cutted fault cells, the cutted fault cells having a reduced height as compared to the height of the non-cutted fault cells, the method further comprising:
detecting hybrid-interface fault cells interfacing with both the fault and the cutted cells;
generating pseudo fault cells representing the fault cells in which the hybrid-interface fault cells are represented by pseudo fault cells having a height limited to the height of the portion of the hybrid-interface fault cells interfacing with the fault; and
the one or more sets of cells classified as relay cell sets are pseudo fault cells.