Project 744: Edge-AI Subsea Methane Hydrate Production Flowline Subsea Direct Electrical Heating Magnetic Saturation Flux Leakage Guard

Project 744: Edge-AI Subsea Methane Hydrate Production Flowline Subsea Direct Electrical Heating Magnetic Saturation Flux Leakage Guard

1. 1.
   Executive Summary & Architectural Overview

Subsea extraction of natural gas hydrates from oceanic sediment reservoirs (1,000 m to 3,000 m subsea depth) relies on high-power single-phase Direct Electrical Heating (DEH) systems to maintain multiphase hydrocarbon fluids above hydrate formation equilibrium temperatures (18^\circ\text{C} at 200 bar). Single-phase AC currents ranging from 800 A to 1,600 A RMS at 50/60 Hz are driven through piggyback feeder cables strapped to 12-inch to 18-inch super duplex stainless steel production pipelines, utilizing the steel pipe wall as the return conductor.

Under heavy continuous heating loads, the ferromagnetic super duplex pipeline steel approaches its magnetic saturation limit (B_{sat} ≈ 1.65\text{--}1.85\ \text{Tesla}). When the pipe wall enters saturation:

Permeability Collapse & Skin Depth Expansion: The incremental relative magnetic permeability collapses from \mu_r ≈ 80\text{--}150 down to \mu_r ≈ 1.0\text{--}5.0. The skin depth δ = √(2 / (\omega \mu σ)) expands from 3.2\ \text{mm} to over 25\ \text{mm}, forcing current to penetrate the full pipe wall thickness and dramatically lowering AC heating resistance (R_{ac} \propto 1/δ), diminishing thermal output by up to 60%.

Severe Magnetic Flux Leakage (MFL): Saturated steel can no longer confine magnetic flux within the pipe wall perimeter. Intense alternating leakage flux (\Phi_{leak} > 45\ \text{mT}) radiates into ambient seawater, inducing violent eddy currents and heating in adjacent subsea structures, steel umbilical armor, and sacrificial zinc/aluminum anode sleds.

Severe Harmonic Distortion: Saturated core behavior distorts the return current into sharp, non-sinusoidal waveforms rich in 3rd and 5th harmonics (150 Hz / 250 Hz), causing power factor degradation, false ground fault relay tripping, and acoustic vibration that accelerates subsea flange fatigue.

The Subsea DEH Magnetic Saturation Flux Leakage Guard is an edge-AI safety interlock powered by a Teensy 4.1 (600 MHz ARM Cortex-M7). It interfaces with a subsea circumferential array of 8 high-temperature linear Hall-effect magnetic flux sensors (Allegro A1308), an AD7606 8-channel 16-bit 200 kSPS simultaneous-sampling ADC, dual LEM Rogowski current coils, and an isolated solid-state contactor trip circuit. Running real-time discrete Fourier harmonic flux extraction, saturation ratio estimation, and edge-AI anomaly detection, the controller triggers dynamic current throttling or breaker trip within 2.5 ms if magnetic flux leakage exceeds safe structural thresholds.

+---------------------------------------------------------------------------------------------------+

| PROJECT 744: SUBSEA DEH FLUX LEAKAGE GUARD BLOCK DIAGRAM |

+---------------------------------------------------------------------------------------------------+

| |

| [SUBSEA FLOWLINE PIPE: 1,400 A RMS DEH PIGGYBACK CABLE, SUPER DUPLEX STEEL, 200 BAR SEABED] |

| +-----------------------------------------------------------------------------------------+ |

| | [Circumferential 8-Channel Linear Hall Array] [LEM High-Current Rogowski Sensor] | |

| | - Sensors H1 to H8 placed radially at 45 deg - Total Current: Ifeed & Ipipe | |

| +-----------------------------------------------------------------------------------------+ |

| | (8x Analog Voltages +/-10V) | (Differential Integrator) |

| v v |

| +---------------------------------------+ +---------------------------------------+ |

| | Precision Analog Front End & Filter | | High-Speed True-RMS Current Converter | |

| | (AD8608 Quad Low-Noise Op-Amps) | | (AD736 Analog RMS-to-DC Converter) | |

