شبیه سازی یک روش صدای الکترومغناطیسی برای اکتشاف نفت

View Categories

شبیه سازی یک روش صدای الکترومغناطیسی برای اکتشاف نفت

10 min read

PDF

شبیه سازی یک روش صدای الکترومغناطیسی برای اکتشاف نفت
معرفی
روش الکترومغناطیسی منبع کنترل شده دریایی (CSEM) برای اکتشاف نفت به عنوان یک تکنیک امیدوارکننده در سال‌های اخیر ظهور کرده است. این مثال یکی از انواع آن را نشان می دهد. این آنتن از یک آنتن دوقطبی افقی متحرک 1  هرتزی استفاده می کند که در ارتفاع 150  متری از کف دریا بکسل می شود. مجموعه ای از گیرنده های کف دریا میدان الکتریکی را در فواصل مختلف دور از آنتن اندازه گیری می کنند. هنگام اندازه گیری در فاصله به اندازه کافی بزرگ، مقداری از انرژی ارسال شده توسط مخزن مقاومتی منعکس/هدایت می شود و سیگنال دریافتی بالاتری نسبت به زمانی که هیچ مخزنی وجود نداشته باشد، می شود.
تعریف مدل
حوزه محاسباتی کره ای به شعاع 5  کیلومتر است. ناحیه بالای این کره نشان دهنده هوا است. در فرکانس پایینی مانند 1  هرتز، مقداری تثبیت عددی در این حوزه مورد نیاز است، بنابراین رسانایی مصنوعی 0.001  S/m مشخص شده است. یک حوزه آب اقیانوسی به عمق 100  متر با رسانایی 3  S/m و گذردهی نسبی 80 در بالای صفحه میانی کره مشخص شده است. در زیر صفحه میانی، رسانایی 1.5  S/m و گذردهی نسبی 30 برای سنگ مشخص شده است. در عمق متوسط ​​250  متری صخره جاسازی شده، بلوکی شکل به عمق 100  متر و مخزن هیدروکربنی به عرض 4  کیلومتر در 1 کیلومتر. رسانایی لایه هیدروکربنی 0.01  S/m و میزان گذردهی آن 4 است . فرستنده به عنوان یک  قطعه جریان AC با دامنه کوتاه 10 کیلو آمپری در 150  متر بالاتر از سطح میانی مدل سازی شده است. در مرزهای کروی خارجی، شرایط مرزی نوع پراکندگی امواج کروی خروجی را جذب می کند. معادله زیر برای بردار میدان الکتریکی E در حوزه محاسباتی حل می‌شود:
که در آن μ r نشان دهنده نفوذپذیری نسبی، j واحد خیالی، σ رسانایی، ω فرکانس زاویه ای، ε r گذردهی نسبی، و ε0 گذردهی فضای آزاد است .
نتایج و بحث
شکل اول زیر یک نمودار برشی از بزرگی میدان الکتریکی را در مقیاس دسی بل نشان می دهد که در آن اثر هدایت کننده لایه هیدروکربنی به وضوح قابل مشاهده است.
شکل 1: بزرگی میدان الکتریکی در یک برش حاوی آنتن (نقطه قرمز) رسم شده است. اثر هدایت لایه هیدروکربنی به وضوح قابل مشاهده است.
این اثر به صورت کمی در نمودار بعدی نشان داده شده است، جایی که مقایسه ای از بزرگی میدان الکتریکی در کف دریا به عنوان تابعی از فاصله برای شبیه سازی با و بدون لایه هیدروکربنی ترسیم شده است. حداکثر درست در زیر آنتن به دست می آید که در آن تفاوت کمی بین دو مورد وجود دارد. هنگامی که از آنتن دور می شوید، هنگامی که لایه هیدروکربنی وجود دارد، قدرت سیگنال به طور قابل توجهی بالاتر است.
شکل 2: بزرگی میدان الکتریکی در کف دریا به عنوان تابعی از فاصله. خط چین نشان دهنده مورد بدون لایه هیدروکربنی است.
نکاتی در مورد پیاده سازی COMSOL
این مثال از فرمول امواج الکترومغناطیسی سه بعدی موجود در ماژول RF استفاده می کند، اگرچه طول موج بسیار بزرگتر از دامنه محاسباتی است. شما همچنین می توانید این مثال را با استفاده از یک فرمول شبه استاتیک در ماژول AC/DC ایجاد کنید.
مسیر کتابخانه برنامه: RF_Module/Scattering_and_RCS/marine_csem
دستورالعمل های مدل سازی
از منوی File ، New را انتخاب کنید .
جدید
در پنجره جدید ، روی  Model  Wizard کلیک کنید .
مدل جادوگر
1
در پنجره Model  Wizard ، روی  3D کلیک کنید .
2
در درخت انتخاب  فیزیک ، فرکانس رادیویی   امواج الکترومغناطیسی،  دامنه فرکانس  (emw) را انتخاب کنید .
3
روی افزودن کلیک کنید .
4
 روی مطالعه کلیک کنید .
