Iteratively Coupled Flow and Geomechanics in Fractured Poroelastic Reservoirs: A Phase Field Fracture Model

<div>Sneddon’s pressurized fracture. (a) Internally pressurized plane-strain fracture is centered in the 2D domain. (b) Phase field contour surface in a 2D mesh (<span class="nowrap"><svg height="8.69875pt" id="M150" style="vertical-align:-0.3499298pt" version="1.1" viewbox="-0.0498162 -8.34882 59.9977 8.69875" width="59.9977pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M456 178V225H360V632H320C217 496 115 347 20 206V178H280V106C280 40 276 34 189 27V0H445V27C364 34 360 39 360 106V178H456ZM280 225H82C149 335 214 431 278 520H280V225Z"></path></g><g transform="matrix(.013,0,0,-0.013,6.24,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z"></path></g><g transform="matrix(.013,0,0,-0.013,12.48,0)"><path d="M797 0V26C739 32 732 36 732 103V296C732 394 682 449 605 449C576 449 550 437 529 423C504 407 475 389 446 366C425 418 382 449 334 449C303 449 279 437 253 421C222 403 201 385 180 371V452C135 432 85 419 41 411V388C99 379 102 374 102 310V103C102 38 93 32 27 26V0H238V26C189 32 180 38 180 103V338C210 363 250 390 289 390C351 390 377 348 377 275V103C377 37 368 32 306 26V0H520V26C465 32 456 38 456 101V296C456 314 455 326 453 338C491 369 529 390 565 390C628 390 653 345 653 274V107C653 36 642 32 583 26V0H797Z"></path></g><g transform="matrix(.013,0,0,-0.013,26.073,0)"><path d="M528 54L331 254L528 455L492 493L294 291L96 493L60 455L257 254L60 54L96 16L294 217L492 16L528 54Z"></path></g><g transform="matrix(.013,0,0,-0.013,36.609,0)"><path d="M456 178V225H360V632H320C217 496 115 347 20 206V178H280V106C280 40 276 34 189 27V0H445V27C364 34 360 39 360 106V178H456ZM280 225H82C149 335 214 431 278 520H280V225Z"></path></g><g transform="matrix(.013,0,0,-0.013,42.849,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z"></path></g><g transform="matrix(.013,0,0,-0.013,49.089,0)"><path d="M797 0V26C739 32 732 36 732 103V296C732 394 682 449 605 449C576 449 550 437 529 423C504 407 475 389 446 366C425 418 382 449 334 449C303 449 279 437 253 421C222 403 201 385 180 371V452C135 432 85 419 41 411V388C99 379 102 374 102 310V103C102 38 93 32 27 26V0H238V26C189 32 180 38 180 103V338C210 363 250 390 289 390C351 390 377 348 377 275V103C377 37 368 32 306 26V0H520V26C465 32 456 38 456 101V296C456 314 455 326 453 338C491 369 529 390 565 390C628 390 653 345 653 274V107C653 36 642 32 583 26V0H797Z"></path></g></svg>).</span> <svg height="12.3916pt" id="M151" style="vertical-align:-3.42948pt" version="1.1" viewbox="-0.0498162 -8.96212 28.6589 12.3916" width="28.6589pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z"></path></g><g transform="matrix(.013,0,0,-0.013,10.978,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z"></path></g><g transform="matrix(.013,0,0,-0.013,22.241,0)"><path d="M384 0V27C293 34 287 42 287 114V635C232 613 172 594 109 583V559L157 557C201 555 205 550 205 499V114C205 42 199 34 109 27V0H384Z"></path></g></svg> represents the preexisting fracture and <svg height="12.3916pt" id="M152" style="vertical-align:-3.42948pt" version="1.1" viewbox="-0.0498162 -8.96212 28.6589 12.3916" width="28.6589pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M542 242C542 369 469 429 379 468C370 472 372 496 376 513L416 681L400 691L344 648L219 51C167 69 113 128 113 217C113 324 166 387 266 428L252 459C124 418 23 336 23 202C23 95 90 13 205 -8L182 -132C165 -222 191 -254 209 -261L215 -258L270 -7C415 -4 542 98 542 242ZM469 232C469 126 403 33 279 39L355 409C399 394 469 331 469 232Z"></path></g><g transform="matrix(.013,0,0,-0.013,10.978,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z"></path></g><g transform="matrix(.013,0,0,-0.013,22.241,0)"><path d="M241 635C89 635 35 457 35 312C35 153 89 -12 240 -12C390 -12 443 166 443 312C443 466 390 635 241 635ZM238 602C329 602 354 454 354 312C354 172 330 22 240 22C152 22 124 173 124 313S148 602 238 602Z"></path></g></svg> represents the elastic medium. Zero displacement Dirichlet boundary conditions are applied for all the domain boundaries. However, the numerical simulation is insensitive to other types of displacement boundary conditions since the computational domain is large relative to the fracture length.</div>

Geofluids

fig7

Figure 7

Figure 7: Iteratively Coupled Flow and Geomechanics in Fractured Poroelastic Reservoirs: A Phase Field Fracture Model 