Sondaj kuyularında çelik muhafaza borusu, çevresindeki çimento kılıfı ve kaya formasyonu birlikte çalışan çok katmanlı bir taşıyıcı ve yalıtım sistemi oluşturur. Muhafaza borusunun kuyu merkezinden uzaklaşması veya çimento ile kaya arasındaki bağın zayıflaması, gerilmeleri belirli bölgelerde yoğunlaştırarak çatlakların daha erken başlamasına, farklı yönlere sapmasına ve akışkan geçişine elverişli hasar yolları oluşmasına neden olabilir.

Araştırmacılar bu süreci incelemek için zayıf malzeme arayüzlerini de kapsayan hibrit bir faz alanı kırılma modeli geliştirdi. Model önce analitik çözümler ve daha önce yayımlanmış karşılaştırma problemleriyle doğrulandı. Ardından muhafaza borusu eksantrisitesi, çimento-kaya arayüzünün dayanımı ve üç boyutlu çatlak etkileşimleri ayrı sayısal deneylerle araştırıldı.

Sonuçlara göre muhafaza borusu merkezdeyken çatlak başlangıç basıncı 71 MPa iken, yüzde 83,33 eksantrisitede bu değer 49,5 MPa’a düştü. Bu yaklaşık yüzde 30’luk bir dayanım kaybına karşılık geliyor. Eksantrisite yüzde 50’ye ulaştığında kaya formasyonu içinde radyal olmayan, eğimli kayma çatlakları da ortaya çıktı. Çimento-kaya arayüzünün çevredeki malzemelerin dayanımının yaklaşık yüzde 30’una düşmesi halinde ise çatlaklar kayaya geçmek yerine arayüz boyunca yön değiştirdi ve çimento içinde daha fazla radyal çatlak oluştu.

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

Çimentolanmış kuyular 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 kuyu mimarisine dayanır. Bu sistemlerde çelik muhafaza borusu mekanik taşıyıcılığı sağlarken çimento kılıfı, boru ile kaya arasındaki boşluğu doldurur ve farklı jeolojik tabakalar arasında akışkan geçişini sınırlar.

Ancak muhafaza borusunun kuyu merkezine tam olarak yerleşmemesi, çimento kalınlığının çevresel yönde değişmesine neden olur. İnce çimento bölgesi daha yüksek gerilmelere maruz kalabilir. Çimento büzülmesi, sondaj çamuru kalıntıları, boşluklar, sıcaklık değişimleri ve yetersiz çimentolama da çimento-kaya veya çelik-çimento arayüzlerini zayıflatabilir.

Oluşabilecek başlıca hasarlar radyal çatlaklar, kayma çatlakları, çelik-çimento ayrılması, çimento-kaya ayrılması ve kuyu eksenine dik disk biçimli kırılmalardır. Bu kusurlar birleştiğinde mikrohalka olarak adlandırılan akışkan kanalları meydana gelebilir ve kuyunun hidrolik yalıtımı bozulabilir.

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

Geleneksel sonlu eleman modellerinde çatlağın izleyeceği yolun önceden tanımlanması gerekebilir. Faz alanı yaklaşımında ise çatlak, keskin bir geometrik yüzey olarak değil, malzemenin sağlam durumdan hasarlı duruma geçtiği yayılı bir bölge olarak temsil edilir.

Modelde d faz alanı değişkeni malzemenin hasar 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 çatlak başlangıcı, dallanması, birleşmesi ve farklı malzemeler arasındaki yön değişimleri önceden çizilmiş bir çatlak hattına ihtiyaç duyulmadan enerji minimizasyonu sonucunda oluşabilir.

Toplam potansiyel enerji; iki malzemenin elastik şekil değiştirme enerjisini, hacim içindeki düzenlileştirilmiş çatlak yüzeyi enerjisini ve malzeme arayüzüne özgü kırılma 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 şekil değiştirme tensörünü, ψ elastik enerji yoğunluğunu, G_c kritik enerji salım hızını, ℓ faz alanı uzunluk ölçeğini ve Γ_int malzeme arayüzünü ifade eder.