5
در درخت انتخاب  مطالعه ، General  Studies>Frequency  Domain را انتخاب کنید .
6
 روی Done کلیک کنید .
هندسه 1
کره 1 (sph1)
1
در نوار ابزار Geometry ، روی  Sphere کلیک کنید .
2
در پنجره تنظیمات برای Sphere ، بخش Size را پیدا کنید .
3
در قسمت متن Radius ، 5e3 را تایپ کنید .
سیلندر 1 (cyl1)
1
در نوار ابزار Geometry ، روی  Cylinder کلیک کنید .
2
در پنجره تنظیمات سیلندر ، بخش اندازه  و  شکل را پیدا کنید .
3
در قسمت متن Radius ، 6e3 را تایپ کنید .
4
در قسمت متن ارتفاع ، 1e3 را تایپ کنید .
تقاطع 1 (int1)
1
در نوار ابزار Geometry ، روی  Booleans  and  Partitions کلیک کنید و Intersection را انتخاب کنید .
2
در پنجره Graphics کلیک کنید و سپس Ctrl+A را فشار دهید تا هر دو شی انتخاب شوند.
کره 2 (sph2)
1
در نوار ابزار Geometry ، روی  Sphere کلیک کنید .
2
در پنجره تنظیمات برای Sphere ، بخش Size را پیدا کنید .
3
در قسمت متن Radius ، 5e3 را تایپ کنید .
بلوک 1 (blk1)
1
در نوار ابزار Geometry ، روی  Block کلیک کنید .
2
در پنجره تنظیمات برای Block ، قسمت Size  and  Shape را پیدا کنید .
3
در قسمت متن Width ، 4e3 را تایپ کنید .
4
در قسمت متن عمق ، 1e3 را تایپ کنید .
5
در قسمت متن ارتفاع ، 100 را تایپ کنید .
6
قسمت Position را پیدا کنید . از لیست پایه ، مرکز را انتخاب کنید .
7
در قسمت متن z ، -250 را تایپ کنید .
بخش خط 1 (ls1)
1
در نوار ابزار Geometry ، روی  More  Primitives کلیک کنید و Line  Segment را انتخاب کنید .
2
در پنجره تنظیمات برای بخش خط  ، بخش Starting Point را پیدا کنید .
3
از لیست Specify ، Coordinates را انتخاب کنید .
4
قسمت Endpoint را پیدا کنید . از لیست Specify ، Coordinates را انتخاب کنید .
5
قسمت Starting  Point را پیدا کنید . در قسمت متن x ، -10 را تایپ کنید  و z را تا 150 تایپ کنید .
6
قسمت Endpoint را پیدا کنید . در قسمت متن x ، 10 را تایپ کنید  و z را تا 150 تایپ کنید .
بخش خط 2 (ls2)
1
در نوار ابزار Geometry ، روی  More  Primitives کلیک کنید و Line  Segment را انتخاب کنید .
2
در پنجره تنظیمات برای بخش خط  ، بخش Starting Point را پیدا کنید .
3
از لیست Specify ، Coordinates را انتخاب کنید .
4
قسمت Endpoint را پیدا کنید . از لیست Specify ، Coordinates را انتخاب کنید .
5
قسمت Starting  Point را پیدا کنید . در قسمت متن x ، -5e3 را تایپ کنید .
6
قسمت Endpoint را پیدا کنید . در قسمت متن x ، 5e3 را تایپ کنید .
7
 روی Build  All  Objects کلیک کنید .
