Бургулоо кудукlarında çelik каптама түтүк, çevresindeki цемент kılıfı ve тоо тектүү формация birlikte çalışan çok katmanlı bir taşıyıcı ve yalıtım sistemi oluşturur. Muhafaza borusunun кудук merkezinden uzaklaşması veya цемент ile kaya arasındaki bağın zayıflaması, чыңалууleri belirli bölgelerde yoğunlaştırarak жаракаların daha erken başlamasına, farklı yönlere sapmasına ve суюктук geçişine elverişli зыян yolları oluşmasına neden olabilir.

Araştırmacılar bu süreci incelemek için zayıf malzeme интерфейсlerini de kapsayan hibrit bir фаза талаасы сынуу modeli geliştirdi. Model önce analitik çözümler ve daha önce yayımlanmış karşılaştırma problemleriyle doğrulandı. Ardından каптама түтүк эксцентриситетsi, цемент-kaya интерфейсünün бекемдикı ve үч өлчөмдүү жарака etkileşimleri ayrı sayısal deneylerle araştırıldı.

Sonuçlara göre каптама түтүк борбордо турганда жарака башталуу басымы 71 MPa iken, yüzde 83,33 эксцентриситетde bu değer 49,5 MPa’a düştü. Bu yaklaşık yüzde 30’luk bir бекемдик kaybına karşılık geliyor. Eksantrisite yüzde 50’ye ulaştığında тоо тектүү формация içinde радиалдык olmayan, eğimli жылма жаракаları da ortaya çıktı. Çimento-kaya интерфейсünün çevredeki malzemelerin бекемдикının yaklaşık yüzde 30’una düşmesi halinde ise жаракаlar kayaya geçmek yerine интерфейс boyunca yön değiştirdi ve цемент içinde daha fazla радиалдык жарака oluştu.

Kuyu bütünlüğü neden önemlidir?

Çimentolanmış кудукlar yalnızca petrol ve doğal gaz üretiminde kullanılmaz. Jeotermal enerji, yer altı hidrojen ve doğal gaz depolama, karbondioksit enjeksiyonu, yer altı atık depolama, yer altı suyu üretimi ve çeşitli madencilik uygulamaları da benzer кудук mimarisine dayanır. Bu sistemlerde çelik каптама түтүк mekanik taşıyıcılığı sağlarken цемент kılıfı, boru ile kaya arasındaki boşluğu doldurur ve farklı jeolojik tabakalar arasında суюктук geçişini sınırlar.

Ancak каптама түтүкnun кудук merkezine tam olarak yerleşmemesi, цемент kalınlığının çevresel yönde değişmesine neden olur. İnce цемент bölgesi daha yüksek чыңалууlere maruz kalabilir. Çimento büzülmesi, бургулоо çamuru kalıntıları, boşluklar, sıcaklık değişimleri ve yetersiz цементlama da цемент-kaya veya çelik-цемент интерфейсlerini zayıflatabilir.

Oluşabilecek başlıca зыянlar радиалдык жаракаlar, жылма жаракаları, çelik-цемент ayrılması, цемент-kaya ayrılması ve кудук eksenine dik disk biçimli сынууlardır. Bu kusurlar birleştiğinde mikrohalka olarak adlandırılan суюктук kanalları meydana gelebilir ve кудукnun hidrolik yalıtımı bozulabilir.

Faz alanı yöntemi çatlağı nasıl temsil ediyor?

Geleneksel чектүү элемент modellerinde çatlağın izleyeceği yolun önceden диагностикаmlanması gerekebilir. Faz alanı yaklaşımında ise жарака, keskin bir geometrik yüzey olarak değil, malzemenin sağlam durumdan зыянlı duruma geçtiği yayılı bir bölge olarak temsil edilir.

Modelde d фаза талаасы değişkeni malzemenin зыян durumunu gösterir. Yaklaşık olarak d = 0 sağlam malzemeyi, d = 1 ise tamamen çatlamış bölgeyi temsil eder. Malzemenin rijitliği şu fonksiyonla azaltılır:

\[ g(d)=(1-d)^2 \]

Bu sayede жарака başlangıcı, dallanması, birleşmesi ve farklı malzemeler arasındaki yön değişimleri önceden çizilmiş bir жарака hattına ihtiyaç duyulmadan enerji minimizasyonu sonucunda oluşabilir.

Toplam potansiyel enerji; iki malzemenin elastik форма değiştirme enerjisini, hacim içindeki düzenlileştirilmiş жарака yüzeyi enerjisini ve malzeme интерфейсüne özgü сынуу enerjisini birlikte içerir:

\[ U = \sum_{m=1}^{2}\int_{\Omega^m}g(d^m)\psi^m(\varepsilon(u^m))\,d\Omega + \sum_{m=1}^{2}\int_{\Omega^m}G_c^m \frac{1}{2} \left[ \frac{(d^m)^2}{\ell} + \ell|\nabla d^m|^2 \right]d\Omega + \int_{\Gamma_{\mathrm{int}}}G_c^{\mathrm{ib}}\tilde{d}^{\,2}\,d\Gamma \]

Burada u yer değiştirme alanını, ε küçük форма değiştirme tensörünü, ψ elastik enerji yoğunluğunu, G_c kritik enerji salım hızını, ℓ фаза талаасы узундук ölçeğini ve Γ_int malzeme интерфейсünü ifade eder.

Arayüz бекемдикı nasıl değiştirildi?

Çimento ile kaya arasındaki bağın бекемдикı κ parametresiyle диагностикаmlandı. Etkin интерфейс сынуу enerjisi şu bağıntıyla ifade edildi:

\[ G_c^{\mathrm{eff}}=G_c^{\mathrm{nb}}(1+\kappa) \]

G_c^nb nötr bağlı интерфейсün сынуу enerjisidir. κ = 0 nötr bağı, κ > 0 güçlendirilmiş bağı, κ < 0 ise zayıflatılmış bağı temsil eder. Örneğin κ = −0,7 olduğunda интерфейсün etkin сынуу enerjisi nötr durumun yüzde 30’una iner.

Model, çekme baskın bölgelerde rijitliği azaltırken басым altında fiziksel olmayan жарака ilerlemesini engelleyen bir çekme-басым enerji ayrımı kullanıyor. Hasarın geri kapanmasını veya iyileşmesini önlemek için geçmişte ulaşılan en yüksek çekme enerjisini saklayan bir geçmiş alanı da hesaba katılıyor.

Muhafaza borusu эксцентриситетsi

Muhafaza borusunun борбордон uzaklığı өлчөмsuz эксцентриситет parametresiyle диагностикаmlandı:

\[ e=\frac{r_e}{R_2-R_1} \]

Burada r_e каптама түтүк merkezi ile кудук merkezi arasındaki uzaklığı, R₁ ve R₂ ise ilgili iç ve dış yarıçapları gösterir. e = 0 eş борборлуу sistemi, yüksek e değerleri ise цемент kılıfının bir tarafta belirgin biçimde inceldiği geometrileri temsil eder.

Model üç aşamada doğrulandı

1. Eksantrik delikli boruda analitik чыңалуу çözümü

İlk doğrulamada iç yüzeyine 30 MPa, dış yüzeyine 50 MPa басым uygulanan eksantrik delikli bir цемент borusu modellenerek teğetsel чыңалууler Jeffery’nin kapalı biçimli analitik çözümüyle karşılaştırıldı. Modelde iç yarıçap 69,5 mm, dış yarıçap 94,5 mm, merkez жылмаsı 10 mm ve эксцентриситет yüzde 40 olarak belirlendi.

Yaklaşık 7,6 milyon dörtgen чектүү элементın kullanıldığı çözüm, analitik чыңалуу dağılımını genel olarak yeniden üretti. En yüksek bağıl hata, yüksek чыңалуу gradyanının bulunduğu 90 derecelik konumda yüzde 11,37 olarak ölçüldü. Araştırmacılar bu farkı kullanılan birinci dereceden dörtgen элементların чыңалуу çözünürlüğüyle ilişkilendirdi.

2. Zayıf dairesel интерфейсe sahip iki malzemeli plaka

İkinci doğrulamada 540 × 540 mm өлчөмlarında, 50 mm başlangıç çatlağına ve zayıf dairesel интерфейсe sahip iki malzemeli bir plaka incelendi. Arayüz сынуу enerjisi 50 N/m ve zayıflatma parametresi −0,875 olarak диагностикаmlandı.

Çatlak önce birinci malzeme içinde açılma modunda ilerledi, dairesel zayıf интерфейсe ulaştığında интерфейс boyunca saptı ve daha sonra ikinci malzemeye yeniden girdi. Elde edilen зыян yolu, yayımlanmış referans simülasyonuyla güçlü niteliksel uyum gösterdi. Çatlağın интерфейсü terk ettiği noktalardaki küçük fark, referans modelde интерфейсün sonlu kalınlıkta üçüncü bir malzeme olarak, bu araştırmada ise sıfır kalınlıklı keskin bir yüzey olarak modellenmesine bağlandı.

