Qazma quyularında çelik qoruyucu boru, çevresindeki sement kılıfı ve qaya formasiyaı birlikte çalışan çok katmanlı bir taşıyıcı ve yalıtım sistemi oluşturur. Muhafaza borusunun quyu merkezinden uzaklaşması veya sement ile kaya arasındaki bağın zayıflaması, gərginlikleri belirli bölgelerde yoğunlaştırarak çatlaqların daha erken başlamasına, farklı yönlere sapmasına ve maye geçişine elverişli zədə 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 sahəsi qırılma modeli geliştirdi. Model önce analitik çözümler ve daha önce yayımlanmış karşılaştırma problemleriyle doğrulandı. Ardından qoruyucu boru eksantrikliksi, sement-kaya arayüzünün möhkəmlikı ve üçölçülü çatlaq etkileşimleri ayrı sayısal deneylerle araştırıldı.

Sonuçlara göre qoruyucu boru mərkəzdə olarkən çatlaq başlanğıc təzyiqi 71 MPa iken, yüzde 83,33 eksantriklikde bu değer 49,5 MPa’a düştü. Bu yaklaşık yüzde 30’luk bir möhkəmlik kaybına karşılık geliyor. Eksantrisite yüzde 50’ye ulaştığında qaya formasiyaı içinde radial olmayan, eğimli sürüşmə çatlaqları da ortaya çıktı. Çimento-kaya arayüzünün çevredeki malzemelerin möhkəmlikının yaklaşık yüzde 30’una düşmesi halinde ise çatlaqlar kayaya geçmek yerine arayüz boyunca yön değiştirdi ve sement içinde daha fazla radial çatlaq oluştu.

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

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

Ancak qoruyucu borunun quyu merkezine tam olarak yerleşmemesi, sement kalınlığının çevresel yönde değişmesine neden olur. İnce sement bölgesi daha yüksek gərginliklere maruz kalabilir. Çimento büzülmesi, qazma çamuru kalıntıları, boşluklar, sıcaklık değişimleri ve yetersiz sementlama da sement-kaya veya çelik-sement arayüzlerini zayıflatabilir.

Oluşabilecek başlıca zədəlar radial çatlaqlar, sürüşmə çatlaqları, çelik-sement ayrılması, sement-kaya ayrılması ve quyu eksenine dik disk biçimli qırılmalardır. Bu kusurlar birleştiğinde mikrohalka olarak adlandırılan maye kanalları meydana gelebilir ve quyunun hidrolik yalıtımı bozulabilir.

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

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

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

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

Arayüz möhkəmlikı nasıl değiştirildi?

Çimento ile kaya arasındaki bağın möhkəmlikı κ parametresiyle diaqnostikamlandı. Etkin arayüz qı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 qı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 qırılma enerjisi nötr durumun yüzde 30’una iner.

Model, çekme baskın bölgelerde rijitliği azaltırken təzyiq altında fiziksel olmayan çatlaq ilerlemesini engelleyen bir çekme-təzyiq 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 eksantrikliksi

Muhafaza borusunun mərkəzdən uzaklığı ölçüsuz eksantriklik parametresiyle diaqnostikamlandı:

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

Burada r_e qoruyucu boru merkezi ile quyu merkezi arasındaki uzaklığı, R₁ ve R₂ ise ilgili iç ve dış yarıçapları gösterir. e = 0 eş mərkəzli sistemi, yüksek e değerleri ise sement 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 gərginlik çözümü

İlk doğrulamada iç yüzeyine 30 MPa, dış yüzeyine 50 MPa təzyiq uygulanan eksantrik delikli bir sement borusu modellenerek teğetsel gərginlikler 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ürüşməsı 10 mm ve eksantriklik yüzde 40 olarak belirlendi.

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

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

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

Ç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 zədə 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, sement ve kaya içeren quyu geometrisi

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

Model; sement içinde radial çatlaqları, çatlaqların qaya formasiyaına geçişini, çelik-sement arayüzündeki mikroçatlaqları ve ayrılmayı yayımlanmış referans çalışmayla uyumlu biçimde yeniden üretti. Çeliğin çok yüksek qırılma enerjisi nedeniyle çatlaq çelik gövdeye giremedi ve çelik-sement sınırı boyunca yön değiştirdi.

Eksantrisite arttıkça quyu möhkəmlikı azaldı

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