8
 روی دکمه Wireframe  Rendering در نوار ابزار Graphics کلیک کنید .
مواد
آب دریا
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Materials راست کلیک کرده و Blank  Material را انتخاب کنید .
2
در پنجره تنظیمات برای مواد ، دریای آب را در قسمت متن برچسب تایپ کنید .
3
فقط دامنه 1 را انتخاب کنید.
4
قسمت محتوای مواد  را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
مجوز نسبی
epsilonr_iso ; epsilonrii = epsilonr_iso، epsilonrij = 0
80
1
پایه ای
نفوذپذیری نسبی
mur_eye ; murii = mur_iso، murij = 0
1
1
پایه ای
رسانایی الکتریکی
sigma_iso ; sigmaii = sigma_iso، sigmaij = 0
3
S/m
پایه ای
سنگ
1
روی Materials کلیک راست کرده و Blank  Material را انتخاب کنید .
2
در پنجره تنظیمات برای مواد ، Rock را در قسمت متن برچسب تایپ کنید .
3
فقط دامنه های 2 و 4 را انتخاب کنید.
4
قسمت محتوای مواد  را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
مجوز نسبی
epsilonr_iso ; epsilonrii = epsilonr_iso، epsilonrij = 0
30
1
پایه ای
نفوذپذیری نسبی
mur_eye ; murii = mur_iso، murij = 0
1
1
پایه ای
رسانایی الکتریکی
sigma_iso ; sigmaii = sigma_iso، sigmaij = 0
1.5
S/m
پایه ای
روغن
1
روی Materials کلیک راست کرده و Blank  Material را انتخاب کنید .
2
در پنجره تنظیمات مواد ، Oil را در قسمت نوشتار Label تایپ کنید .
3
فقط دامنه 4 را انتخاب کنید.
توجه داشته باشید که این انتخاب تا حدی تنظیمات محدوده هندسی را که برای سنگ تعریف کردید، لغو می کند. این به شما امکان می‌دهد بعداً با غیرفعال کردن Oil Node به راحتی مدل را برای مواردی که ذخیره هیدروکربنی وجود ندارد اجرا کنید.
4
قسمت محتوای مواد  را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
مجوز نسبی
epsilonr_iso ; epsilonrii = epsilonr_iso، epsilonrij = 0
4
1
پایه ای
نفوذپذیری نسبی
mur_eye ; murii = mur_iso، murij = 0
1
1
پایه ای
رسانایی الکتریکی
sigma_iso ; sigmaii = sigma_iso، sigmaij = 0
0.01
S/m
پایه ای
هوا
1
روی Materials کلیک راست کرده و Blank  Material را انتخاب کنید .
2
در پنجره تنظیمات برای مواد ، Air را در قسمت متن برچسب تایپ کنید .
3
فقط دامنه 3 را انتخاب کنید.
4
قسمت محتوای مواد  را پیدا کنید . در جدول تنظیمات زیر را وارد کنید:
 