Arayüz dayanımı nasıl değiştirildi?

Çimento ile kaya arasındaki bağın dayanımı κ parametresiyle tanımlandı. Etkin arayüz kırılma enerjisi şu bağıntıyla ifade edildi:

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

G_c^nb nötr bağlı arayüzün kırılma enerjisidir. κ = 0 nötr bağı, κ > 0 güçlendirilmiş bağı, κ < 0 ise zayıflatılmış bağı temsil eder. Örneğin κ = −0,7 olduğunda arayüzün etkin kırılma enerjisi nötr durumun yüzde 30’una iner.

Model, çekme baskın bölgelerde rijitliği azaltırken basınç altında fiziksel olmayan çatlak ilerlemesini engelleyen bir çekme-basınç 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 eksantrisitesi

Muhafaza borusunun merkezden uzaklığı boyutsuz eksantrisite parametresiyle tanımlandı:

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

Burada r_e muhafaza borusu merkezi ile kuyu merkezi arasındaki uzaklığı, R₁ ve R₂ ise ilgili iç ve dış yarıçapları gösterir. e = 0 eş merkezli sistemi, yüksek e değerleri ise çimento 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 gerilme çözümü

İlk doğrulamada iç yüzeyine 30 MPa, dış yüzeyine 50 MPa basınç uygulanan eksantrik delikli bir çimento borusu modellenerek teğetsel gerilmeler 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 kayması 10 mm ve eksantrisite yüzde 40 olarak belirlendi.

Yaklaşık 7,6 milyon dörtgen sonlu elemanın kullanıldığı çözüm, analitik gerilme dağılımını genel olarak yeniden üretti. En yüksek bağıl hata, yüksek gerilme 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 elemanların gerilme çözünürlüğüyle ilişkilendirdi.

2. Zayıf dairesel arayüze sahip iki malzemeli plaka

İkinci doğrulamada 540 × 540 mm boyutlarında, 50 mm başlangıç çatlağına ve zayıf dairesel arayüze sahip iki malzemeli bir plaka incelendi. Arayüz kırılma enerjisi 50 N/m ve zayıflatma parametresi −0,875 olarak tanımlandı.

Çatlak önce birinci malzeme içinde açılma modunda ilerledi, dairesel zayıf arayüze ulaştığında arayüz boyunca saptı ve daha sonra ikinci malzemeye yeniden girdi. Elde edilen hasar yolu, yayımlanmış referans simülasyonuyla güçlü niteliksel uyum gösterdi. Çatlağın arayüzü terk ettiği noktalardaki küçük fark, referans modelde arayüzü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, çimento ve kaya içeren kuyu geometrisi

Üçüncü doğrulamada çelik muhafaza borusu, çimento kılıfı ve Castlegate kumtaşı formasyonundan oluşan gerçekçi bir kuyu kesiti modellendi. Simetri nedeniyle geometrinin dörtte biri kullanıldı ve yaklaşık 700 bin dörtgen elemanla çözüldü.

Model; çimento içinde radyal çatlakları, çatlakların kaya formasyonuna geçişini, çelik-çimento arayüzündeki mikroçatlakları ve ayrılmayı yayımlanmış referans çalışmayla uyumlu biçimde yeniden üretti. Çeliğin çok yüksek kırılma enerjisi nedeniyle çatlak çelik gövdeye giremedi ve çelik-çimento sınırı boyunca yön değiştirdi.

Eksantrisite arttıkça kuyu dayanımı azaldı

Araştırmacılar yüzde 0 ile yüzde 83,33 arasında yedi farklı eksantrisite değeri kullandı. Çatlak başlangıcındaki iç kuyu basıncı aşağıdaki gibi değişti:

Eksantrisite yüzde 0’dan yüzde 83,33’e yükseldiğinde çatlak başlangıç basıncı 21,5 MPa azaldı. Eş merkezli geometride teğetsel gerilme çimento çevresinde büyük ölçüde düzgün dağılırken, eksantrik geometrilerde en ince çimento bölgesinde yoğunlaştı. İlk çatlak da sürekli olarak bu ince bölgede başladı.