3. Çelik, цемент ve kaya içeren кудук geometrisi

Üçüncü doğrulamada çelik каптама түтүк, цемент kılıfı ve Castlegate kumtaşı formasyonundan oluşan gerçekçi bir кудук kesiti modellendi. Simetri nedeniyle geometrinin dörtte biri kullanıldı ve yaklaşık 700 bin dörtgen элементla çözüldü.

Model; цемент içinde радиалдык жаракаları, жаракаların тоо тектүү формацияna geçişini, çelik-цемент интерфейсündeki mikroжаракаları ve ayrılmayı yayımlanmış referans çalışmayla uyumlu biçimde yeniden üretti. Çeliğin çok yüksek сынуу enerjisi nedeniyle жарака çelik gövdeye giremedi ve çelik-цемент sınırı boyunca yön değiştirdi.

Eksantrisite arttıkça кудук бекемдикı azaldı

Araştırmacılar yüzde 0 ile yüzde 83,33 arasında yedi farklı эксцентриситет değeri kullandı. Çatlak başlangıcındaki iç кудук basıncı aşağıdaki gibi değişti:

Eksantrisite yüzde 0’dan yüzde 83,33’e yükseldiğinde жарака башталуу басымы 21,5 MPa azaldı. Eş борборлуу geometride teğetsel чыңалуу цемент çevresinde büyük ölçüde düzgün dağılırken, eksantrik geometrilerde en ince цемент bölgesinde yoğunlaştı. İlk жарака da sürekli olarak bu ince bölgede başladı.

Basınç 97 MPa’a çıkarıldığında eş борборлуу кудукda simetrik жарака deseni korunurken, эксцентриситет arttıkça зыян belirgin biçimde asimetrik hale geldi. İlk çatlağın ardından daha kalın цемент bölgelerinde saat yönünün tersine ilerleyen sırayla yeni жаракаlar oluştu.

Yüzde 50 эксцентриситет yeni bir сынуу türünü ortaya çıkardı

Eksantrisite yüzde 50’ye ulaştığında kaya içinde радиалдык doğrultudan sapan eğimli жаракаlar gözlendi. Yaklaşık 10 derecelik azimut bölgesinde başlayan bu жаракаlar, yüksek yerel жылма чыңалууleriyle ilişkiliydi. Eş борборлуу sistemde görülmeyen bu сынуу modu, geometrik kusurun yalnızca сынуу basıncını değil çatlağın türünü ve yönünü de değiştirebildiğini gösterdi.

Zayıf интерфейс çatlağın yolunu değiştirdi

Çimento-kaya интерфейсü için κ değeri −0,5 ile −0,9 arasında değiştirildi. Arayüz zayıflamasının belirgin etkisi κ ≤ −0,7 düzeyinde ortaya çıktı. Bu koşul, интерфейсün nötr bağlı duruma göre yüzde 30 veya daha düşük etkin сынуу enerjisine sahip olması anlamına geliyor.

Daha güçlü интерфейсlerde цемент içinde başlayan радиалдык жарака doğrudan тоо тектүү формацияna geçti. Zayıf интерфейсlerde ise жарака kaya içine girmek yerine цемент-kaya sınırı boyunca yön değiştirdi. Çatlağın интерфейсde tutulması, цемент kılıfındaki чыңалууnin hemen boşalmasını engelledi ve yeni радиалдык жаракаların oluşmasına ortam hazırladı.

Arayüzün yüzde 80 zayıflatıldığı κ = −0,8 durumunda yük-yer değiştirme eğrilerinde art arda oluşan çok sayıda сынуу sıçraması görüldü. Nötr bağlı durumda ise iki радиалдык жарака yaklaşık eş zamanlı başladı ve eğride tek baskın сынуу noktası oluştu.

En zayıf интерфейс her zaman kayanın en erken çatladığı durum olmadı. κ = −0,8 için жарака тоо тектүү формацияna 123,5 MPa’da geçerken, κ = −0,9 için geçiş 144 MPa’a kadar gecikti. Bunun nedeni yüzde 90 zayıflatılmış интерфейсün daha erken ve daha kapsamlı biçimde ayrılması, цемент içindeki чыңалууyi yerel olarak rahatlatması ve kaya sistemini mekanik olarak kısmen ayırmasıydı.