ویژگی
متغیر
ارزش
واحد
گروه اموال
مجوز نسبی
epsilonr_iso ; epsilonrii = epsilonr_iso، epsilonrij = 0
1
1
پایه ای
نفوذپذیری نسبی
mur_eye ; murii = mur_iso، murij = 0
1
1
پایه ای
رسانایی الکتریکی
sigma_iso ; sigmaii = sigma_iso، sigmaij = 0
0.001
S/m
پایه ای
امواج الکترومغناطیسی، دامنه فرکانس (EMW)
شرایط مرزی پراکندگی 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی امواج الکترومغناطیسی ،  دامنه فرکانس (emw) کلیک راست کرده و Scattering Boundary Condition را انتخاب کنید .
2
در پنجره تنظیمات برای شرایط مرزی Scattering  ، بخش انتخاب مرز را پیدا کنید .
3
از لیست انتخاب ، همه  مرزها را انتخاب کنید .
4
قسمت Scattering  Boundary  Condition را پیدا کنید . از لیست نوع موج پراکنده  ، موج کروی را انتخاب کنید .
جریان لبه 1
1
در نوار ابزار Physics ، روی  Edges کلیک کنید و Edge  Current را انتخاب کنید .
اکنون می خواهید جریانی را روی لبه بسیار کوتاه نزدیک مرکز هندسه اعمال کنید. اگر چند بار زوم کنید این کار آسان تر خواهد شد.
2
 روی دکمه Zoom  In در نوار ابزار Graphics کلیک کنید .
3
فقط Edge 17 را انتخاب کنید.
4
در پنجره Model  Builder ، روی Edge  Current  1 کلیک کنید .
5
در پنجره تنظیمات Edge  Current ، قسمت Edge  Current را پیدا کنید .
6
در قسمت متن 0 ، 1e4 را تایپ کنید .
7
 روی دکمه Zoom  Extents در نوار ابزار Graphics کلیک کنید .
مش 1
1
در پنجره Model  Builder ، در قسمت Component   (comp1) روی Mesh  1 کلیک کنید .
2
در پنجره تنظیمات مش ، قسمت Physics-Controlled  Mesh را پیدا کنید .
3
در جدول، تیک Use for Waves Electromagnetic ,  Frequency Domain (emw) را پاک کنید .
4
از لیست اندازه عنصر  ، Fine را انتخاب کنید .
5
 روی ساخت  همه کلیک کنید .
