Upscaled Lattice Boltzmann Method for Simulations of Flows in Heterogeneous Porous Media

<div>Comparisons of <svg height="10.2124pt" id="M449" style="vertical-align:-3.42943pt" version="1.1" viewbox="-0.0498162 -6.78297 7.83752 10.2124" width="7.83752pt" 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="M570 304C570 398 525 448 414 448C385 448 343 445 312 434L329 511L321 518C297 504 262 482 244 460L233 411C195 397 159 381 128 358L135 332C160 347 189 360 224 373L111 -147C97 -210 84 -218 17 -231L13 -257L254 -247L259 -218L233 -216C183 -212 177 -202 189 -142L218 -1C238 -10 266 -12 283 -12C351 3 429 48 483 105C543 168 570 242 570 304ZM482 289C482 161 380 33 304 33C278 33 248 51 233 69L303 396C326 400 352 403 369 403C428 403 482 380 482 289Z" id="g113-113"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="593" vert-adv-y="593"></glyph.data></g></svg>, <svg height="6.1673pt" id="M450" style="vertical-align:-0.2063904pt" version="1.1" viewbox="-0.0498162 -5.96091 7.11985 6.1673" width="7.11985pt" 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="M515 96L502 119C471 88 431 62 423 62C416 62 411 70 417 101C440 223 469 341 497 448H486L412 422L380 277C330 188 210 57 152 57C137 57 126 69 139 124L195 366C210 431 205 448 182 448C155 448 89 413 23 350L36 326C73 354 103 376 112 376C118 376 118 365 113 340L61 118C54 90 52 68 52 51C52 0 75 -12 98 -12S151 -3 181 17C242 58 305 119 362 193H364L345 104C323 3 339 -12 359 -12C390 -12 464 35 515 96Z" id="g113-118"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="538" vert-adv-y="538"></glyph.data></g></svg>, and <svg height="6.02377pt" id="M451" style="vertical-align:-0.2063999pt" version="1.1" viewbox="-0.0498162 -5.81737 5.98461 6.02377" width="5.98461pt" 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="M219 86C216 168 211 250 206 337C201 410 189 448 163 448C131 448 79 396 43 344L60 322C91 359 110 375 118 375S132 358 136 298C141 238 152 81 155 -12H182C242 62 331 177 390 258C435 321 451 360 451 391C450 424 432 448 408 448C390 448 372 435 366 419C362 410 362 401 366 394C373 383 376 367 376 350C376 283 262 138 221 86H219Z" id="g185-40"></path><glyph.data ascent="989" descent="-303" horiz-adv-x="447" vert-adv-y="447"></glyph.data></g></svg> between two fine-grid simulations using <svg height="14.7318pt" id="M452" style="vertical-align:-3.2957pt" version="1.1" viewbox="-0.0498162 -11.4361 54.5286 14.7318" width="54.5286pt" 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="M213 67C239 136 257 225 255 300C253 353 242 447 164 447C112 447 54 389 41 319L67 309C90 362 117 377 134 377C184 377 189 277 190 239C193 153 164 43 146 -12H176C278 73 450 253 450 390C450 439 410 463 396 463C384 463 380 459 382 449C385 440 389 415 389 402C389 284 312 146 216 67H213Z" id="g185-47"></path><glyph.data ascent="989" descent="-303" horiz-adv-x="454" vert-adv-y="454"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,5.902,3.132)"><path d="M384 109C349 76 310 59 268 59C197 59 123 113 120 245C240 248 367 259 384 262C403 266 407 275 407 296C407 376 338 451 254 451H253C201 451 148 425 105 378C65 334 39 272 39 203C39 88 112 -12 235 -12C266 -12 339 5 402 84L384 109ZM230 411C285 411 319 365 319 313C319 298 313 293 294 293C236 291 180 290 125 290C141 370 185 411 230 411Z" id="g58-99"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="434" vert-adv-y="434"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,9.851,3.132)"><path d="M100 439H52L25 408L28 396H100V106C100 41 92 35 31 30V0H250V29C188 37 181 41 181 106V396H372V107C372 40 365 37 309 29V0H537V29C461 37 453 40 453 106V396H564C576 403 579 428 574 439H453V469C452 558 462 604 474 623S500 656 527 656C562 656 589 634 608 615C619 604 628 607 639 615C647 622 658 633 661 641C666 653 664 664 656 674C637 696 609 708 572 710C528 703 483 682 439 638C418 661 380 685 320 688C269 683 230 666 190 633C139 593 122 550 113 525S100 476 100 451V439ZM372 439H181V461C181 530 194 571 211 598C226 621 248 643 287 643C343 643 379 607 400 578C391 561 384 546 381 529C376 510 372 486 372 457V439Z" id="g58-26"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="572" vert-adv-y="572"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,20.061,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,31.324,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" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,37.564,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" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,43.804,-5.741)"><path d="M556 236V289H56V236H556Z" id="g54-33"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,49.364,-5.741)"><path d="M158 548H390L417 615L410 623H122L83 318C105 326 143 337 185 337C296 337 350 275 350 188C350 116 308 42 225 42C164 42 122 74 100 93C90 101 82 99 72 92C60 82 51 68 50 59C48 46 52 38 66 24C82 9 125 -12 172 -12C225 -11 292 15 346 59C408 108 437 166 437 226C437 309 371 397 242 397C214 397 170 382 133 369L158 548Z" id="g50-54"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg> and 0 m<sup>2</sup> s<sup>−1</sup>, respectively, <svg height="15.6668pt" id="M453" style="vertical-align:-3.429399pt" version="1.1" viewbox="-0.0498162 -12.2374 102.436 15.6668" width="102.436pt" xmlns="http://www.w3.org/2000/svg" xmlns:xlink="http://www.w3.org/1999/xlink"><g transform="matrix(.013,0,0,-0.013,8.338,-2.897)"><path d="M-21 591C-45 610 -87 648 -131 700L-155 684L-108 616H-397V566H-108L-155 498L-131 482C-87 534 -45 572 -21 591Z" id="g117-14"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="0" vert-adv-y="0"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,0,0)"><path d="M673 297H417L411 270C517 261 522 258 506 183L487 92C476 38 428 19 360 19C207 19 123 132 123 287C123 472 243 631 441 631C557 631 617 583 617 469L647 472C650 545 655 604 658 632C625 643 554 666 467 666C201 666 23 502 23 278C23 90 156 -16 339 -16C410 -16 501 6 569 24C566 43 567 70 573 101L589 183C604 259 607 262 667 271L673 297Z" id="g113-72"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="684" vert-adv-y="684"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,12.557,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,23.82,0)"><path d="M300 -147C201 -63 143 98 143 270S200 602 300 686L282 710C136 610 70 450 70 271V270C70 89 136 -72 282 -170L300 -147Z" id="g113-41"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,28.318,0)"><path d="M319 325C317 349 306 409 297 431C277 440 250 449 209 449C117 449 57 389 57 319C57 243 122 209 182 182C232 159 261 135 261 91C261 48 227 21 190 21C130 21 85 79 68 145L41 140C41 104 51 36 58 22C75 7 121 -12 172 -12C252 -12 337 35 337 126C337 195 286 231 210 262C166 281 126 304 126 348C126 388 152 417 191 417C240 417 274 378 294 318L319 325Z" id="g121-113"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="367" vert-adv-y="367"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,33.089,0)"><path d="M135 536C164 536 186 560 186 587C186 617 164 639 136 639C109 639 85 617 85 587C85 560 109 536 135 536ZM252 0V26C188 32 181 38 181 106V451C138 433 90 420 39 412V388C99 379 102 374 102 312V106C102 38 95 32 32 26V0H252Z" id="g121-103"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="268" vert-adv-y="268"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,36.573,0)"><path d="M524 0V26C466 32 460 36 460 104V297C460 393 411 449 331 449C302 449 276 437 248 419C223 402 201 387 181 372V451C137 432 90 420 42 411V388C96 378 102 374 102 310V104C102 38 97 33 29 26V0H246V26C187 32 181 36 181 104V339C211 365 250 390 290 390C357 390 381 345 381 276V109C381 40 374 32 315 26V0H524Z" id="g121-108"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="547" vert-adv-y="547"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,45.921,0)"><path d="M574 449L545 460C526 432 516 430 487 430C404 430 311 435 226 435C104 435 56 379 25 341L43 318C81 354 121 372 181 372C161 246 87 53 23 3L30 -12C48 -12 88 -4 113 11C157 75 207 248 232 371L386 367L326 109C321 86 318 66 318 50C318 4 339 -12 366 -12C410 -12 461 21 505 69L492 96C467 79 434 60 418 60C406 60 400 78 411 147C422 217 439 300 457 366C487 366 524 367 536 370C547 385 558 408 574 449Z" id="g113-238"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="591" vert-adv-y="591"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,53.474,0)"><path d="M536 404C536 423 520 448 491 448C445 448 398 404 308 283L286 338C255 416 242 448 217 448C182 448 138 402 92 341L111 321C149 368 169 378 178 378C188 378 198 363 214 320L254 212C185 117 138 65 107 65C98 65 85 69 82 75C79 82 73 86 65 86C44 86 23 60 23 39C23 7 44 -12 71 -12C119 -12 168 33 265 177L306 61C321 17 347 -12 373 -12C413 -12 465 33 507 96L491 119C459 84 432 60 413 60C395 60 378 92 358 148L321 250C341 279 369 310 389 332C417 363 439 382 456 382C466 382 475 376 481 368C486 361 492 358 496 358C513 358 536 381 536 404Z" id="g113-121"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="559" vert-adv-y="559"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,60.742,0)"><path d="M95 130C70 130 46 113 46 88C46 72 54 64 59 64C93 55 121 33 121 -3C121 -41 93 -68 44 -88L55 -117C117 -98 186 -56 186 22C186 91 131 130 95 130Z" id="g113-45"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="228" vert-adv-y="228"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,65.885,0)"><path d="M319 325C317 349 306 409 297 431C277 440 250 449 209 449C117 449 57 389 57 319C57 243 122 209 182 182C232 159 261 135 261 91C261 48 227 21 190 21C130 21 85 79 68 145L41 140C41 104 51 36 58 22C75 7 121 -12 172 -12C252 -12 337 35 337 126C337 195 286 231 210 262C166 281 126 304 126 348C126 388 152 417 191 417C240 417 274 378 294 318L319 325Z" id="g121-113"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="367" vert-adv-y="367"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,70.656,0)"><path d="M135 536C164 536 186 560 186 587C186 617 164 639 136 639C109 639 85 617 85 587C85 560 109 536 135 536ZM252 0V26C188 32 181 38 181 106V451C138 433 90 420 39 412V388C99 379 102 374 102 312V106C102 38 95 32 32 26V0H252Z" id="g121-103"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="268" vert-adv-y="268"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,74.14,0)"><path d="M524 0V26C466 32 460 36 460 104V297C460 393 411 449 331 449C302 449 276 437 248 419C223 402 201 387 181 372V451C137 432 90 420 42 411V388C96 378 102 374 102 310V104C102 38 97 33 29 26V0H246V26C187 32 181 36 181 104V339C211 365 250 390 290 390C357 390 381 345 381 276V109C381 40 374 32 315 26V0H524Z" id="g121-108"></path><glyph.data ascent="989" descent="-360" horiz-adv-x="547" vert-adv-y="547"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,83.488,0)"><path d="M574 449L545 460C526 432 516 430 487 430C404 430 311 435 226 435C104 435 56 379 25 341L43 318C81 354 121 372 181 372C161 246 87 53 23 3L30 -12C48 -12 88 -4 113 11C157 75 207 248 232 371L386 367L326 109C321 86 318 66 318 50C318 4 339 -12 366 -12C410 -12 461 21 505 69L492 96C467 79 434 60 418 60C406 60 400 78 411 147C422 217 439 300 457 366C487 366 524 367 536 370C547 385 558 408 574 449Z" id="g113-238"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="591" vert-adv-y="591"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,90.128,0)"><path d="M556 393C556 426 537 448 514 448C496 448 478 435 471 420C468 414 462 401 466 394C473 382 480 368 480 346C480 268 392 143 338 67H336C329 163 319 253 309 330S286 448 254 448C215 448 160 383 127 337L143 311C167 344 200 373 208 373S222 365 229 320C246 214 264 66 268 -20C219 -83 131 -171 17 -239L25 -261L137 -235C248 -110 273 -79 335 8C384 77 481 215 520 287C544 332 556 363 556 393Z" id="g113-122"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="579" vert-adv-y="579"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,97.657,0)"><path d="M275 270C275 450 212 609 64 710L45 686C145 604 203 442 203 270S147 -63 45 -147L64 -170C213 -68 275 89 275 270Z" id="g113-42"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="346" vert-adv-y="346"></glyph.data></g></svg> m s<sup>−2</sup>, and <svg height="14.7318pt" id="M454" style="vertical-align:-3.2957pt" version="1.1" viewbox="-0.0498162 -11.4361 64.6857 14.7318" width="64.6857pt" 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="M484 406C484 429 465 448 438 448C400 448 352 405 302 356L156 213L191 373C195 391 197 404 197 416C197 436 191 448 178 448C151 448 83 406 23 350L37 328C67 350 101 369 109 369C114 369 115 362 111 344L29 -5L34 -12C63 -5 86 -1 111 3C124 73 137 138 147 161C165 180 189 202 206 217C215 194 249 111 275 57C304 -1 319 -12 342 -12C368 -12 405 -5 471 92L454 115C427 80 395 55 384 55C374 55 364 58 341 106L262 274C282 296 367 378 414 378C427 378 440 370 444 367C449 363 454 362 462 368C476 379 484 395 484 406Z" id="g113-232"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="507" vert-adv-y="507"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,6.266,3.132)"><path d="M396 114C349 71 317 59 274 59C216 59 123 104 123 241C123 342 180 400 246 400C283 400 317 386 348 359C355 352 361 348 366 348C378 348 401 371 402 394C402 405 398 413 385 424C369 437 337 451 294 451H293C254 451 193 434 141 394C73 342 39 277 39 198C39 90 114 -12 243 -12C303 -12 371 31 415 90L396 114Z" id="g58-97"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="432" vert-adv-y="432"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,10.197,3.132)"><path d="M262 451C169 451 38 375 38 210C38 98 121 -12 261 -12C360 -12 481 64 481 227C481 351 386 451 262 451ZM249 414C337 414 385 320 385 204C385 70 333 25 272 25C189 25 134 117 134 241C134 355 189 414 249 414Z" id="g58-109"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="519" vert-adv-y="519"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,14.81,3.132)"><path d="M532 0V30C473 38 467 38 467 107V298C467 395 416 451 336 451C307 451 279 439 251 420C226 404 205 389 185 374V453C142 434 94 422 43 413V386C98 378 104 373 104 309V107C104 41 100 38 31 30V0H251V30C191 35 185 40 185 107V337C216 364 255 389 295 389C360 389 386 347 386 276V110C386 44 380 35 321 30V0H532Z" id="g58-108"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="557" vert-adv-y="557"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,19.797,3.132)"><path d="M327 325C324 350 312 408 304 433C283 442 256 451 213 451C120 451 59 391 59 320C59 243 125 210 188 181C237 158 266 136 266 93C266 52 231 24 195 24C134 24 93 81 74 148L42 143C42 106 51 37 59 23C78 8 125 -12 177 -12C257 -12 345 35 345 127C345 198 292 234 216 264C171 283 131 305 131 348C131 387 156 416 195 416C243 416 277 378 297 317L327 325Z" id="g58-113"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="376" vert-adv-y="376"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,23.173,3.132)"><path d="M303 37L293 66C280 59 256 49 231 49S173 63 173 143V396H281C294 404 298 428 287 439H173V575L158 577L92 509V439H46L17 408L22 396H92V107C92 27 128 -12 192 -12C203 -12 218 -9 237 2L303 37Z" id="g58-114"></path><glyph.data ascent="1007" descent="-369" horiz-adv-x="312" vert-adv-y="312"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,30.217,0)"><path d="M535 323V373H52V323H535ZM535 138V188H52V138H535Z" id="g117-34"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="587" vert-adv-y="587"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,41.48,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" id="g113-50"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.013,0,0,-0.013,47.72,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" id="g113-49"></path><glyph.data ascent="3473" descent="-2876" horiz-adv-x="480" vert-adv-y="480"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,53.96,-5.741)"><path d="M556 236V289H56V236H556Z" id="g54-33"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="611" vert-adv-y="611"></glyph.data></g><g transform="matrix(.0091,0,0,-0.0091,59.52,-5.741)"><path d="M453 623H65C60 580 56 530 46 472H80C103 538 110 548 175 548H389C310 380 195 168 90 0L98 -12L175 -2C273 206 366 408 462 610L453 623Z" id="g50-56"></path><glyph.data ascent="3443" descent="-2856" horiz-adv-x="487" vert-adv-y="487"></glyph.data></g></svg> m<sup>2</sup>.</div>

Geofluids

fig7

Figure 7

Figure 7: Upscaled Lattice Boltzmann Method for Simulations of Flows in Heterogeneous Porous Media 