Bu sonuç, интерфейс zayıflığının кудук davranışını tek yönlü ve doğrusal biçimde değiştirmediğini gösteriyor. Daha zayıf bağ, интерфейс ayrılmasını kolaylaştırırken bazı koşullarda çatlağın kayaya geçişini geciktirebilir; buna karşılık цемент içinde daha fazla зыян ve olası sızıntı yolu oluşturabilir.

Üç өлчөмlu model düzlem форма değiştirme varsayımının göremediği süreçleri gösterdi

Son aşamada 100 mm derinliğe sahip үч өлчөмдүү bir кудук modeli çözüldü. Çimento-kaya интерфейсü yüzde 90 zayıflatıldı. Başlangıç ağı 724 bin üç doğrusal sekiz düğümlü элементdan oluşuyordu. Hasar ilerledikçe uyarlamalı ağ inceltme etkinleştirildi.

İlk жаракаların çevresinde oluşan чыңалуу gölgeleri, daha sonraki жаракаların tüm кудук derinliği boyunca aynı anda başlamasını engelledi. Bu nedenle жаракаlar farklı derinliklerde çekirdeklendi, eksenden saptı ve интерфейс boyunca birleşti. Bu davranışların эки өлчөмдүү düzlem форма değiştirme modelleriyle temsil edilmesi mümkün değildir.

Hesaplama başlangıçta yaklaşık 400 bin serbestlik derecesine sahipken uyarlamalı ağ inceltme sonucunda 113 milyon serbestlik derecesine ulaştı. Simülasyon, 2,45 GHz AMD 64 çekirdekli işlemcilerden oluşan bir kümede yaklaşık 20 bin çekirdek-saat gerektirdi. Araştırmacılara göre aynı çözümün tekdüze ince ağla gerçekleştirilmesi yaklaşık 12 milyar serbestlik derecesi gerektirecekti.

Mühendislik açısından ne anlama geliyor?

Çalışma, кудук merkezleme kalitesinin yalnızca цемент tüketimi veya montaj toleransı açısından değerlendirilmemesi gerektiğini gösteriyor. Muhafaza borusunun борбордон жылмаsı, çatlağın başladığı basıncı azaltabiliyor ve belirli bir eşikten sonra kaya içinde yeni жылма сынууsı biçimleri oluşturabiliyor.

Çimento-kaya bağının бекемдикı da yalnızca “güçlü” veya “zayıf” şeklinde değerlendirilmemeli. Arayüz бекемдикı; çatlağın kayaya geçip geçmeyeceğini, интерфейс boyunca sapıp sapmayacağını, цемент içinde kaç yeni жарака oluşacağını ve зыянın derinlik boyunca nasıl birleşeceğini belirleyebiliyor.

Model; цементlama tasarımı, каптама түтүк merkezleyici yerleşimi, цемент bağ loglarının yorumlanması ve кудук басым sınırlarının değerlendirilmesi için bir hesaplama altyapısı sunuyor. Ancak sonuçlar doğrudan saha güvenlik katsayısı olarak kullanılmadan önce gerçek кудук geometrileri, malzeme özellikleri, sıcaklık geçmişi, gözenek basıncı ve суюктук etkileriyle yeniden kalibre edilmelidir.

Çalışmanın Yöntemi ve Bulguları

Araştırma tasarımı

Araştırma fiziksel deney yerine doğrulanmış hesaplamalı deneyler kullanıyor. Hibrit фаза талаасы сынуу modeli, zayıf интерфейсleri içeren çok malzemeli чектүү элемент sistemlerine uygulanıyor. Çözüm, yer değiştirme ve зыян alanlarının sırayla çözüldüğü kademeli bir iteratif algoritmayla gerçekleştiriliyor.

• İki өлчөмlu analizlerde düzlem форма değiştirme varsayımı kullanıldı.
• İki өлчөмlu ağlarda çift doğrusal dörtgen, үч өлчөмдүү ağda üç doğrusal altı yüzlü элементlar kullanıldı.
• Ardışık iki iterasyondaki düğüm зыян vektörlerinin L2 norm farkı için yakınsama toleransı 10⁻² olarak belirlendi.
• Görselleştirmelerde зыян değişkeni 0,9’dan büyük olan элементlar çatlamış bölge olarak kabul edildi.
• Model deal.II açık kaynaklı чектүү элемент kütüphanesinde uygulandı.
• Üç өлчөмlu çözümde paralel hesaplama ve uyarlamalı ağ inceltme kullanıldı.