Basınç 97 MPa’a çıkarıldığında eş merkezli kuyuda simetrik çatlak deseni korunurken, eksantrisite arttıkça hasar belirgin biçimde asimetrik hale geldi. İlk çatlağın ardından daha kalın çimento bölgelerinde saat yönünün tersine ilerleyen sırayla yeni çatlaklar oluştu.

Yüzde 50 eksantrisite yeni bir kırılma türünü ortaya çıkardı

Eksantrisite yüzde 50’ye ulaştığında kaya içinde radyal doğrultudan sapan eğimli çatlaklar gözlendi. Yaklaşık 10 derecelik azimut bölgesinde başlayan bu çatlaklar, yüksek yerel kayma gerilmeleriyle ilişkiliydi. Eş merkezli sistemde görülmeyen bu kırılma modu, geometrik kusurun yalnızca kırılma basıncını değil çatlağın türünü ve yönünü de değiştirebildiğini gösterdi.

Zayıf arayüz çatlağın yolunu değiştirdi

Çimento-kaya arayüzü 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, arayüzün nötr bağlı duruma göre yüzde 30 veya daha düşük etkin kırılma enerjisine sahip olması anlamına geliyor.

Daha güçlü arayüzlerde çimento içinde başlayan radyal çatlak doğrudan kaya formasyonuna geçti. Zayıf arayüzlerde ise çatlak kaya içine girmek yerine çimento-kaya sınırı boyunca yön değiştirdi. Çatlağın arayüzde tutulması, çimento kılıfındaki gerilmenin hemen boşalmasını engelledi ve yeni radyal çatlakları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 kırılma sıçraması görüldü. Nötr bağlı durumda ise iki radyal çatlak yaklaşık eş zamanlı başladı ve eğride tek baskın kırılma noktası oluştu.

En zayıf arayüz her zaman kayanın en erken çatladığı durum olmadı. κ = −0,8 için çatlak kaya formasyonuna 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ış arayüzün daha erken ve daha kapsamlı biçimde ayrılması, çimento içindeki gerilmeyi yerel olarak rahatlatması ve kaya sistemini mekanik olarak kısmen ayırmasıydı.

Bu sonuç, arayüz zayıflığının kuyu davranışını tek yönlü ve doğrusal biçimde değiştirmediğini gösteriyor. Daha zayıf bağ, arayüz ayrılmasını kolaylaştırırken bazı koşullarda çatlağın kayaya geçişini geciktirebilir; buna karşılık çimento içinde daha fazla hasar ve olası sızıntı yolu oluşturabilir.

Üç boyutlu model düzlem şekil değiştirme varsayımının göremediği süreçleri gösterdi

Son aşamada 100 mm derinliğe sahip üç boyutlu bir kuyu modeli çözüldü. Çimento-kaya arayüzü yüzde 90 zayıflatıldı. Başlangıç ağı 724 bin üç doğrusal sekiz düğümlü elemandan oluşuyordu. Hasar ilerledikçe uyarlamalı ağ inceltme etkinleştirildi.

İlk çatlakların çevresinde oluşan gerilme gölgeleri, daha sonraki çatlakların tüm kuyu derinliği boyunca aynı anda başlamasını engelledi. Bu nedenle çatlaklar farklı derinliklerde çekirdeklendi, eksenden saptı ve arayüz boyunca birleşti. Bu davranışların iki boyutlu düzlem şekil 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, kuyu merkezleme kalitesinin yalnızca çimento tüketimi veya montaj toleransı açısından değerlendirilmemesi gerektiğini gösteriyor. Muhafaza borusunun merkezden kayması, çatlağın başladığı basıncı azaltabiliyor ve belirli bir eşikten sonra kaya içinde yeni kayma kırılması biçimleri oluşturabiliyor.

