Bu çalışma, yenilenebilir enerji sistemleriyle birlikte giderek daha önemli hale gelen yer altı enerji depolama ve özellikle tuz mağaralarında hidrojen depolama konusunu jeomekanik güvenlik açısından incelemektedir. Güneş ve rüzgâr gibi yenilenebilir kaynaklar sürekli aynı miktarda enerji üretmez. Bu nedenle fazla elektriğin yeşil hidrojene dönüştürülmesi ve bu hidrojenin büyük hacimlerde güvenli şekilde depolanması, enerji sistemlerinin geleceği açısından kritik görülmektedir. Tuz mağaraları bu amaç için güçlü adaylardır; çünkü tuz kayası çok düşük geçirgenliğe, kendi kendini kısmen iyileştirebilme yeteneğine ve uygun mekanik davranışa sahiptir.

Ancak araştırmanın temel uyarısı şudur: Tuz mağaraları geleneksel olarak yılda bir veya az sayıda yavaş dolum-boşaltım döngüsüyle çalıştırıldığında güvenli kalabilir; fakat yenilenebilir enerji dalgalanmalarını dengelemek için daha sık hidrojen enjeksiyonu ve üretimi yapılırsa, mağara çevresindeki tuz kayası çok daha ağır bir termo-mekanik yorulmaya maruz kalır. Gaz enjeksiyonu sırasında sıkışma nedeniyle mağara içinde sıcaklık artar. Gaz üretimi sırasında ise genleşme ve Joule-Thomson etkisi nedeniyle sıcaklık hızla düşer. Tuz kayasının ısıl iletkenliği düşük olduğu için bu sıcaklık değişimleri mağara duvarına yakın dar bir bölgede yoğunlaşır ve burada tekrarlayan ısıl gerilmeler oluşur.

Araştırmacılar bu problemi açıklamak için klasik Burgers sürünme modelini geliştirmiştir. Yeni model, sıcaklıkla elastik modülün zayıflamasını, Arrhenius yasasıyla viskozite değişimini, gerilme artışının sürünmeyi hızlandırmasını ve mikro çatlakların birikerek hasarı büyütmesini birlikte ele alır. Model; laboratuvar ölçekli tuz mağarası deneyleri, akustik emisyon ölçümleri, gaz sızıntısı takibi ve 30 yıllık sayısal simülasyonlarla doğrulanmaya çalışılmıştır.

Çalışmanın ana sonucu oldukça nettir: Yılda 1 enjeksiyon-üretim döngüsü, incelenen JK8-6 Jintan tuz mağarası için 30 yıl sonunda güvenli görünmektedir. Yılda 2 döngüde yerel yapısal kararsızlık riski belirginleşmektedir. Yılda 3 döngüde ise kümülatif hasar faktörü geniş bölgelerde 0,5 mühendislik eşik değerini aşmakta ve makroskobik çatlak ağı, sızdırmazlık kaybı ve hidrojen göçü riski ciddi biçimde artmaktadır.

Bu bulgu, hidrojen depolama sistemlerinin yalnızca hacim, basınç ve ekonomik kapasite bakımından değil, işletme frekansı bakımından da tasarlanması gerektiğini gösterir. Aynı mağara, yılda bir çevrimle güvenli çalışabilirken, yılda üç çevrimle uzun vadede riskli hale gelebilir. Bu nedenle tuz mağarası hidrojen depolaması, yenilenebilir enerji şebekelerine bağlanırken “daha sık çalıştırmak daha iyi esneklik sağlar” mantığıyla değil, kayaç yorulması, termal şok, sürünme deformasyonu ve sızdırmazlık kaybı birlikte değerlendirilerek yönetilmelidir.

Çalışmanın temel problemi nedir?

Bu araştırmanın merkezinde, enerji dönüşümünün çok kritik ama çoğu zaman gözden kaçan bir sorunu vardır: Büyük miktarda yeşil hidrojen üretmek kadar, bu hidrojeni uzun süre güvenli şekilde depolamak da zordur. Rüzgâr ve güneş enerjisi elektrik üretiminde hızla büyümektedir; fakat bu kaynaklar doğaları gereği dalgalıdır. Güneş geceleri üretmez, bulutlu günlerde azalır; rüzgâr ise her zaman aynı hızda esmez. Enerji sisteminin dengede kalabilmesi için fazla üretimin depolanması ve ihtiyaç anında tekrar kullanılabilmesi gerekir.

Yeşil hidrojen bu noktada önemli bir çözüm olarak görülür. Fazla yenilenebilir elektrik, elektroliz yoluyla hidrojene dönüştürülebilir. Hidrojen daha sonra sanayide, ulaşımda, enerji üretiminde veya sezonluk enerji depolamada kullanılabilir. Fakat hidrojenin büyük hacimlerde depolanması için yüzeyde tank kurmak çoğu zaman yeterli ve ekonomik değildir. Bu nedenle yer altı depolama, özellikle de tuz mağaraları, büyük ölçekli hidrojen depolama için güçlü bir jeolojik seçenek olarak öne çıkar.

Tuz mağaraları doğal veya yapay olarak oluşturulan büyük yer altı boşluklarıdır. Tuz kayası çok düşük geçirgenliğe sahiptir; yani gazın kayaç içinden kaçması zordur. Ayrıca tuz kayası zamanla plastik/sürünme davranışı gösterdiği için mikro çatlakları kapatabilme ve kendi sızdırmazlığını kısmen koruyabilme özelliğine sahiptir. Bu özellikler tuz mağaralarını doğal gaz, basınçlı hava ve hidrojen depolama için değerli hale getirir.

Ancak çalışma, yenilenebilir enerjiyle entegre hidrojen depolamada geleneksel işletme koşullarının değiştiğini vurgular. Klasik gaz depolama mağaraları genellikle sezonluk çalışır: yılın belli döneminde doldurulur, belli döneminde boşaltılır. Böyle bir işletme modunda basınç ve sıcaklık değişimleri daha yavaş ve daha az sıklıkta gerçekleşir. Yenilenebilir enerji dalgalanmalarını dengelemek için kullanılan hidrojen depolarında ise daha sık enjeksiyon ve üretim gerekebilir. Bu da mağara çevresindeki tuz kayasını yüksek frekanslı yükleme-boşaltma ve ısıl genişleme-büzülme döngülerine maruz bırakır.

Çalışmanın temel sorusu bu noktada ortaya çıkar: Tuz mağarası hidrojen depolama sistemi, yenilenebilir enerji şebekesinin sık dalgalanmalarını karşılamak için yılda birkaç kez doldurulup boşaltılırsa, mağara çevresindeki kayaç 30 yıl boyunca sızdırmaz ve yapısal olarak güvenli kalabilir mi?

Yüksek frekanslı işletme neden daha tehlikeli olabilir?

Yüksek frekanslı enjeksiyon-üretim, mağara çevresindeki kayaçta iki temel etkiyi aynı anda büyütür: mekanik yorulma ve termal yorulma.

Mekanik tarafta, mağara içindeki gaz basıncı enjeksiyon sırasında artar, üretim sırasında azalır. Basınç arttığında mağara duvarı dışarı doğru zorlanır; basınç düştüğünde çevre kayaç üzerindeki destek azalır. Bu sürekli yükleme-boşaltma, tuz kayasının sürünme deformasyonunu ve gerilme yeniden dağılımını etkiler. Tek bir çevrimde küçük görünen deformasyonlar, onlarca yıl boyunca biriktiğinde mağara hacminde küçülme, çevre kayaçta çatlak gelişimi ve yerel kararsızlık oluşturabilir.

Termal tarafta, gaz davranışı daha karmaşıktır. Enjeksiyon sırasında yüksek basınçlı gaz mağaraya dolarken sıkışma ve ısı birikimi meydana gelir. Bu süreç mağara içi sıcaklığını yükseltir. Üretim sırasında gaz genişler ve Joule-Thomson etkisiyle soğuma ortaya çıkar. Bu kez mağara içi sıcaklığı hızla düşer. Yani mağara duvarı bir çevrimde önce ısınır, sonra soğur.