Ana sayısal deneylerde kullanılan malzemeler

Temel bulgular

• Muhafaza borusu эксцентриситетsi yüzde 0’dan yüzde 83,33’e çıktığında жарака башталуу басымы 71 MPa’dan 49,5 MPa’a düştü.
• İlk жарака, eksantrik кудукlarda sürekli olarak цемент kılıfının en ince bölgesinde başladı.
• Eksantrisite yüzde 50 ve üzerine çıktığında тоо тектүү формацияnda yüksek жылма чыңалууsinin bulunduğu bölgelerde eğimli жаракаlar oluştu.
• Çimento-kaya интерфейсünün etkin бекемдикı nötr durumun yüzde 30’una veya altına indiğinde жарака geçişinden çok интерфейс boyunca sapma ve ayrılma baskın hale geldi.
• Zayıf интерфейс, цемент içinde ardışık yeni радиалдык жаракаların oluşmasını destekledi.
• En fazla zayıflatılmış интерфейс, kayanın en erken çatladığı durum olmadı; erken tam ayrılma bazı bölgelerde чыңалуу rahatlaması oluşturarak kaya geçişini geciktirdi.
• Üç өлчөмlu model, farklı derinliklerde жарака başlangıcını, чыңалуу gölgesini, eksenden sapan жарака өсүүsini ve интерфейс boyunca birleşmeyi ortaya çıkardı.

Sınırlılıklar

• Çalışma hakem değerlendirmesinden geçmemiş bir preprinttir.
• Sonuçlar fiziksel saha deneyi değil, doğrulanmış sayısal model sonuçlarıdır.
• Model doğrusal izotropik elastisite ve kırılgan фаза талаасы varsayımlarına dayanıyor.
• Arayüzde yalnızca açılma baskın зыян modelleniyor; sürtünmeli жылма hesaba katılmıyor.
• Akışkan akışı, gözenek elastisitesi, termal etkiler ve hidro-mekanik bağlantı modele dahil edilmedi.
• İki өлчөмlu sonuçların çoğu düzlem форма değiştirme varsayımına dayanıyor.
• Üç өлчөмlu model 100 mm derinliğinde bir gösterim problemi olup gerçek saha ölçeğini doğrudan temsil etmiyor.
• Analitik чыңалуу doğrulamasındaki en yüksek bağıl hata yüzde 11,37’dir.
• Hesaplanan басым eşikleri seçilen geometriye, malzeme özelliklerine ve фаза талаасы узундук ölçeğine bağlıdır; bütün кудукlar için evrensel sınırlar olarak yorumlanmamalıdır.
• Olası sızıntı yolları ve sürekli muhafaza basıncı riski mekanik жарака geometrisinden çıkarılmaktadır; çalışmada doğrudan суюктук taşınımı simüle edilmemiştir.

Kaynak ve Yöntem Notu

Özgün çalışma: “Phase-field analysis of fracture in heterogeneous wellbore systems: effects of casing eccentricity and cement-formation interface strength.”

Yazarlar: Tharunsarathy Sachithanantham, Wasim Niyaz Munshi, Chandrasekhar Annavarapu ve Birendra Jha.

Sorumlu yazar: Chandrasekhar Annavarapu. Eş katkı beyanı bulunmamaktadır.

Kurumlar: Hindistan Teknoloji Enstitüsü Madras, İnşaat Mühendisliği Bölümü; Güney Kaliforniya Üniversitesi, Mork Ailesi Kimya Mühendisliği ve Malzeme Bilimi Bölümü.

Yıl ve kaynak türü: 2026 tarihli özgün sayısal mühendislik araştırması ve preprint.

Hakemlik durumu: Çalışma henüz hakem değerlendirmesinden geçmemiştir. Sonuçlar yayımlanmış nihai dergi makalesi gibi değerlendirilmemelidir.

Platform: SSRN.

DOI:10.2139/ssrn.6947717

Resmî çalışma bağlantısı:SSRN çalışma kaydı

Bu içerikteki bilimsel yöntem, sayısal değerler, denklemler, malzeme özellikleri, görsel yorumları, bulgular ve sınırlılıklar yalnızca incelenen çalışmaya dayanmaktadır. Dış kaynaklar yalnızca yazar kimliği, kurum, platform, yayımlanma tarihi ve DOI bilgilerinin bibliyografik doğrulanması için kullanılmış; çalışmada bulunmayan bilimsel sonuç eklenmemiştir.