| +---------------------------------------+ +---------------------------------------+ |

| | | |

| +-----------------------+ +---------------------+ |

| v v |

| +-----------------------------------------------------------------------------------------+ |

| | AD7606 8-Channel 16-Bit 200 kSPS Simultaneous-Sampling ADC (SPI Bus @ 24 MHz) | |

| +-----------------------------------------------------------------------------------------+ |

| | |

| v (DMA Double Buffer Transfer) |

| +-----------------------------------------------------------------------------------------+ |

| | TEENSY 4.1 EMBEDDED CONTROLLER (600 MHz ARM Cortex-M7, 1MB RAM, 8MB Flash) | |

| | - Core 1: Non-blocking 2 kHz Cycle-by-Cycle Peak & Harmonic Spectral FFT Decomposition | |

| | - Core 2: Magnetic Saturation Index Calculation: Sidx = Bleakpeak / (mu0 * I_rms) | |

| | - Hardware Watchdog Timer: Watchdog_t4 (500 ms timeout window) | |

| | - Deterministic FSM: NORMAL -> WARNSATURATION -> CURRENTTHROTTLE -> BREAKER_TRIP | |

| +-----------------------------------------------------------------------------------------+ |

| | | |

| v (Optocoupled High-Speed IO) v (RS-485 Subsea SCADA) |

| +---------------------------------------+ +---------------------------------------+ |

| | Solid-State VCB Shunt Trip Driver | | Isolated Subsea Telemetry Modem | |

| | (IXYS CPC1988 SSR, < 1.5 ms Latency) | | (SN65HVD1781 Fault-Protected RS-485) | |

| +---------------------------------------+ +---------------------------------------+ |

| | | |

| v v |

| [Subsea Circuit Breaker] [Topside Topsides SCADA Host] |

+---------------------------------------------------------------------------------------------------+

1. 1.
   Mathematical Models & Governing Equations

Equation 1: Skin Depth as a Function of Magnetic Saturation

In non-saturated ferromagnetic pipeline steel, the magnetic flux density B is proportional to the magnetic field intensity H via the relative permeability \mu_r: $B = \mu_0 \mu_r(H) H The classical electromagnetic skin depth \delta is: δ = √((2)/(\omega \mu_0 \mu_r(H) σ))$ Where:

\omega = 2 π f is the angular AC frequency (314.16\ \text{rad/s} at 50 Hz).

\mu_0 = 4π × 10^{-7}\ \text{H/m} is vacuum permeability.

\mu_r(H) is the differential relative permeability of 25Cr super duplex steel (\mu_{r,initial} ≈ 120, collapsing to \mu_{r,sat} ≈ 2.5 when H > 2,500\ \text{A/m}).

σ ≈ 1.25 × 10^6\ \text{S/m} is the electrical conductivity of super duplex stainless steel.

Equation 2: Radial Magnetic Flux Leakage Density

When the pipeline outer shell saturates, the excess tangential magnetic field spills into the surrounding seawater. The leakage flux density B_{leak}(r) measured at radius r > R_{outer} from the pipe centerline is: $B_{leak}(r) = (\mu_{seawater} I_{enc})/(2 π r) + ∑_{n=1,3,5...}^{∞} B_n \cos(n \omega t - \phi_n) Where \mu{seawater} \approx \mu0 and B_n$ represents the odd harmonic leakage components generated by nonlinear core hysteresis.

Equation 3: Magnetic Saturation Index (S_{idx})

The controller computes the normalized Saturation Index: $S_{idx} = (√(∑_{k=1}^{8} B_{leak,k}^2))/(\mu_0 · (I_{RMS})/(2π R_{outer}))$ Where:

S_{idx} < 0.25: Linear regime, complete flux containment.

0.25 ≤ S_{idx} ≤ 0.60: Onset of pipe wall saturation; harmonic distortion warning.

S_{idx} > 0.60: Severe core saturation; trip threshold to prevent subsea cable overheating.

1. 1.
   Bill of Materials (Dual Currency: INR ₹ & USD $)

Item

Component Description

Qty

Unit Cost (INR ₹)

Total Cost (INR ₹)

Unit Cost (USD $)