Tuz kayasının ısıl iletkenliği düşük olduğu için bu sıcaklık değişimi mağaradan uzak bölgelere hızla yayılamaz. Isı alışverişi büyük ölçüde mağara duvarına yakın ince bir termal etkilenmiş bölgede yoğunlaşır. Bu dar bölgede kayaç sürekli ısınıp soğuduğunda, tekrarlayan ısıl genleşme ve büzülme gerilmeleri doğar. Bu gerilmeler mikro çatlakları başlatabilir, var olan çatlakları büyütebilir ve zamanla sızdırmazlık bariyerini zayıflatabilir.

Bu yüzden yüksek frekanslı işletme, yalnızca “aynı işlemi daha sık yapmak” değildir. Her çevrim, mağara duvarında küçük bir mekanik ve termal hasar izi bırakabilir. Çevrim sayısı arttıkça bu izler birleşebilir ve yapısal güvenliği belirleyen kritik hasar seviyelerine ulaşabilir.

Kavramsal enerji sistemi şeması ne anlatıyor?

Çalışmanın başlangıcındaki kavramsal enerji sistemi şeması, yenilenebilir elektrikten yeşil hidrojene ve oradan yer altı depolamaya uzanan zinciri göstermektedir. Rüzgâr ve güneş gibi kesintili kaynaklardan gelen elektrik, güç şebekesiyle ilişkilendirilir. Fazla elektrik, elektroliz sistemlerinde hidrojen üretmek için kullanılır. Üretilen hidrojen büyük hacimli ve uzun süreli depolama amacıyla yer altı tuz mağaralarına gönderilir. Daha sonra bu hidrojen ihtiyaç halinde enerji sistemine geri dönebilir.

Bu şema, çalışmanın neden yalnızca kaya mekaniği değil, enerji sistemi güvenliğiyle de ilgili olduğunu gösterir. Eğer hidrojen depolama mağarası güvenli çalışmazsa, yenilenebilir enerji entegrasyonu için gerekli esneklik de zayıflar. Yani yer altındaki bir kayaç hasarı, yüzeydeki enerji şebekesinin güvenilirliğini dolaylı olarak etkileyebilir.

Periyodik enjeksiyon-üretim mekanizması nasıl açıklanıyor?

Çalışmada verilen işletme mekanizması şeması, gaz enjeksiyonu ve üretimi sırasında mağara duvarında oluşan hasar yollarını özetler. Enjeksiyon aşamasında gaz basıncı ve sıcaklık yükselir. Üretim aşamasında basınç azalır, gaz genişler ve mağara içi sıcaklık düşer. Bu çevrim, mağara duvarının yakınında tekrar eden termal gerilmeler ve mekanik yük değişimleri oluşturur.

Şemada özellikle üç risk öne çıkar:

• Hacim küçülmesi: Tuz kayasının sürünme davranışı nedeniyle mağara hacmi zamanla azalabilir.
• Çekme ve kesme hasarı: Termal şok ve gerilme yeniden dağılımı yerel çekme çatlakları veya kesme hasarı oluşturabilir.
• Sızdırmazlık kaybı: Mikro çatlaklar birleşerek makroskobik çatlak ağına dönüşürse hidrojen göçü ve gaz kaçışı riski artabilir.

Bu noktada çalışmanın önemli bir mesajı vardır: Tuz mağarasında güvenlik yalnızca mağaranın çöküp çökmediğiyle ölçülemez. Mağara geometrik olarak ayakta kalıyor gibi görünse bile, çevre kayaçtaki mikro çatlak ağı sızdırmazlığı bozarsa depolama güvenliği tehlikeye girebilir. Hidrojen molekülünün küçük olması ve sızdırmazlık gereksiniminin yüksek olması, bu riski daha hassas hale getirir.

Klasik Burgers modeli nedir?

Tuz kayası uzun süre yük altında kaldığında anında kırılmak yerine zamanla şekil değiştirir. Bu davranışa sürünme denir. Yer altı tuz mağaralarında sürünme çok önemlidir; çünkü mağara çevresindeki tuz kayası yıllar boyunca yavaşça hareket eder, mağara hacmini küçültür ve gerilme alanını yeniden düzenler.

Klasik Burgers modeli, sürünme davranışını açıklamak için kullanılan viskoelastik bir modeldir. Model, bir Maxwell elemanı ile bir Kelvin elemanının seri bağlanmasıyla oluşturulur. Maxwell elemanı, anlık elastik tepki ve uzun süreli viskoz akışı temsil eder. Kelvin elemanı ise gecikmeli elastik deformasyonu açıklar.

Sabit gerilme altındaki klasik Burgers sürünme bağıntısı şöyle yazılabilir:

\[ \varepsilon(t)= \frac{\sigma}{E_1} + \frac{\sigma}{\eta_1}t + \frac{\sigma}{E_2} \left( 1-e^{-\frac{E_2}{\eta_2}t} \right) \]

Bu formülde \(\varepsilon(t)\), zamana bağlı sürünme şekil değiştirmesidir. \(\sigma\), uygulanan gerilmedir. \(E_1\), Maxwell gövdesinin anlık elastik modülüdür ve yükleme anındaki geri kazanılabilir elastik deformasyonu belirler. \(\eta_1\), Maxwell viskozite katsayısıdır ve uzun vadeli kalıcı viskoz akışı kontrol eder. \(E_2\), Kelvin gövdesinin gecikmeli elastik modülüdür. \(\eta_2\), Kelvin gövdesinin viskozite katsayısıdır ve gecikmeli deformasyonun hızını belirler.

Formüldeki ilk terim anlık elastik deformasyonu, ikinci terim zamanla biriken kalıcı viskoz deformasyonu, üçüncü terim ise zamanla yavaşlayan gecikmeli elastik deformasyonu temsil eder. Bu nedenle Burgers modeli, tuz kayasının başlangıç deformasyonu, yavaşlayan sürünme ve kararlı sürünme davranışını açıklamak için kullanışlıdır.

Klasik Burgers modeli neden yeterli görülmemiş?

Klasik Burgers modelinin önemli bir sınırlılığı vardır: Model parametreleri çoğu zaman sabit kabul edilir. Oysa yer altı hidrojen depolamasında mağara çevresindeki sıcaklık, basınç ve gerilme durumu sürekli değişir. Enjeksiyon sırasında sıcaklık ve basınç artar; üretim sırasında düşer. Bu değişimler tuz kayasının elastik modülünü, viskozitesini ve hasar gelişimini etkiler.

Başka bir deyişle, kayaç her çevrimde aynı malzeme gibi davranmaz. Sıcaklık arttıkça elastik modül zayıflayabilir. Gerilme yükseldikçe viskoz akış hızlanabilir. Mikro çatlaklar biriktikçe efektif taşıyıcı alan azalır ve aynı dış gerilme daha yüksek iç gerilme etkisi yaratır. Bu nedenle klasik model, yüksek frekanslı hidrojen depolama işletmelerindeki termo-mekanik yorulmayı tam yakalamakta yetersiz kalabilir.

Bu çalışma tam olarak bu sınırlılığı gidermek için Burgers modelini sıcaklık, gerilme ve hasar evrimiyle birlikte genişletmektedir.

Sıcaklık elastik modülü nasıl etkiliyor?

Çalışmada tuz kayasının elastik modülünün sıcaklıkla doğrusal değil, üstel biçimde azaldığı kabul edilmiştir. Bu davranış şu MathJax formülüyle ifade edilir:

\[ E_i(T)=E_{i0}\exp(-\beta_E \Delta T) \]

\[ \Delta T=T-T_0 \]

Burada \(E_i(T)\), sıcaklık \(T\) altındaki elastik modülü gösterir. \(E_{i0}\), referans sıcaklıktaki elastik modüldür. \(T_0\), referans sıcaklıktır. \(\Delta T\), mevcut sıcaklık ile referans sıcaklık arasındaki farktır. \(\beta_E\), elastik modülün sıcaklığa duyarlılığını gösteren zayıflama katsayısıdır.

Bu formülün anlamı şudur: Sıcaklık arttıkça tuz kayasının elastik direnci azalır. Yani kayaç aynı gerilme altında daha fazla şekil değiştirmeye yatkın hale gelir. Gaz enjeksiyonu sırasında mağara duvarı ısındığında, bu sıcaklık zayıflaması sürünme deformasyonunu hızlandırabilir.

Viskozite neden Arrhenius yasasıyla değiştiriliyor?