Eksantrisite yüzde 0’dan yüzde 83,33’e yükseldiğinde çatlaq başlanğıc təzyiqi 21,5 MPa azaldı. Eş mərkəzli geometride teğetsel gərginlik sement çevresinde büyük ölçüde düzgün dağılırken, eksantrik geometrilerde en ince sement bölgesinde yoğunlaştı. İlk çatlaq da sürekli olarak bu ince bölgede başladı.

Basınç 97 MPa’a çıkarıldığında eş mərkəzli quyuda simetrik çatlaq deseni korunurken, eksantriklik arttıkça zədə belirgin biçimde asimetrik hale geldi. İlk çatlağın ardından daha kalın sement bölgelerinde saat yönünün tersine ilerleyen sırayla yeni çatlaqlar oluştu.

Yüzde 50 eksantriklik yeni bir qırılma türünü ortaya çıkardı

Eksantrisite yüzde 50’ye ulaştığında kaya içinde radial doğrultudan sapan eğimli çatlaqlar gözlendi. Yaklaşık 10 derecelik azimut bölgesinde başlayan bu çatlaqlar, yüksek yerel sürüşmə gərginlikleriyle ilişkiliydi. Eş mərkəzli sistemde görülmeyen bu qırılma modu, geometrik kusurun yalnızca qı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 qırılma enerjisine sahip olması anlamına geliyor.

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

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

Bu sonuç, arayüz zayıflığının quyu 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 sement içinde daha fazla zədə ve olası sızıntı yolu oluşturabilir.

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

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

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

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

Model; sementlama tasarımı, qoruyucu boru merkezleyici yerleşimi, sement bağ loglarının yorumlanması ve quyu təzyiq 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 quyu geometrileri, malzeme özellikleri, sıcaklık geçmişi, gözenek basıncı ve maye 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 sahəsi qırılma modeli, zayıf arayüzleri içeren çok malzemeli sonlu element sistemlerine uygulanıyor. Çözüm, yer değiştirme ve zədə alanlarının sırayla çözüldüğü kademeli bir iteratif algoritmayla gerçekleştiriliyor.

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

Ana sayısal deneylerde kullanılan malzemeler

Temel bulgular

• Muhafaza borusu eksantrikliksi yüzde 0’dan yüzde 83,33’e çıktığında çatlaq başlanğıc təzyiqi 71 MPa’dan 49,5 MPa’a düştü.
• İlk çatlaq, eksantrik quyularda sürekli olarak sement kılıfının en ince bölgesinde başladı.
• Eksantrisite yüzde 50 ve üzerine çıktığında qaya formasiyaında yüksek sürüşmə gərginliksinin bulunduğu bölgelerde eğimli çatlaqlar oluştu.
• Çimento-kaya arayüzünün etkin möhkəmlikı nötr durumun yüzde 30’una veya altına indiğinde çatlaq geçişinden çok arayüz boyunca sapma ve ayrılma baskın hale geldi.
• Zayıf arayüz, sement içinde ardışık yeni radial çatlaqları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 gərginlik rahatlaması oluşturarak kaya geçişini geciktirdi.
• Üç ölçülu model, farklı derinliklerde çatlaq başlangıcını, gərginlik gölgesini, eksenden sapan çatlaq böyüməsini 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 sahəsi varsayımlarına dayanıyor.
• Arayüzde yalnızca açılma baskın zədə modelleniyor; sürtünmeli sürüşmə hesaba katılmıyor.
• Akışkan akışı, gözenek elastisitesi, termal etkiler ve hidro-mekanik bağlantı modele dahil edilmedi.
• İki ölçülu sonuçların çoğu düzlem forma değiştirme varsayımına dayanıyor.
• Üç ölçülu model 100 mm derinliğinde bir gösterim problemi olup gerçek saha ölçeğini doğrudan temsil etmiyor.
• Analitik gərginlik doğrulamasındaki en yüksek bağıl hata yüzde 11,37’dir.
• Hesaplanan təzyiq eşikleri seçilen geometriye, malzeme özelliklerine ve faz sahəsi uzunluq ölçeğine bağlıdır; bütün quyular için evrensel sınırlar olarak yorumlanmamalıdır.
• Olası sızıntı yolları ve sürekli muhafaza basıncı riski mekanik çatlaq geometrisinden çıkarılmaktadır; çalışmada doğrudan maye 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.