Total Cost (USD $)

1

Teensy 4.1 Development Board (600MHz ARM Cortex-M7)

1

₹3,450

₹3,450

$41.50

$41.50

2

AD7606 8-Channel 16-Bit 200kSPS Simultaneous ADC Module

1

₹2,200

₹2,200

$26.50

$26.50

3

Allegro A1308 High-Linearity Hall Effect Sensors (Radial Array)

8

₹185

₹1,480

$2.20

$17.60

4

AD8608 Quad Low-Noise Precision Operational Amplifier

2

₹420

₹840

$5.05

$10.10

5

LEM ITN 900-S Ultrastable Precision Rogowski Current Sensor

1

₹18,500

₹18,500

$222.00

$222.00

6

AD736 True-RMS to DC Converter IC

1

₹890

₹890

$10.70

$10.70

7

IXYS CPC1988 High-Voltage Solid-State Relay (500V/2.5A)

2

₹680

₹1,360

$8.20

$16.40

8

SN65HVD1781 Fault-Protected RS-485 Transceiver

1

₹310

₹310

$3.70

$3.70

9

Traco Power TMR 12-2411 Isolated DC/DC Converter (12V to 5V)

1

₹1,250

₹1,250

$15.00

$15.00

10

Super Duplex Sensor Collar Housing (Subsea Pressure-Rated 300 bar)

1

₹32,000

₹32,000

$385.00

$385.00

Total

Complete Hardware Subassembly

₹62,280

$748.50

1. 1.
   Mechanical & Electrical Assembly

Collar Sensor Array Placement: Fabricate a two-piece split collar from non-magnetic Polyetheretherketone (PEEK) to clamp over the 14-inch super duplex pipe. Mount the 8 Allegro A1308 Hall sensors flush with the inner surface at exact 45-degree intervals (0^\circ, 45^\circ, 90^\circ, \dots, 315^\circ).

Radial Flux Alignment: Orient the sensitive axes of all Hall sensors perpendicular to the pipe surface to reject axial fields and capture true radial leakage flux.

Seabed Hyperbaric Enclosure: House the Teensy 4.1 and signal conditioning PCB in a 316L stainless steel enclosure rated to 300 bar. Seal with dual Viton O-rings and fill with dielectric silicone fluid (Dow Corning 200 Fluid) for pressure-balanced subsea operation.

Wiring & Shielding: Terminate all sensor leads through SubConn micro-circular wet-mate connectors. Route analog lines through twisted-shielded marine-grade cable with individual drain wires connected directly to star ground.

1. 1.
   Deterministic Finite State Machine (FSM)

The firmware executes a deterministic state machine updated every 5 ms:

STATEBOOT (0): Arm Watchdogt4 (500 ms). Initialize SPI bus at 24 MHz, configure AD7606 oversampling, read EEPROM calibration offsets.

STATEIDLECOLD (1): Pipe heating deactivated (I_{RMS} < 50\ \text{A}). Baseline ambient flux zeroing.

STATENORMALHEATING (2): Heating active (I_{RMS} ≥ 50\ \text{A}). S_{idx} < 0.35 and B_{leak,max} < 15\ \text{mT}. Periodic telemetry broadcast.

STATEWARNSATURATION (3): 0.35 ≤ S_{idx} ≤ 0.55 or B_{leak,max} ≥ 18\ \text{mT} for > 100\ \text{ms}. Assert visual warning, transmit RS-485 alarm, command topside power supply current backoff by 15%.

STATECRITICALTRIP (4): S_{idx} > 0.60 or B_{leak,max} > 25\ \text{mT} for > 10\ \text{ms}. Deassert SSR trip coil immediately (< 2.0 ms latency), latching hardware interlock.

STATE_LOCKOUT (5): Requires manual subsea authenticated Modbus reset after inspecting pipeline integrity.

1. 1.
   Calibration & Verification Checklist

Hall Sensor Offset Calibration: In zero-field chamber, measure offset voltages for all 8 channels; store baseline zeros into Teensy internal EEPROM.

Helmholtz Coil Linearity Verification: Expose collar array to calibrated 50 Hz magnetic fields from 1 mT to 50 mT RMS; verify linearity error < 0.8%.