Viskozite, malzemenin uzun vadeli akmaya veya sürünmeye karşı direncidir. Tuz kayası sıcaklık yükseldikçe daha kolay sürünebilir. Bu davranış, birçok termal aktivasyonlu süreçte olduğu gibi Arrhenius tipi bir ilişkiyle açıklanabilir.

Çalışmada viskozite katsayısı sıcaklık ve gerilme etkisiyle şu şekilde değiştirilmiştir:

\[ \eta_j(T,\sigma)= \frac{ \eta_{j0} \exp\left(\frac{Q_j}{R(T+273)}\right) }{ 1+\alpha_\sigma \sigma } \]

Burada \(\eta_j(T,\sigma)\), sıcaklık ve gerilme altında düzeltilmiş viskozite katsayısıdır. \(\eta_{j0}\), referans durumdaki viskozite katsayısıdır. \(Q_j\), viskoz akışın aktivasyon enerjisidir. \(R\), ideal gaz sabitidir. \(T+273\), sıcaklığın Kelvin ölçeğine çevrilmiş halidir. \(\alpha_\sigma\), gerilme pekiştirme veya gerilme duyarlılık katsayısıdır. \(\sigma\), uygulanan gerilme düzeyini ifade eder.

Formül iki temel etkiyi aynı anda taşır. Sıcaklık yükseldiğinde viskoz akış hızlanır. Gerilme arttığında paydadaki \(1+\alpha_\sigma\sigma\) terimi büyür ve etkili viskozite azalır. Viskozite azaldığında kayaç daha kolay sürünür. Böylece formül, “daha yüksek sıcaklık ve daha yüksek gerilme, daha hızlı sürünme deformasyonu üretir” mesajını matematiksel olarak ifade eder.

Hasar değişkeni D ne anlama gelir?

Çalışmanın en kritik taraflarından biri, mikro çatlakların birikimini temsil eden sürekli hasar değişkenidir. Bu değişken \(D\) ile gösterilir ve genel olarak \(0 \leq D \leq 1\) aralığında düşünülür. \(D=0\), malzemenin hasarsız veya sağlam olduğunu; \(D\) büyüdükçe mikro çatlakların arttığını; \(D\) çok yüksek değerlere ulaştığında malzemenin taşıma ve sızdırmazlık yeteneğini kaybettiğini ifade eder.

Sürekli hasar mekaniğine göre mikro çatlaklar arttıkça malzemenin etkin taşıyıcı alanı azalır. Bu durumda aynı dış gerilme, hasarlı malzeme içinde daha yüksek etkin gerilme gibi davranır:

\[ \tilde{\sigma}=\frac{\sigma}{1-D} \]

Burada \(\tilde{\sigma}\), hasar nedeniyle büyütülmüş etkin gerilmedir. \(\sigma\), dışarıdan uygulanan makroskobik gerilmedir. \(D\), hasar değişkenidir. \(D\) arttıkça \(1-D\) küçülür ve etkin gerilme büyür. Bu, hasarın kendi kendini hızlandıran yapısını açıklar: Hasar arttıkça malzeme daha fazla zorlanır; daha fazla zorlandıkça hasar daha hızlı büyür.

Hasar evrimi nasıl modelleniyor?

Çalışmada hasarın zamanla büyümesi için sıcaklık, yükleme frekansı, eşdeğer gerilme ve kendi kendini hızlandırma etkisini içeren bir denklem kullanılmıştır:

\[ \frac{dD}{dt} = C f \left[ \frac{\sigma_{eq}/\sigma_c}{1-D} \right]^m \exp\left[ -\frac{Q_D}{R(T+273)} \right] \]

Bu formülde \(\frac{dD}{dt}\), hasarın zamanla artış hızıdır. \(C\), deneysel olarak kalibre edilen hasar katsayısıdır. \(f\), alternatif yükleme frekansıdır; yani enjeksiyon-üretim çevrimlerinin sıklığını temsil eder. \(\sigma_{eq}\), üç boyutlu gerilme durumunu tek bir eşdeğer gerilme büyüklüğüne çeviren değerdir. \(\sigma_c\), karakteristik dayanımdır. \(m\), hasar evriminin doğrusal olmayan duyarlılığını belirleyen üsteldir. \(Q_D\), hasar aktivasyon enerjisidir. \(R\), ideal gaz sabitidir. \(T\), sıcaklıktır.

Bu denklem çalışmanın kalbidir. Çünkü üç önemli fikri aynı anda taşır:

• İşletme frekansı \(f\) arttıkça hasar daha hızlı birikir.
• Eşdeğer gerilme \(\sigma_{eq}\) arttıkça mikro çatlak gelişimi hızlanır.
• Hasar \(D\) büyüdükçe \(1-D\) küçülür ve sonraki hasar daha da hızlanır.

Bu nedenle yüksek frekanslı hidrojen depolamada risk doğrusal değildir. Yılda 3 çevrim, yılda 1 çevrimin yalnızca üç katı etki üretmeyebilir; hasar kendi kendini hızlandırdığı için uzun vadeli sonuç çok daha ağır olabilir.

Üç boyutlu gerilme durumu nasıl hesaba katılıyor?

Yer altındaki tuz kayası tek eksenli basit bir gerilme altında değildir. Üzerinde örtü yükü, yatay jeolojik gerilmeler, mağara içi gaz basıncı, sıcaklık kaynaklı gerilmeler ve tabaka etkileri birlikte çalışır. Bu nedenle model, tek eksenli gerilme yerine Mohr-Coulomb yaklaşımına dayalı bir eşdeğer gerilme kullanır:

\[ \sigma_{eq} = \frac{1}{3}I_1\sin\varphi + \sqrt{J_2} \left( \cos\theta - \frac{\sin\theta\sin\varphi}{\sqrt{3}} \right) - c\cos\varphi \]

Burada \(\sigma_{eq}\), hasar ve akma değerlendirmesinde kullanılan eşdeğer gerilmedir. \(I_1\), gerilme tensörünün birinci değişmezidir ve ortalama/normal gerilme durumuyla ilişkilidir. \(J_2\), deviatorik gerilme değişmezidir ve şekil değiştirici/kesme etkisini temsil eder. \(\theta\), Lode açısıdır ve üç boyutlu gerilme yolunun karakterini ifade eder. \(\varphi\), içsel sürtünme açısıdır. \(c\), kohezyondur.

Bu formülün pratik anlamı şudur: Mağara çevresindeki kayaç yalnızca basınca göre değil, üç boyutlu gerilme birleşimine göre değerlendirilir. Bu daha gerçekçi bir yaklaşımdır; çünkü tuz mağaralarında hasar çoğu zaman tek bir yöndeki yükten değil, basınç, kesme, çekme ve termal gerilme etkilerinin birleşiminden doğar.

Deney düzeneği nasıl kurulmuş?

Modelin yalnızca teorik kalmaması için laboratuvar ölçekli fiziksel simülasyon deneyleri yapılmıştır. Deneylerde Pakistan’dan alınan doğal tuz kayası örneği kullanılmıştır. Numunenin çapı 200 mm, yüksekliği 400 mm olarak verilmiştir. Numunenin içinde küçük ölçekli bir tuz mağarası oluşturulmuştur. Oluşturulan mağaranın maksimum çapı 62,58 mm, yüksekliği 75,12 mm, hacmi ise 158,64 mL olarak belirtilmiştir.

Bu deney düzeneği, gerçek tuz mağarası işletmesini laboratuvar ölçeğinde taklit etmeyi amaçlar. Numuneye dışarıdan eksenel ve çevresel basınç uygulanır. Mağaranın içine yüksek basınçlı gaz enjekte edilir. Gaz basıncı döngüsel olarak artırılıp azaltılır. Böylece gerçek depolama mağarasında meydana gelen enjeksiyon-üretim çevrimleri küçük ölçekte simüle edilir.

Deney boyunca eksenel deformasyon, radyal deformasyon ve gaz sızıntısı ölçülmüştür. Ayrıca mikro çatlak gelişimini takip etmek için akustik emisyon sistemi kullanılmıştır. Akustik emisyon, kaya içinde çatlak oluşumu veya mikro kayma gibi olaylar sırasında açığa çıkan elastik dalgaları algılar. Bu nedenle laboratuvar ölçeğinde hasarın ne zaman ve hangi hızla geliştiğini anlamak için güçlü bir araçtır.

Deney koşulları nelerdir?

