Method and device for determining damage to a geological model

A device for determining damage to a geological model, said geological model comprising a plurality of nodes and a plurality of links, each link from the plurality of links linking together two nodes from the plurality of nodes, the device comprising a processor, wherein the processor performs the following:

• for each link from the plurality of links:
• determining a force applied to said link,
• if said force is greater than a force threshold value:
• for each node of the two nodes from the plurality of nodes linked together by said link, creating a contact sphere associated with said node, and
• eliminating said link,
• determining a plurality of contacts between spheres comprising at least one contact between spheres, each contact between spheres associating two contact spheres from the plurality of contact spheres, wherein
• a first plurality of nodes from the plurality of node of the geological model is associated with a first geological phase modeled by the geological model,
• a second plurality of nodes from the plurality of nodes of the geological model is associated with a second geological phase modeled by the geological model,
• a plurality of first geological phase links from the plurality of links of the geological model consists of links from the plurality of links linking together two nodes from the first plurality of nodes, with each link from the plurality of first geological phase links is associated a first geological phase rigidity value kϕi and a first geological phase force threshold value fcϕ1 determined on the basis of an elastic modulus value of the first geological phase Eϕ1 modeled by the geological model and a tensile strength threshold value of the first geological phase σϕ1 modeled by the geological model,
• a plurality of second geological phase links from the plurality of links of the geological model consists of links from the plurality of links linking together two nodes from the second plurality of nodes, with each link from the plurality of second geological phase links is associated a second geological phase rigidity value kϕ2 and a second geological phase force threshold value fcϕ2 determined on the basis of an elastic modulus values Eϕ2 and of the second geological phase and a tensile strength threshold value σϕ2 of the second geological phase modeled by the geological model,
• a plurality of interface links from the plurality of links of the geological model consists of links from the plurality of links linking together a node of the first plurality of nodes and a node from the second plurality of nodes, with each link from the plurality of interface links is associated an interface rigidity value kinterface and an interface force threshold value fcinterface, and
• wherein the first geological phase rigidity value kϕ1 associated with one of the first geological phase link is different to the second geological phase rigidity value kϕ2 associated with one of the second geological phase link and/or the first geological phase force threshold value fcϕ1 associated with one of the first geological phase link is different to the second geological phase force threshold value fcϕ2 associated with one of the second geological phase link.

A computer program product comprising instructions for the implementation, when this program is executed by a processor, for determining damage to a geological model, said geological model comprising a plurality of nodes and a plurality of links, each link from the plurality of links linking together two nodes from the plurality of nodes, comprising the following steps:

• for each link from the plurality of links:
• determining a force applied to said link,
• if said force is greater than a force threshold value:
• for each node of the two nodes from the plurality of nodes linked together by said link, creating a contact sphere associated with said node, and
• eliminating said link,
• determining a plurality of contacts between spheres comprising at least one contact between spheres, each contact between spheres associating two contact spheres from the plurality of contact spheres, wherein
• a first plurality of nodes from the plurality of node of the geological model is associated with a first geological phase modeled by the geological model,
• a second plurality of nodes from the plurality of nodes of the geological model is associated with a second geological phase modeled by the geological model,
• a plurality of first geological phase links from the plurality of links of the geological model consists of links from the plurality of links linking together two nodes from the first plurality of nodes, with each link from the plurality of first geological phase links is associated a first geological phase rigidity value kϕ1 and a first geological phase force threshold value fcϕ1 determined on the basis of an elastic modulus value of the first geological phase Eϕ1 modeled by the geological model and a tensile strength threshold value of the first geological phase σϕ1 modeled by the geological model,
• a plurality of second geological phase links from the plurality of links of the geological model consists of links from the plurality of links linking together two nodes from the second plurality of nodes, with each link from the plurality of second geological phase links is associated a second geological phase rigidity value kϕ2 and a second geological phase force threshold value fcϕ2 determined on the basis of elastic modulus values Eϕ2 of the second geological phase and tensile strength threshold values σϕ2 of the second geological phase modeled by the geological model,
• a plurality of interface links from the plurality of links of the geological model consists of links from the plurality of links linking together a node of the first plurality of nodes and a node from the second plurality of nodes, with each link from the plurality of interface links is associated an interface rigidity value kinterface and an interface force threshold value fcinterface, and
• wherein the first geological phase rigidity value kϕ1 associated with one of the first geological phase link is different to the second geological phase rigidity value kϕ2 associated with one of the second geological phase link and/or the first geological phase force threshold value fcϕ1 associated with one of the first geological phase link is different to the second geological phase force threshold value fcϕ2 associated with one of the second geological phase link.