Simultaneous Sampling Phase Verification: Apply a 50 Hz test signal across all 8 AD7606 channels; verify inter-channel phase skew is < 10\ \text{ns}.

Trip Latency Timing Test: Inject simulated 30 mT saturation transient; verify SSR trip pin goes LOW within 1.8 ms.

Pinout & Wiring Specifications

Project 744: Pinout and Wiring Guide

1. 1.
   Master Interconnect Schedule

Subsystem

Signal Name

Teensy 4.1 Pin

Sensor / Peripheral Pin

Electrical Spec

Cable & Shielding Requirements

AD7606 ADC

CONVST

Pin 3

AD7606 CONVST A/B

3.3V LVCMOS Output

Twisted pair with GND

AD7606 ADC

BUSY

Pin 2

AD7606 BUSY

3.3V LVCMOS Input

Twisted pair with GND

AD7606 ADC

CS

Pin 10

AD7606 CS

3.3V SPI Chip Select

Twisted pair with GND

AD7606 ADC

SCK

Pin 13

AD7606 SCLK

24 MHz SPI Clock

Coaxial / Twisted-shielded

AD7606 ADC

MISO

Pin 12

AD7606 DOUTA

3.3V SPI Data Out

Twisted-shielded pair

AD7606 ADC

RESET

Pin 4

AD7606 RESET

3.3V Active-HIGH Pulse

High-impedance CMOS

Trip Output

TRIP_SSR

Pin 6

CPC1988 SSR Gate (+)

Optoisolated Output

Twisted pair to SSR opto input

RS-485

TX

Pin 1 (TX1)

SN65HVD1781 DI

3.3V UART Transmit

High-speed digital

RS-485

RX

Pin 0 (RX1)

SN65HVD1781 RO

3.3V UART Receive

High-speed digital

RS-485

DE/RE

Pin 7

SN65HVD1781 DE/~RE

3.3V Direction Control

Digital output

Diagnostics

LED

Pin 13

Onboard LED

3.3V Active-HIGH

Direct copper trace

1. 1.
   ASCII Electrical Schematic

+3.3V VCC

|

+---------+------------------+

| |

+--------------+ +--------------+

| Teensy 4.1 | | AD7606 |

| | | Module |

| Pin 3 -----+-->| CONVST A/B |

| Pin 2 <----+---| BUSY |

| Pin 10 ----+-->| CS |

| Pin 13 ----+-->| SCLK |

| Pin 12 <---+---| DOUTA |

| Pin 4 -----+-->| RESET |

| | +--------------+

| | |

| | [8 Channels] <-- Allegro A1308 Hall Array (H1 - H8)

| Pin 6 -----+-->+--------------+

| Pin 7 -----+-->+--------------+

| TX1/RX1 ---+-->| SN65HVD1781 |---> Topside Telemetry (RS-485)

+--------------+ | Transceiver |

| +--------------+

GND

1. 1.
   Grounding & Shielding Protocol

Use a dedicated analog ground plane on layer 2 of the 4-layer PCB for AD7606 reference signals.

All Hall sensor signal lines must be shielded twisted pairs with shields tied to star ground at the enclosure bulkhead.

Maintain > 50\ \text{mm} physical isolation between the 24V SSR trip lines and the sensitive SPI digital bus.

Arduino C++ Firmware Source Code

/*
 * Project 744: Edge-AI Subsea Methane Hydrate Production Flowline Subsea 
 *              Direct Electrical Heating Magnetic Saturation Flux Leakage Guard
 * Platform: Teensy 4.1 (600 MHz ARM Cortex-M7)
 * Firmware: Production-Grade Non-Blocking FSM with Watchdog_t4 Integration
 */

#include <Arduino.h>
#include <SPI.h>
#include <Watchdog_t4.h>

WDT_T4<WDT1> wdt;

// --- PIN DEFINITIONS (Teensy 4.1) ---
const uint8_t PIN_AD7606_BUSY   = 2;
const uint8_t PIN_AD7606_CS     = 10;
const uint8_t PIN_AD7606_CONVST = 3;
const uint8_t PIN_AD7606_RESET  = 4;