Çalışmada üç ana deney koşulu oluşturulmuştur:

Her deney 5 çevrimden oluşmuştur. Her çevrim 12 saat sürmüş; düşük basınç, yüksek basınç ve basınç ayarlama aşamalarını içermiştir. Gaz olarak helyum kullanılmıştır. Maksimum gaz basıncı 17 MPa, minimum gaz basıncı 7 MPa’dır. Deneylerden sonra önceki hasarın sonraki deneyleri etkilememesi için mağaraya doygun tuzlu su enjekte edilmiş ve 7 MPa basınç altında 3 gün bekletilerek tuz kayasının kendi kendini iyileştirme etkisinden yararlanılmıştır.

Deney sonuçları modelin hangi yönlerini destekliyor?

Deney sonuçları üç ana açıdan değerlendirilmiştir: makroskobik sürünme deformasyonu, akustik emisyonla hasar gelişimi ve gaz sızdırmazlık performansı.

İlk olarak, radyal sürünme deformasyonu belirgin farklılık göstermiştir. Temel koşulda 5 çevrim sonunda radyal şekil değiştirme %0,0159 seviyesine çıkmıştır. Yüksek sıcaklık koşulunda bu değer %0,0245 olmuştur. Yüksek deviatorik gerilme koşulunda ise %0,0388 seviyesine ulaşmıştır.

Bu sonuç şunu gösterir: Sıcaklık artışı tuz kayasının sürünmesini hızlandırır; fakat yüksek deviatorik gerilme daha güçlü bir deformasyon hızlandırıcı etkiye sahiptir. Bu bulgu, modelde hem sıcaklık yumuşaması hem de eşdeğer gerilme sürüş teriminin bulunmasını destekler.

İkinci olarak, akustik emisyonlardan çıkarılan kümülatif hasar değişkeni doğrusal değil, S biçimli bir gelişim göstermiştir. Başta hızlı bir artış, sonra daha kararlı bir gelişim, daha sonra ise hızlanan hasar evrimi gözlenmiştir. Özellikle yüksek deviatorik gerilme koşulunda 5. çevrimde hasar değişkeni D = 0,076 seviyesine ulaşmış ve kararsız yayılım eğilimi göstermiştir. Bu davranış, hasar denklemindeki \(1-D\) paydasının temsil ettiği kendi kendini hızlandırma mekanizmasıyla uyumludur.

Üçüncü olarak, gaz sızıntısı hasar değişkeniyle paralel gelişmiştir. Temel koşulda ve yüksek sıcaklık koşulunda sızıntı artışı daha sınırlı ve daha düzenlidir. Yüksek deviatorik gerilme koşulunda ise ilk çevrimlerde yavaş birikimden sonra 5. çevrimde sızıntı oranı keskin biçimde artarak %0,105 seviyesine ulaşmıştır. Bu, mikro çatlakların birleşerek sızdırmazlık bariyerini bozma eğilimini gösterir.

Bu üç kanıt birlikte değerlendirildiğinde, modelin yalnızca deformasyonu değil, hasar birikimi ve sızdırmazlık kaybı yönünü de temsil edebildiği görülmektedir. Kalibrasyon sonucunda modelin genel uyum katsayısının \(R^2 > 0,92\) seviyesine ulaşması, deney verileriyle güçlü bir uyum sağlandığını göstermektedir.

Kalibre edilen model parametreleri ne anlama geliyor?

Çalışmada termo-mekanik hasarlı Burgers modeli için temel parametreler deney verileriyle kalibre edilmiştir. Öne çıkan bazı değerler şunlardır:

Bu parametreler, modelin yalnızca teorik bir denklem olmadığını, deney verisiyle sayısal olarak çalıştırılabilir hale getirildiğini gösterir. Özellikle \(m=3,5\) değeri, hasar evriminin doğrusal değil, gerilme ve hasar durumuna güçlü biçimde duyarlı olduğunu ifade eder.

Beş boyutlu stabilite değerlendirme sistemi nedir?

Çalışmanın önemli katkılarından biri, tuz mağarası güvenliğini tek bir ölçüte indirgememesidir. Bir mağaranın güvenli olup olmadığını yalnızca deplasmanla, yalnızca gerilmeyle veya yalnızca hasarla değerlendirmek yetersiz olabilir. Bu nedenle araştırma beş boyutlu bir değerlendirme sistemi kurmuştur:

• Toplam deplasman,
• Çekme riski indeksi,
• Dinamik dayanım rezerv faktörü,
• Eşdeğer viskoplastik şekil değiştirme,
• Kümülatif hasar faktörü.

Bu yapı, makroskobik geometri değişimi, çekme çatlağı riski, kesme dayanımı, sürünme akışı ve mikro çatlak birikimini aynı anda değerlendirir. Böylece mağara yalnızca “çöktü mü?” sorusuyla değil, “sızdırmazlığı ve uzun vadeli yapısal bütünlüğü korunuyor mu?” sorusuyla incelenir.

Toplam deplasman nasıl hesaplanıyor?

Toplam deplasman, kayaç içindeki bir noktanın üç uzaysal yöndeki hareketinin birleşik büyüklüğüdür:

\[ U_{\mathrm{total}}=\sqrt{u^2+v^2+w^2} \]

Burada \(u\), \(v\) ve \(w\), sırasıyla \(X\), \(Y\) ve \(Z\) yönlerindeki deplasman bileşenleridir. \(U_{\mathrm{total}}\), bu üç bileşenin toplam geometrik etkisini verir.

Makroskobik güvenlik açısından deplasmanın mağara maksimum çapının %10’unu aşmaması ve hacim küçülmesinin başlangıç hacminin %30’unu geçmemesi önerilir. Ancak çalışma, yalnızca mutlak değere değil, zaman içindeki değişim hızına da dikkat eder. Deplasman hızı giderek artıyorsa, kayaç hızlanan sürünme evresine girmiş olabilir ve bu, kararsızlığın yaklaştığını gösterir.

Çekme riski indeksi neyi ölçer?

Tuz kayası basınca karşı görece dayanıklı, çekmeye karşı ise zayıf bir kayaçtır. Çekme dayanımı genellikle 1-2 MPa gibi düşük seviyelerdedir. Yüksek frekanslı üretim sırasında gaz basıncı ve sıcaklık düştüğünde, mağara duvarında yerel çekme gerilmeleri oluşabilir. Bu nedenle çekme hasarı sızdırmazlık için kritik bir risktir.

Çekme riski indeksi şu şekilde tanımlanır:

\[ Risk_{\mathrm{Tensile}}=\sigma_1-\sigma_t(T) \]

Burada \(Risk_{\mathrm{Tensile}}\), çekme çatlağı riskini gösterir. \(\sigma_1\), birinci asal gerilmedir ve çekme pozitif kabul edilir. \(\sigma_t(T)\), sıcaklık etkisi dikkate alınmış çekme dayanımıdır.

Eğer \(Risk_{\mathrm{Tensile}}<0\) ise kayaç güvenli kabul edilir; çünkü mevcut çekme gerilmesi çekme dayanımının altındadır veya kayaç basınç altındadır. Eğer \(Risk_{\mathrm{Tensile}}>0\) ise yerel çekme çatlağı riski vardır. Tuz mağaralarında bu risk yalnızca mekanik kırılma anlamına gelmez; çatlaklar birleşirse gaz kaçışı ve hidrojen göçü için yol oluşturabilir.

Dinamik dayanım rezerv faktörü neyi gösterir?

Tuz kayası derin yer altında üç eksenli basınç altında bulunur. Bu nedenle kesme dayanımı çevresel basınca ve içsel sürtünme davranışına bağlıdır. Çalışmada dinamik dayanım rezerv faktörü, kaya kütlesinin mevcut gerilme durumunda kesme hasarına karşı ne kadar güvenlik payı taşıdığını ölçmek için kullanılır:

\[ SF= \frac{ \sigma_c(T)+ \sigma_3 \tan^2 \left( 45^\circ+\frac{\varphi(T)}{2} \right) }{ \sigma_1 } \]

Burada \(SF\), dinamik dayanım rezerv faktörüdür. \(\sigma_1\) ve \(\sigma_3\), maksimum ve minimum asal gerilmelerdir. \(\sigma_c(T)\), sıcaklık etkisi altındaki tek eksenli basınç dayanımıdır. \(\varphi(T)\), sıcaklığa bağlı içsel sürtünme açısıdır.