مطالعه 1
مرحله 1: دامنه فرکانس
1
در پنجره Model  Builder ، در بخش مطالعه  1 ، روی مرحله  1:  دامنه فرکانس  کلیک کنید .
2
در پنجره تنظیمات دامنه فرکانس  ، بخش تنظیمات مطالعه را پیدا کنید .
3
در قسمت متن فرکانس ، 1[Hz] را تایپ کنید .
4
در نوار ابزار صفحه اصلی ،  روی محاسبه کلیک کنید .
نتایج
به بازتولید شکل 1 ادامه دهید .
چند برش
1
در پنجره Model  Builder ، گره Results>Electric  Field  (emw) را گسترش دهید .
2
روی Results>Electric  Field  (emw)>Multislice کلیک راست کرده و Delete را انتخاب کنید .
برش 1
1
در پنجره Model  Builder ، روی Electric  Field  (emw) کلیک راست کرده و Slice را انتخاب کنید .
2
در پنجره تنظیمات برای Slice ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، 20*log10(emw.normE) را تایپ کنید .
4
چک باکس Description را انتخاب کنید . در فیلد متن مرتبط، dB را وارد کنید Electric field norm .
5
قسمت Plane  Data را پیدا کنید . از لیست هواپیما ، ZX-planes را انتخاب کنید .
6
در قسمت متن Planes ، 1 را تایپ کنید .
7
در نوار ابزار Electric Field (emw) ، روی  Plot کلیک کنید .
Cut Line 3D 1
1
در نوار ابزار نتایج ، بر روی  Cut  Line  3D کلیک کنید .
2
در پنجره تنظیمات برای Cut  Line  3D ، بخش Line  Data را پیدا کنید .
3
در ردیف  1 ، X را روی -5e3 تنظیم کنید .
4
در ردیف  2 ، X را روی 5e3 تنظیم کنید .
گروه طرح 1 بعدی 2
1
در نوار ابزار نتایج ، روی  1D  Plot  Group کلیک کنید .
2
در پنجره Settings for 1D  Plot  Group ، بخش Data را پیدا کنید .
3
از لیست Dataset ، Cut  Line  3D  1 را انتخاب کنید .
4
برای گسترش بخش عنوان کلیک کنید . از لیست نوع عنوان  ، دستی را انتخاب کنید .
5
در قسمت متن عنوان ، هنجار میدان الکتریکی، مقیاس گزارش را تایپ کنید .
نمودار خطی 1
1
روی 1D  Plot  Group  کلیک راست کرده و Line  Graph را انتخاب کنید .
2
در پنجره تنظیمات برای نمودار خط  ، بخش y-Axis Data را پیدا کنید .
3
در قسمت متن Expression ، 20*log10(emw.normE) را تایپ کنید .
4
روی Replace  Expression در گوشه سمت راست بالای بخش x-Axis  Data کلیک کنید . از منو، Component   (comp1)>Geometry>Coordinate>x  –  x-coordinate را انتخاب کنید .
5
در نوار ابزار 1D Plot Group 2 ، روی  Plot کلیک کنید .
نموداری که ایجاد کردید باید از خط ثابت در شکل 2 تقلید کند . به مقایسه این نتیجه با مدلی بدون ذخیره هیدروکربنی ادامه دهید.
مواد
روغن (mat3)
در پنجره Model  Builder ، در قسمت Component   (comp1)>Materials روی Oil  (mat3) کلیک راست کرده و Disable را انتخاب کنید .
مطالعه 1
راه حل 1 (sol1)
1
در پنجره Model  Builder ، گره Study  1>Solver  Configurations را گسترش دهید .
2
روی Solution   (sol1) کلیک راست کرده و Disable را انتخاب کنید .
به این ترتیب یک گره حل کننده جداگانه برای مورد دوم دریافت می کنید.
3
در نوار ابزار صفحه اصلی ،  روی محاسبه کلیک کنید .
نتایج
چند برش
1
در پنجره Model  Builder ، گره Electric  Field  (emw)  1 را گسترش دهید ، سپس روی Multislice کلیک کنید .
2
در پنجره تنظیمات برای Multislice ، بخش Expression را پیدا کنید .
3
در قسمت متن Expression ، 20*log10(emw.normE) را تایپ کنید .
4
بخش Multiplane  Data را پیدا کنید . زیربخش X-planes را پیدا کنید . در قسمت متن Planes عدد 0 را تایپ کنید .
5
زیربخش Z-planes را پیدا کنید . در قسمت متن Planes عدد 0 را تایپ کنید .
6
در نوار ابزار Electric Field (emw) 1 ، روی  Plot کلیک کنید .
Cut Line 3D 2
1
در نوار ابزار نتایج ، بر روی  Cut  Line  3D کلیک کنید .
2
در پنجره تنظیمات برای Cut  Line  3D ، قسمت Data را پیدا کنید .
3
از فهرست مجموعه داده ، مطالعه  1/راه حل   (sol2) را انتخاب کنید .
4
قسمت Line  Data را پیدا کنید . در ردیف  1 ، X را روی -5e3 تنظیم کنید .
5
در ردیف  2 ، X را روی 5e3 تنظیم کنید .
نمودار خط 2
1
در پنجره Model  Builder ، روی 1D  Plot  Group  2 کلیک راست کرده و Line  Graph را انتخاب کنید .
2
در پنجره تنظیمات برای نمودار خط  ، بخش داده را پیدا کنید .
3
از لیست Dataset ، Cut  Line  3D  2 را انتخاب کنید .
4
قسمت y-Axis  Data را پیدا کنید . در قسمت متن Expression ، 20*log10(emw.normE) را تایپ کنید .
5
روی Replace  Expression در گوشه سمت راست بالای بخش x-Axis  Data کلیک کنید . از منو، Component   (comp1)>Geometry>Coordinate>x  –  x-coordinate را انتخاب کنید .
6
برای گسترش بخش Coloring  and  Style کلیک کنید . زیربخش Line  style را پیدا کنید . از لیست Line ، Dashed را انتخاب کنید .
7
در نوار ابزار 1D Plot Group 2 ، روی  Plot کلیک کنید .