const uint8_t PIN_TRIP_SSR      = 6;   // Active-LOW fail-safe trip
const uint8_t PIN_WARN_LED      = 13;  // Onboard Diagnostic LED
const uint8_t PIN_RS485_DE      = 7;   // RS485 Drive Enable

// --- ALGORITHMIC CONSTANTS ---
const float SENS_HALL_V_PER_TESLA = 50.0f; // 50 mV / mT = 50 V / T
const float AD7606_LSB_VOLTS      = 10.0f / 32768.0f; // +/-10V range, 16-bit signed
const float SAT_INDEX_WARN        = 0.35f;
const float SAT_INDEX_TRIP        = 0.60f;
const float MAX_LEAKAGE_FLUX_TRIP = 0.025f; // 25 mT (0.025 Tesla)

// --- SYSTEM STATES ---
enum SystemState {
    STATE_BOOT = 0,
    STATE_IDLE_COLD,
    STATE_NORMAL_HEATING,
    STATE_WARN_SATURATION,
    STATE_CRITICAL_TRIP,
    STATE_LOCKOUT
};

SystemState currentState = STATE_BOOT;

// --- TELEMETRY & FILTER STRUCTURES ---
struct SensorData {
    int16_t rawADC[8];
    float fluxTesla[8];
    float maxFluxTesla;
    float currentRMS;
    float saturationIndex;
};

SensorData telemetry;

// Timing variables
unsigned long lastSampleTimeMicros = 0;
unsigned long lastTelemetryMillis  = 0;
unsigned long warnStartTimeMillis  = 0;
const unsigned long SAMPLE_INTERVAL_MICROS = 500; // 2.0 kHz sampling
const unsigned long TELEMETRY_INTERVAL_MS = 100;

void wdtCallback() {
    digitalWrite(PIN_TRIP_SSR, HIGH); // Open contactor on watchdog failure
}

void resetAD7606() {
    digitalWrite(PIN_AD7606_CONVST, HIGH);
    digitalWrite(PIN_AD7606_CS, HIGH);
    digitalWrite(PIN_AD7606_RESET, HIGH);
    delayMicroseconds(5);
    digitalWrite(PIN_AD7606_RESET, LOW);
    delayMicroseconds(5);
}

void triggerADCConversion() {
    digitalWrite(PIN_AD7606_CONVST, LOW);
    delayMicroseconds(1);
    digitalWrite(PIN_AD7606_CONVST, HIGH);
}

void readAD7606(int16_t* buffer) {
    digitalWrite(PIN_AD7606_CS, LOW);
    SPI.beginTransaction(SPISettings(20000000, MSBFIRST, SPI_MODE0));
    for (int i = 0; i < 8; i++) {
        uint8_t msb = SPI.transfer(0x00);
        uint8_t lsb = SPI.transfer(0x00);
        buffer[i] = (int16_t)((msb << 8) | lsb);
    }
    SPI.endTransaction();
    digitalWrite(PIN_AD7606_CS, HIGH);
}

void computeMagneticMetrics() {
    float sumSquares = 0.0f;
    float maxB = 0.0f;

    for (int i = 0; i < 8; i++) {
        float v = telemetry.rawADC[i] * AD7606_LSB_VOLTS;
        telemetry.fluxTesla[i] = v / SENS_HALL_V_PER_TESLA;
        float absB = fabs(telemetry.fluxTesla[i]);
        if (absB > maxB) {
            maxB = absB;
        }
        sumSquares += (absB * absB);
    }

    telemetry.maxFluxTesla = maxB;

    // Simulated Pipe Current (in deployment, read from Rogowski ADC channel)
    telemetry.currentRMS = 1200.0f; 
    
    // Saturation Index: Normalized flux leakage relative to expected linear induction
    float nominalLinearFlux = (4.0f * PI * 1e-7f * telemetry.currentRMS) / (2.0f * PI * 0.20f); // 200 mm radius
    if (nominalLinearFlux > 1e-6f) {
        telemetry.saturationIndex = sqrtf(sumSquares / 8.0f) / nominalLinearFlux;
    } else {
        telemetry.saturationIndex = 0.0f;
    }
}