Değerlendirme aralıkları şöyledir:

• \(SF>1,5\): uzun vadeli güvenli bölge,
• \(1,0 \leq SF \leq 1,5\): kritik güvenlik bölgesi,
• \(SF<1,0\): plastik akma veya kesme hasarı riski.

Bu indeks, mağara çevresindeki kayaçta kesme hasarına ne kadar yaklaşıldığını gösterir. Özellikle yılda 2 ve 3 çevrim senaryolarında kritik bölgelerin genişlemesi, uzun vadeli işletme açısından önemli bir uyarıdır.

Eşdeğer viskoplastik şekil değiştirme nasıl tanımlanır?

Tuz kayası üç boyutlu karmaşık bir gerilme alanında sürünür. Bu nedenle yalnızca tek yöndeki deformasyonu izlemek yeterli değildir. Çalışmada Von Mises benzeri eşdeğer şekil değiştirme kullanılmıştır:

\[ \varepsilon_{vp} = \sqrt{ \frac{2}{3} e^{vp}_{ij}e^{vp}_{ij} } \]

Burada \(\varepsilon_{vp}\), eşdeğer viskoplastik şekil değiştirmedir. \(e^{vp}_{ij}\), ortalama hacimsel şekil değiştirme çıkarıldıktan sonra kalan deviatorik viskoplastik şekil değiştirme tensörüdür:

\[ e^{vp}_{ij} = \varepsilon^{vp}_{ij} - \varepsilon^{vp}_m\delta_{ij} \]

\[ \varepsilon^{vp}_m= \frac{1}{3}\varepsilon^{vp}_{kk} \]

Burada \(\delta_{ij}\), Kronecker deltasıdır. Bu formül, üç boyutlu şekil değiştirmeyi tek bir skaler değere indirger. Böylece laboratuvar deneyleriyle ve mühendislik eşikleriyle daha kolay karşılaştırılabilir.

Çalışmada güvenlik aralıkları şu şekilde değerlendirilir:

• \(\varepsilon_{vp}<3\%\): güvenli veya erken kararlı deformasyon bölgesi,
• \(3\% \leq \varepsilon_{vp} \leq 10\%\): dikkat gerektiren belirgin sürünme bölgesi,
• \(\varepsilon_{vp}>10\%\): hızlanan sürünme veya reolojik hasar riski.

Kümülatif hasar faktörü neden en kritik göstergelerden biri?

Kümülatif hasar faktörü \(D\), mikro çatlakların gelişimini ve kayaç sızdırmazlığının bozulma riskini temsil eder. Çalışmada hasar için şu eşikler kullanılmıştır:

• \(D<0,1\): kayaç büyük ölçüde sağlam kabul edilir.
• \(0,1 \leq D \leq 0,5\): hasar başlama ve gelişme bölgesidir.
• \(D>0,5\): mühendislik açısından sızdırmazlık kaybı ve makro çatlak ağı riski ciddi hale gelir.

Matematiksel olarak \(D\) 1’e kadar çıkabilir; ancak mühendislik uygulamasında \(D>0,5\) çok kritik bir eşiği temsil eder. Çünkü bu seviyede mikro çatlakların birleşerek makroskobik çatlak ağına dönüşme olasılığı yüksektir. Hidrojen depolama açısından bu, yalnızca kayaç hasarı değil, doğrudan gaz göçü ve sızdırmazlık güvenliği problemidir.

Sayısal model hangi gerçek mağaraya uygulanmış?

Çalışmada 30 yıllık sayısal analizler Çin’in Jintan bölgesindeki JK8-6 tuz mağarası verileri kullanılarak yapılmıştır. Model, sonar mağara ölçüm verilerine ve yerel stratigrafik özelliklere dayanan üç boyutlu rafine bir jeomekanik modeldir.

Hesaplama bölgesi 1500 m × 1500 m × 1000 m boyutlarında bir blok olarak kurulmuştur. Tuz tabakası yaklaşık 995-1195 m derinlik aralığında yer alır ve 200 m kalınlığındadır. Tuz tabakasının üstünde ve altında yoğun çamurtaşı katmanları bulunur. İncelenen mağaranın yüksekliği yaklaşık 89,7 m, maksimum çapı ise yaklaşık 69,2 m olarak verilmiştir.

Model boyutlarının mağara boyutunun 10 katından fazla seçilmesi, sınır etkilerinin sonucu bozmasını önlemek içindir. Ayrıca gerçek mağara geometrisindeki aşırı düzensizliklerin sayısal tekillik yaratmaması için mağara sınırları yumuşatılmıştır. Bu, gerçek mühendislik verisi ile hesaplanabilir sayısal model arasında denge kuran bir yaklaşımdır.

Üç işletme senaryosu nasıl tanımlanmış?

30 yıllık simülasyonlarda üç farklı enjeksiyon-üretim frekansı incelenmiştir:

Bu senaryolar, yenilenebilir enerji sistemleriyle birlikte tuz mağarası hidrojen depolamanın nasıl değişebileceğini doğrudan test eder. Geleneksel doğal gaz depolaması için yılda bir çevrim yeterli olabilirken, hidrojenin şebeke dengeleme amacıyla daha sık kullanılması gerekebilir. Çalışma bu işletme sıklığının jeomekanik bedelini sayısallaştırmaya çalışır.

Mağara içi sıcaklık alanı nasıl değişiyor?

Termal analizlerde enjeksiyon ve üretim sırasında mağara içinde güçlü sıcaklık dalgalanmaları oluştuğu görülmektedir. Gaz enjeksiyonu sırasında sıkışma nedeniyle mağara içindeki sıcaklık yükselir ve merkez bölgede sıcaklık yaklaşık 68 °C seviyesine ulaşabilir. Gaz üretimi sırasında ise genleşme ve Joule-Thomson etkisiyle sıcaklık hızla düşer; merkez bölgede sıcaklık 45 °C seviyesinin altına inebilir.

Bu sıcaklık dalgalanması yalnızca gaz içinde kalmaz. Mağara duvarındaki tuz kayasıyla ısı alışverişi oluşur. Fakat tuz kayasının düşük ısıl iletkenliği nedeniyle bu ısı değişimi uzak bölgelere yayılmaz; mağara duvarına yakın ince termal etkilenmiş bölgede yoğunlaşır. Bu dar alanda tekrarlayan sıcak-soğuk döngüsü, çevre kayaçta güçlü alternatif termal gerilmeler üretir.

Bu mekanizma, çalışmanın hasar yorumunda merkezi rol oynar. Yani uzun vadeli hasarın ana tetikleyicilerinden biri yalnızca basınç çevrimi değil, basınç çevrimiyle birlikte gelen tekrarlı sıcaklık şokudur.

Yılda 1 çevrim sonucunda mağara güvenli mi?

Geleneksel sezonluk işletme olan yılda 1 çevrim senaryosunda JK8-6 tuz mağarası 30 yıl sonunda genel olarak güvenli görünmektedir. Deformasyon zamana bağlı olarak kararlı ve yaklaşık doğrusal bir sürünme davranışı sergiler. 30 yıl sonunda maksimum yatay deplasman 2,37 m, maksimum zemin oturması 1,85 m olarak hesaplanmıştır. Hacim küçülme oranı ise %4,02 seviyesindedir.

Bu değerler güvenlik sınırlarının oldukça altındadır. Mağara maksimum çapının %10’u ve hacim küçülmesi için %30 sınırı dikkate alındığında, yılda 1 çevrim işletme mağara geometrisi açısından güvenli kalmaktadır.

Diğer göstergeler de aynı yöndedir. Çekme riski genel olarak sıfır veya negatiftir; yani çevre kayaç büyük ölçüde basınç altında kalır. Dinamik dayanım rezerv faktörü çoğu bölgede \(SF>1,5\) seviyesindedir. Eşdeğer viskoplastik şekil değiştirme uzak alanlarda genellikle %3’ün altındadır. Kümülatif hasar faktörü ise çoğu bölgede \(D<0,1\) seviyesinde kalır ve \(D>0,5\) olan bağlantılı tehlikeli bölgeler oluşmaz.