Çimento-kaya bağının dayanımı da yalnızca “güçlü” veya “zayıf” şeklinde değerlendirilmemeli. Arayüz dayanımı; çatlağın kayaya geçip geçmeyeceğini, arayüz boyunca sapıp sapmayacağını, çimento içinde kaç yeni çatlak oluşacağını ve hasarın derinlik boyunca nasıl birleşeceğini belirleyebiliyor.

Model; çimentolama tasarımı, muhafaza borusu merkezleyici yerleşimi, çimento bağ loglarının yorumlanması ve kuyu basınç 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 kuyu geometrileri, malzeme özellikleri, sıcaklık geçmişi, gözenek basıncı ve akışkan 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 faz alanı kırılma modeli, zayıf arayüzleri içeren çok malzemeli sonlu eleman sistemlerine uygulanıyor. Çözüm, yer değiştirme ve hasar alanlarının sırayla çözüldüğü kademeli bir iteratif algoritmayla gerçekleştiriliyor.

• İki boyutlu analizlerde düzlem şekil değiştirme varsayımı kullanıldı.
• İki boyutlu ağlarda çift doğrusal dörtgen, üç boyutlu ağda üç doğrusal altı yüzlü elemanlar kullanıldı.
• Ardışık iki iterasyondaki düğüm hasar vektörlerinin L2 norm farkı için yakınsama toleransı 10⁻² olarak belirlendi.
• Görselleştirmelerde hasar değişkeni 0,9’dan büyük olan elemanlar çatlamış bölge olarak kabul edildi.
• Model deal.II açık kaynaklı sonlu eleman kütüphanesinde uygulandı.
• Üç boyutlu çözümde paralel hesaplama ve uyarlamalı ağ inceltme kullanıldı.

Ana sayısal deneylerde kullanılan malzemeler

Temel bulgular

• Muhafaza borusu eksantrisitesi yüzde 0’dan yüzde 83,33’e çıktığında çatlak başlangıç basıncı 71 MPa’dan 49,5 MPa’a düştü.
• İlk çatlak, eksantrik kuyularda sürekli olarak çimento kılıfının en ince bölgesinde başladı.
• Eksantrisite yüzde 50 ve üzerine çıktığında kaya formasyonunda yüksek kayma gerilmesinin bulunduğu bölgelerde eğimli çatlaklar oluştu.
• Çimento-kaya arayüzünün etkin dayanımı nötr durumun yüzde 30’una veya altına indiğinde çatlak geçişinden çok arayüz boyunca sapma ve ayrılma baskın hale geldi.
• Zayıf arayüz, çimento içinde ardışık yeni radyal çatlakların oluşmasını destekledi.
• En fazla zayıflatılmış arayüz, kayanın en erken çatladığı durum olmadı; erken tam ayrılma bazı bölgelerde gerilme rahatlaması oluşturarak kaya geçişini geciktirdi.
• Üç boyutlu model, farklı derinliklerde çatlak başlangıcını, gerilme gölgesini, eksenden sapan çatlak büyümesini ve arayüz 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 faz alanı varsayımlarına dayanıyor.
• Arayüzde yalnızca açılma baskın hasar modelleniyor; sürtünmeli kayma hesaba katılmıyor.
• Akışkan akışı, gözenek elastisitesi, termal etkiler ve hidro-mekanik bağlantı modele dahil edilmedi.
• İki boyutlu sonuçların çoğu düzlem şekil değiştirme varsayımına dayanıyor.
• Üç boyutlu model 100 mm derinliğinde bir gösterim problemi olup gerçek saha ölçeğini doğrudan temsil etmiyor.
• Analitik gerilme doğrulamasındaki en yüksek bağıl hata yüzde 11,37’dir.
• Hesaplanan basınç eşikleri seçilen geometriye, malzeme özelliklerine ve faz alanı uzunluk ölçeğine bağlıdır; bütün kuyular için evrensel sınırlar olarak yorumlanmamalıdır.
• Olası sızıntı yolları ve sürekli muhafaza basıncı riski mekanik çatlak geometrisinden çıkarılmaktadır; çalışmada doğrudan akışkan 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.