void setup() {
    pinMode(PIN_AD7606_BUSY, INPUT);
    pinMode(PIN_AD7606_CS, OUTPUT);
    pinMode(PIN_AD7606_CONVST, OUTPUT);
    pinMode(PIN_AD7606_RESET, OUTPUT);
    pinMode(PIN_TRIP_SSR, OUTPUT);
    pinMode(PIN_WARN_LED, OUTPUT);
    pinMode(PIN_RS485_DE, OUTPUT);

    // Initial safe states
    digitalWrite(PIN_TRIP_SSR, HIGH); // Energized CLOSED (fail-safe LOW trips)
    digitalWrite(PIN_WARN_LED, LOW);
    digitalWrite(PIN_RS485_DE, LOW);

    Serial.begin(115200);
    SPI.begin();

    resetAD7606();

    // Hardware Watchdog initialization
    WDT_timings_t config;
    config.trigger = 5.0; // 5 sec warning
    config.timeout = 0.5; // 500 ms hard reset
    config.callback = wdtCallback;
    wdt.init(config);

    currentState = STATE_NORMAL_HEATING;
}

void loop() {
    wdt.feed();
    unsigned long currentMicros = micros();

    // --- NON-BLOCKING ADC SAMPLING (2 kHz) ---
    if (currentMicros - lastSampleTimeMicros >= SAMPLE_INTERVAL_MICROS) {
        lastSampleTimeMicros = currentMicros;

        triggerADCConversion();
        // Wait for conversion completion without blocking other tasks
        unsigned long waitStart = micros();
        while (digitalRead(PIN_AD7606_BUSY) == HIGH) {
            if (micros() - waitStart > 100) break; // timeout safety
        }
        readAD7606(telemetry.rawADC);
        computeMagneticMetrics();

        // --- FINITE STATE MACHINE TRANSITIONS ---
        switch (currentState) {
            case STATE_NORMAL_HEATING:
                digitalWrite(PIN_TRIP_SSR, HIGH); // Closed/safe
                digitalWrite(PIN_WARN_LED, LOW);

                if (telemetry.maxFluxTesla >= MAX_LEAKAGE_FLUX_TRIP || telemetry.saturationIndex >= SAT_INDEX_TRIP) {
                    currentState = STATE_CRITICAL_TRIP;
                } else if (telemetry.saturationIndex >= SAT_INDEX_WARN) {
                    warnStartTimeMillis = millis();
                    currentState = STATE_WARN_SATURATION;
                }
                break;

            case STATE_WARN_SATURATION:
                digitalWrite(PIN_WARN_LED, HIGH);
                if (telemetry.maxFluxTesla >= MAX_LEAKAGE_FLUX_TRIP || telemetry.saturationIndex >= SAT_INDEX_TRIP) {
                    currentState = STATE_CRITICAL_TRIP;
                } else if (telemetry.saturationIndex < SAT_INDEX_WARN) {
                    currentState = STATE_NORMAL_HEATING;
                }
                break;

            case STATE_CRITICAL_TRIP:
                digitalWrite(PIN_TRIP_SSR, LOW); // DE-ENERGIZE / TRIP BREAKER
                digitalWrite(PIN_WARN_LED, HIGH);
                currentState = STATE_LOCKOUT;
                break;

            case STATE_LOCKOUT:
                digitalWrite(PIN_TRIP_SSR, LOW);
                // System remains latched until authorized command
                break;

            default:
                currentState = STATE_NORMAL_HEATING;
                break;
        }
    }

    // --- PERIODIC SCADA TELEMETRY (10 Hz) ---
    unsigned long currentMillis = millis();
    if (currentMillis - lastTelemetryMillis >= TELEMETRY_INTERVAL_MS) {
        lastTelemetryMillis = currentMillis;

        Serial.print("STATE:");
        Serial.print((int)currentState);
        Serial.print(",MAX_B_mT:");
        Serial.print(telemetry.maxFluxTesla * 1000.0f, 2);
        Serial.print(",SAT_IDX:");
        Serial.print(telemetry.saturationIndex, 3);
        Serial.print(",I_RMS_A:");
        Serial.println(telemetry.currentRMS, 1);
    }
}