Bu sonuç, incelenen mağara için geleneksel düşük frekanslı işletmenin uzun vadeli sızdırmazlık ve yapısal bütünlük bakımından güvenli olduğunu gösterir.

Yılda 2 çevrimde ne değişiyor?

Yılda 2 enjeksiyon-üretim çevrimi yapıldığında mağaranın uzun vadeli stabilitesi belirgin biçimde zayıflar. 30 yıl sonunda maksimum zemin oturması 3,36 m, maksimum yatay deplasman 2,65 m seviyesine çıkar. Hacim küçülme oranı ise %14,32 olur.

Bu hacim küçülmesi hâlâ %30 sınırının altındadır; ancak yılda 1 çevrim senaryosundaki %4,02 değerine göre çok büyük bir artış gösterir. Bu, çevrim sayısı arttıkça mağara çevresindeki sürünme alanının hızla sıkıştığını ve etkin depolama hacminin azaldığını gösterir.

Çok boyutlu stabilite göstergeleri daha ciddi uyarılar verir. Çekme çatlağı riski özellikle ara tabaka bölgelerinde artar. Dinamik dayanım rezerv faktöründe \(1,0 \leq SF \leq 1,5\) kritik geçiş alanları daha kalın ve daha sürekli hale gelir. Eşdeğer viskoplastik şekil değiştirme ve kümülatif hasar bölgeleri genişler. Mağara duvarında mikro çatlak gelişim alanı büyür ve \(0,1 \leq D \leq 0,5\) hasar evrim bölgesi daha belirgin hale gelir.

Bu senaryo için çalışmanın yorumu dikkatli ve dengelidir: Kısa süreli işletme kontrol edilebilir olabilir; ancak 30 yıl boyunca bu frekansta çalışmak yerel yapısal kararsızlık riskini artırır. Yani yılda 2 çevrim, hemen felaket anlamına gelmez; fakat uzun vadeli işletme periyodu sıkı biçimde kontrol edilmelidir.

Yılda 3 çevrim neden kritik hale geliyor?

Yılda 3 çevrim senaryosu, yüksek frekanslı işletmenin tehlikeli sınırını göstermektedir. 30 yıl sonunda maksimum zemin oturması 4,54 m, maksimum yatay deplasman 3,58 m seviyesine çıkar. Hacim küçülme oranı %19,39 olur.

İlk bakışta bu hacim küçülmesi hâlâ %30 sınırının altında görünür. Ancak çalışma, yalnızca geometrik küçülmeye bakmanın yetersiz olduğunu açıkça gösterir. Çünkü asıl kritik bozulma hasar, çekme riski, kesme dayanımı ve sızdırmazlık göstergelerinde ortaya çıkar.

Bu senaryoda ara tabaka bölgelerinde çekme yırtılması çok daha şiddetli hale gelir. Dinamik dayanım rezerv faktörü düşer ve \(SF<1,0\) bölgeleri genişler. Bu, kaya kütlesinin bazı alanlarda plastik akma veya kesme hasarı eşiğine geldiğini gösterir. Eşdeğer viskoplastik şekil değiştirme, mağara boynu, alt köşe ve yan duvarlarda %10 eşiğini aşar; bu da yerel reolojik hasar riskini gösterir.

En kritik sonuç ise kümülatif hasar faktöründedir. Geniş kaya kütlesi bölgelerinde \(D>0,5\) eşiği aşılır. Bu, mikro çatlakların geniş ölçekte birleşerek makroskobik çatlak ağı oluşturabileceği ve mağaranın sızdırmazlık performansının ciddi biçimde düşebileceği anlamına gelir. Hidrojen depolama açısından bu sonuç çok önemlidir; çünkü bu durumda risk yalnızca yapısal deformasyon değil, gaz sızıntısı ve hidrojen göçüdür.

Çalışmanın vardığı sonuç nettir: Yılda 3 çevrim, acil veya kısa süreli olağanüstü durumlarda düşünülebilir; fakat 30 yıl boyunca standart işletme modu olarak uygulanmamalıdır. Uzun vadede genel kararsızlık, hasar ve sızdırmazlık kaybı riski çok yüksektir.

Üç senaryonun sonuçları nasıl karşılaştırılmalı?

Bu tablo, işletme frekansının doğrusal olmayan etkisini gösterir. Hacim küçülmesi 1 çevrimden 2 çevrime geçildiğinde çok büyük artar. 3 çevrimde ise geometrik sınırlar hâlâ tamamen aşılmamış olsa bile, hasar ve sızdırmazlık açısından kritik eşikler aşılır. Bu nedenle mağara güvenliğinde yalnızca hacim küçülmesi veya deplasman yetmez; hasar faktörü ve çatlak gelişimi de ana karar ölçütü olmalıdır.

Bu çalışma gündelik yaşam açısından nasıl anlaşılmalı?

Bu çalışma teknik olarak derin tuz mağarası, jeomekanik modelleme ve hidrojen depolama üzerine olsa da, gündelik bir benzetmeyle anlaşılabilir. Bir metal teli sürekli ileri geri bükerseniz, tek seferde kırılmaz. Fakat yeterince çok büküldüğünde yorulur ve kopar. Tuz mağarası da benzer şekilde tek bir enjeksiyon veya üretim çevriminde hemen hasar görmeyebilir. Ancak sık tekrarlanan basınç ve sıcaklık değişimleri, mağara duvarındaki kayaçta yavaş yavaş yorulma hasarı biriktirebilir.

Bir diğer benzetme, sıcak-soğuk şokudur. Cam bir bardağa çok sıcak su koyup hemen ardından soğuk su verirseniz, sıcaklık farkı malzemede gerilme oluşturur. Tuz mağarası duvarı da gaz enjeksiyonunda ısınır, gaz üretiminde soğur. Bu süreç yılda bir kez olduğunda kayaç toparlanabilir veya hasar sınırlı kalabilir. Ancak bu çevrim sıklaşırsa, aynı bölgede sürekli termal gerilme oluşur ve mikro çatlaklar büyüyebilir.

Bu nedenle hidrojen depolama güvenliği, yalnızca mağaranın ne kadar büyük olduğu veya ne kadar basınca dayanabildiğiyle değil, nasıl işletildiğiyle de ilgilidir. Aynı mağara, yavaş ve seyrek işletmede güvenli; hızlı ve sık işletmede riskli olabilir.

Enerji dönüşümü açısından çalışmanın önemi nedir?

Enerji dönüşümü çoğu zaman güneş paneli, rüzgâr türbini, elektrolizör ve hidrojen üretimi üzerinden konuşulur. Fakat bu çalışma, yer altı depolama altyapısının da aynı derecede kritik olduğunu hatırlatır. Yeşil hidrojen ekonomisi yalnızca hidrojen üretmekten ibaret değildir. Üretilen hidrojen güvenli, ekonomik ve uzun süreli biçimde depolanmalıdır.

Tuz mağaraları büyük ölçekli hidrojen depolama için en umut verici seçeneklerden biridir. Ancak yenilenebilir enerji dalgalanmalarını hızlı biçimde dengelemek için bu mağaralar daha sık çalıştırılmak istenirse, kaya mekaniği sınırları enerji sistemi işletme stratejisini belirleyen bir faktör haline gelir. Yani şebeke operatörü “daha sık doldur-boşalt” dediğinde, yer altındaki tuz kayasının yorulma sınırı da hesaba katılmalıdır.

Bu çalışma bu açıdan önemli bir köprü kurar: Enerji sistemi esnekliği ile jeomekanik stabilite arasındaki bağlantıyı sayısallaştırır. Sonuçlar, hidrojen depolama mağaralarının yenilenebilir enerjiye destek verebileceğini; ancak yüksek frekanslı işletmenin sınırsız güvenli kabul edilemeyeceğini gösterir.

Geçmiş, bugün ve gelecek açısından çalışmanın anlamı

Geçmiş açısından bakıldığında, tuz mağaraları uzun yıllardır doğal gaz depolama ve bazı enerji depolama uygulamaları için kullanılmıştır. Bu işletmeler genellikle sezonluk pik dengeleme mantığına dayanmıştır. Böyle bir kullanımda çevrim sayısı az, basınç ve sıcaklık değişimleri daha yavaş ve uzun aralıklıdır. Bu nedenle klasik stabilite modelleri birçok uygulamada yeterli kabul edilebilmiştir.

Bugün açısından bakıldığında, yenilenebilir enerji oranı arttıkça depolama altyapısının çalışma mantığı değişmektedir. Hidrojen depoları yalnızca kış-yaz arası sezonluk depolama için değil, daha sık güç dalgalanmalarını dengelemek için de kullanılmak istenebilir. Bu durum, tuz mağarasının daha fazla çevrime maruz kalması anlamına gelir. Çalışma, bu yeni işletme koşullarının klasik güvenlik varsayımlarını zorladığını göstermektedir.

Gelecek açısından bakıldığında, hidrojen ekonomisinin güvenli büyümesi için her tuz mağarasının kendi jeolojik koşulları, sıcaklık-basınç döngüleri, çevrim frekansı, hasar birikimi ve sızdırmazlık güvenliğiyle birlikte modellenmesi gerekecektir. Sabit bir “tuz mağarası güvenlidir” veya “tuz mağarası risklidir” yargısı yeterli değildir. Güvenlik, mağaranın şekline, derinliğine, tabaka yapısına, işletme basıncına, sıcaklık dalgalanmasına ve çevrim sayısına bağlıdır.

Çalışmanın güçlü yönleri nelerdir?

Çalışmanın en güçlü yönlerinden biri, yüksek frekanslı hidrojen depolama işletmesini klasik kaya mekaniği sınırlarının ötesinde ele almasıdır. Model yalnızca elastik deformasyonu değil; sürünme, sıcaklık yumuşaması, viskozite değişimi, gerilme güçlenmesi ve hasar evrimini birlikte dikkate alır.

İkinci güçlü yön, teorik modelin laboratuvar deneyleriyle doğrulanmaya çalışılmasıdır. Küçük ölçekli tuz mağarası fiziksel simülasyonu, yüksek basınçlı gaz çevrimleri, akustik emisyon izleme ve gaz sızıntısı ölçümleri, modelin yalnızca matematiksel bir öneri olarak kalmasını engeller.

Üçüncü güçlü yön, güvenliği beş boyutlu değerlendirmesidir. Deplasman, çekme riski, dayanım rezervi, viskoplastik şekil değiştirme ve hasar faktörü birlikte kullanıldığı için mağara stabilitesi daha bütüncül okunur. Bu özellikle hidrojen depolamada önemlidir; çünkü sızdırmazlık kaybı, yapısal çökmeden önce kritik hale gelebilir.

Dördüncü güçlü yön, gerçek bir saha mağarası olan JK8-6 Jintan verileriyle 30 yıllık sayısal analiz yapılmasıdır. Böylece çalışma, yalnızca idealize bir model değil, mühendislik açısından anlamlı bir uzun dönem senaryo karşılaştırması sunar.

Beşinci güçlü yön, sonuçların işletme politikası açısından net bir mesaj vermesidir. Yılda 1 çevrim güvenli, yılda 2 çevrim dikkatli kontrol gerektiren, yılda 3 çevrim ise uzun dönem standart işletme için ciddi riskli görünmektedir. Bu tür eşikler, hidrojen depolama planlaması için pratik değer taşır.

Çalışmanın sınırlılıkları nelerdir?

İlk ve en önemli sınırlılık, çalışmanın hakem değerlendirmesinden geçmemiş bir ön baskı olmasıdır. Bu nedenle sonuçlar değerli bir bilimsel araştırma çıktısı olarak görülmeli; fakat nihai hakemli yayın sonucu gibi kesin kabul edilmemelidir.

İkinci sınırlılık, deneylerin laboratuvar ölçeğinde yapılmış olmasıdır. Laboratuvar numunesi 200 mm çapında ve 400 mm yüksekliğindedir. Gerçek tuz mağarası ise onlarca metre boyutundadır ve çok daha karmaşık jeolojik tabakalarla çevrilidir. Laboratuvar sonuçları mekanizmayı anlamak için güçlüdür; ancak saha ölçeğine doğrudan birebir aktarılmamalıdır.

Üçüncü sınırlılık, deneylerde helyum kullanılmasıdır. Hidrojenin moleküler davranışı, sızma eğilimi, kimyasal etkileşimleri ve operasyonel güvenlik koşulları helyumdan farklı olabilir. Sayısal analiz hidrojen depolama bağlamına yöneliktir; ancak gerçek hidrojen saha testleri ayrıca gerekir.

Dördüncü sınırlılık, modelin belirli bir mağara geometrisi ve Jintan JK8-6 sahası için uygulanmış olmasıdır. Farklı tuz tabakaları, farklı derinlikler, farklı tabaka arakatkıları, farklı mağara şekilleri ve farklı basınç aralıkları farklı sonuçlar doğurabilir. Bu nedenle her saha için yerel kalibrasyon yapılmalıdır.

Beşinci sınırlılık, 30 yıllık analizlerin belirli çevrim frekanslarıyla sınırlı olmasıdır. Gerçek işletmede çevrimler düzensiz olabilir; bazı dönemlerde sık, bazı dönemlerde seyrek çalışılabilir. Ayrıca acil şebeke desteği, piyasa koşulları ve hidrojen talebi işletme profilini değiştirebilir. Bu nedenle gelecekte düzensiz ve stokastik işletme senaryolarının da incelenmesi gerekir.

Bu çalışma ne söylüyor, ne söylemiyor?

Bu çalışma, tuz mağaralarında hidrojen ve yer altı enerji depolamanın yüksek potansiyel taşıdığını kabul eder; fakat güvenliğin işletme frekansına çok bağlı olduğunu gösterir. Yılda 1 çevrim gibi geleneksel işletme koşullarında incelenen mağara 30 yıl boyunca güvenli ve kontrollü kalmaktadır.

Çalışma, yılda 2 çevrimde risklerin belirgin şekilde arttığını söyler. Bu durumda kısa süreli işletme hâlâ kontrol edilebilir olabilir; ancak 30 yıllık uzun vadeli işletmede yerel kararsızlık riski büyür. Bu nedenle işletme periyodu, hasar izleme ve güvenlik sınırları dikkatle yönetilmelidir.

Çalışma, yılda 3 çevrimde yüksek frekanslı işletmenin uzun vadeli standart mod olarak güvenli olmadığını gösterir. Hasar faktörünün 0,5 eşiğini aşması, mikro çatlakların makro çatlak ağına dönüşme ve sızdırmazlık kaybı riskinin çok ciddi hale geldiğini gösterir.

Buna karşılık çalışma, tüm tuz mağaralarının yılda 3 çevrimde mutlaka başarısız olacağını söylemez. Sonuçlar incelenen model, parametreler, JK8-6 saha geometrisi ve kullanılan işletme koşulları için geçerlidir. Farklı mağaralarda farklı sonuçlar elde edilebilir. Ancak ana bilimsel uyarı güçlüdür: Hidrojen depolama mağaralarında işletme frekansı, yapısal güvenliğin ana tasarım parametrelerinden biri olmalıdır.

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

Bu araştırma, tuz mağaralarında yer altı enerji/hidrojen depolamanın yüksek frekanslı enjeksiyon-üretim operasyonları altındaki termo-mekanik hasar ve uzun vadeli stabilite davranışını inceleyen modelleme, deney ve sayısal simülasyon çalışmasıdır.

Çalışmanın yöntemi üç ana aşamadan oluşur:

• Termo-mekanik hasarlı Burgers sürünme modelinin geliştirilmesi,
• Laboratuvar ölçekli tuz mağarası fiziksel simülasyonu ve akustik emisyon doğrulaması,
• JK8-6 Jintan tuz mağarası için 30 yıllık üç işletme frekansı simülasyonu.

Klasik Burgers sürünme modeli şu şekilde ifade edilir:

\[ \varepsilon(t)= \frac{\sigma}{E_1} + \frac{\sigma}{\eta_1}t + \frac{\sigma}{E_2} \left( 1-e^{-\frac{E_2}{\eta_2}t} \right) \]

Model, sıcaklık etkisini elastik modül zayıflamasıyla dahil eder:

\[ E_i(T)=E_{i0}\exp(-\beta_E \Delta T) \]

\[ \Delta T=T-T_0 \]

Viskozite katsayısı sıcaklık ve gerilme etkisiyle şu şekilde düzeltilir:

\[ \eta_j(T,\sigma)= \frac{ \eta_{j0} \exp\left(\frac{Q_j}{R(T+273)}\right) }{ 1+\alpha_\sigma \sigma } \]

Hasar gelişimi ise şu denklemle temsil edilir:

\[ \frac{dD}{dt} = C f \left[ \frac{\sigma_{eq}/\sigma_c}{1-D} \right]^m \exp\left[ -\frac{Q_D}{R(T+273)} \right] \]

Bu denklemde \(D\) hasar değişkenidir, \(f\) çevrim frekansını, \(\sigma_{eq}\) eşdeğer gerilmeyi, \(\sigma_c\) karakteristik dayanımı, \(C\) ve \(m\) deneysel hasar parametrelerini, \(Q_D\) hasar aktivasyon enerjisini ifade eder. Denklem, yüksek frekanslı çevrimlerin, yüksek gerilmenin ve sıcaklığın hasarı hızlandırdığını; hasar büyüdükçe sonraki hasarın kendi kendini hızlandırdığını gösterir.

Stabilite değerlendirmesinde kullanılan temel göstergeler şunlardır:

\[ U_{\mathrm{total}}=\sqrt{u^2+v^2+w^2} \]

\[ Risk_{\mathrm{Tensile}}=\sigma_1-\sigma_t(T) \]

\[ SF= \frac{ \sigma_c(T)+ \sigma_3 \tan^2 \left( 45^\circ+\frac{\varphi(T)}{2} \right) }{ \sigma_1 } \]

\[ \varepsilon_{vp} = \sqrt{ \frac{2}{3} e^{vp}_{ij}e^{vp}_{ij} } \]

\[ e^{vp}_{ij} = \varepsilon^{vp}_{ij} - \varepsilon^{vp}_m\delta_{ij} \]

\[ \varepsilon^{vp}_m= \frac{1}{3}\varepsilon^{vp}_{kk} \]

Deneysel doğrulamada doğal tuz kayası numunesi kullanılmıştır. Numune çapı 200 mm, yüksekliği 400 mm’dir. Numune içinde maksimum çapı 62,58 mm, yüksekliği 75,12 mm ve hacmi 158,64 mL olan küçük bir tuz mağarası oluşturulmuştur. Deneylerde gaz basıncı 7-17 MPa arasında çevrilmiş ve 5 çevrim uygulanmıştır.

Üç deney koşulu şöyledir:

• T1 temel grup: 23 MPa eksenel/çevresel basınç, 50 °C sıcaklık.
• T2 yüksek deviatorik gerilme grubu: 28,8 MPa eksenel/çevresel basınç, 50 °C sıcaklık.
• T3 yüksek sıcaklık grubu: 23 MPa eksenel/çevresel basınç, 65 °C sıcaklık.

Deney bulguları:

• 5 çevrim sonunda radyal şekil değiştirme T1’de %0,0159, T3’te %0,0245, T2’de %0,0388 olmuştur.
• Yüksek sıcaklık sürünmeyi hızlandırmış, yüksek deviatorik gerilme ise daha güçlü deformasyon artışı oluşturmuştur.
• Akustik emisyon verilerinden çıkarılan hasar değişkeni S biçimli ve doğrusal olmayan bir artış göstermiştir.
• Yüksek deviatorik gerilme koşulunda 5. çevrimde hasar değişkeni \(D=0,076\) seviyesine ulaşmış ve kararsız hasar gelişimi eğilimi göstermiştir.
• Gaz sızıntısı hasar gelişimiyle uyumlu ilerlemiştir.
• Yüksek deviatorik gerilme koşulunda 5. çevrimde kümülatif sızıntı oranı %0,105 seviyesine çıkmıştır.
• Model parametreleri deney verileriyle kalibre edilmiş ve genel uyum \(R^2>0,92\) seviyesinde verilmiştir.

30 yıllık sayısal analizlerde JK8-6 Jintan tuz mağarası kullanılmıştır. Mağara yüksekliği yaklaşık 89,7 m, maksimum çapı yaklaşık 69,2 m’dir. Tuz tabakası 995-1195 m derinlik aralığında ve yaklaşık 200 m kalınlığındadır.

Simülasyon bulguları:

• Enjeksiyon sırasında mağara içi sıcaklık merkezde yaklaşık 68 °C’ye kadar yükselmiştir.
• Üretim sırasında gaz genleşmesi ve Joule-Thomson etkisiyle sıcaklık merkezde 45 °C’nin altına düşmüştür.
• Isı değişimi mağara duvarına yakın ince termal etkilenmiş bölgede yoğunlaşmıştır.
• Bu dar bölgede tekrar eden ısıl gerilme mikro çatlak başlatma ve hasar hızlandırma mekanizmasının ana fiziksel nedeni olarak yorumlanmıştır.

Üç işletme senaryosunun 30 yıllık sonuçları:

• 1 çevrim/yıl: Maksimum yatay deplasman 2,37 m, zemin oturması 1,85 m, hacim küçülmesi %4,02. Çekme riski düşük, \(SF>1,5\) çoğu bölgede korunmuş, \(D<0,1\) baskın. Genel durum güvenli ve kontrollü.
• 2 çevrim/yıl: Zemin oturması 3,36 m, yatay deplasman 2,65 m, hacim küçülmesi %14,32. Kritik dayanım ve hasar bölgeleri genişlemiş. Kısa süreli işletme kontrol edilebilir olsa da 30 yıllık işletme yerel kararsızlık riski taşır.
• 3 çevrim/yıl: Zemin oturması 4,54 m, yatay deplasman 3,58 m, hacim küçülmesi %19,39. \(SF<1,0\) bölgeleri genişlemiş, \(\varepsilon_{vp}>10\%\) alanları oluşmuş ve geniş bölgelerde \(D>0,5\) eşiği aşılmıştır. Uzun vadeli standart işletme için ciddi yapısal kararsızlık ve sızdırmazlık kaybı riski vardır.

Genel olarak çalışma, tuz mağaralarında hidrojen depolamanın güvenliğinin yalnızca depolama basıncı veya mağara geometrisiyle değil, işletme frekansı ve termal-mekanik yorulma birikimiyle de belirlendiğini göstermektedir. Yüksek frekanslı işletme, yenilenebilir enerji entegrasyonu için cazip olsa da, uzun vadeli kayaç hasarı ve hidrojen sızdırmazlığı açısından dikkatle sınırlandırılmalıdır.

Kaynak ve Yöntem Notu

Bu içerik, Jiasong Chen, Zhekang Ding, Tongtao Wang, Youqiang Liao ve Tao He tarafından hazırlanan “Thermo-Mechanical Damage and Stability of Underground Energy Storage Under High-Frequency Operations” başlıklı araştırma metnine dayanarak hazırlanmıştır. Yazar kurumları arasında Chongqing University, PipeChina Energy Storage Technology Co., University of Chinese Academy of Sciences ve Chinese Academy of Sciences bünyesindeki State Key Laboratory of Geomechanics and Geotechnical Engineering yer almaktadır.

Çalışma metninde açıkça “This preprint research paper has not been peer reviewed” ifadesi bulunmaktadır. Ayrıca SSRN elektronik kopya bilgisi yer almaktadır. Bu nedenle bu içerikte çalışma hakem değerlendirmesinden geçmemiş preprint / ön baskı araştırma metni olarak değerlendirilmiştir.

Bu makale hazırlanırken araştırmanın özet, giriş, termo-mekanik hasarlı Burgers modeli, sıcaklık ve gerilme düzeltmeleri, hasar evrim denklemi, Mohr-Coulomb temelli eşdeğer gerilme yaklaşımı, laboratuvar fiziksel simülasyon deneyleri, akustik emisyon doğrulaması, gaz sızıntısı sonuçları, model kalibrasyonu, beş boyutlu stabilite değerlendirme sistemi, JK8-6 Jintan tuz mağarası için 30 yıllık sayısal simülasyonları ve sonuç bölümü temel alınmıştır.

Çalışmada yer almayan saha ölçeğinde kesin hidrojen kaçışı ölçümü, tüm tuz mağaraları için evrensel güvenlik sınırı veya doğrudan ticari işletme standardı gibi iddialar eklenmemiştir. Bulgular, kullanılan model varsayımları, laboratuvar doğrulama koşulları, JK8-6 mağarası geometrisi ve 1-3 çevrim/yıl işletme senaryoları çerçevesinde değerlendirilmelidir.