NCS-v1-2d-base
This is the 2D variant of the NCS-model, a seismic foundation model trained on a large share of full-stack seismic cubes from the Norwegian Continental Shelf (NCS) available through the public DISKOS database. This model has been developed by the Norwegian Computing Center (NR) in collaboration with the industry partners Equinor ASA and AkerBP ASA.
Model Description
NCS-v1-2d-base uses the standard ViT-MAE approach applied to seismic data. It processes individual seismic slices (inline or crossline sections) and uses standard 2D sine-cosine positional encodings.
Usage
NCS-v1-2d-base has been designed to produce features that can be used for fine-tuning on dowsntream tasks such as seismic facies classification, salt body segmentation, geological structure detection (e.g., injectites, faults), content-based seismic image retrieval, horizon and event tracking.
3D context can be incorporated during inference via late fusion: inline and crossline embeddings are computed separately by the NCS encoder and concatenated at each patch location.
How to Use
Loading the Model
Install the NCS package from this repository before running the example below. The standalone
ViTModel loading example only requires transformers.
1from transformers import ViTModel
2
3model = ViTModel.from_pretrained("NorskRegnesentralSTI/NCS-v1-2d-base", add_pooling_layer=False)
Feature Extraction
1import torch
2
3# Input: batch of 3-channel (1 channel repeated) seismic crops, normalized by cube std and clamped to 3σ
4pixel_values = torch.randn(1, 3, 224, 224)
5
6with torch.no_grad():
7 outputs = model(pixel_values=pixel_values)
8
9# CLS token (image-level feature)
10cls_features = outputs.last_hidden_state[:, 0, :] # shape: (B, 768)
11
12# Patch-level features
13patch_features = outputs.last_hidden_state[:, 1:, :] # shape: (B, 196, 768)
Inference on Seismic Volumes
For running inference over full seismic volumes (SEG-Y / SGZ), use the NCS inference pipeline:
1uv run scripts/inference.py \
2 --model-path NorskRegnesentralSTI/NCS-v1-2d-base \
3 --input-path /path/to/volume.segy \
4 --output-path ./features_2d.zarr \
5 --direction dir0 \
6 --densify 1 \
7 --num-overlap-patches 7 \
8 --overlap-filter ramp \
9 --batch-size 32 \
10 --device cuda:0 \
11 --dtype float16
Training Details
Pretraining Data
The model was pretrained on seismic reflection data from the Norwegian Continental Shelf (NCS), sourced from the DISKOS national data repository. The training corpus consists of 829 full-stack time and depth migrated 3D seismic cubes (~27 TB), spanning diverse geological settings, acquisition vintages, and processing generations across the NCS.
Preprocessing
- Seismic amplitudes are standardized per-cube to unit variance.
- Values are clipped at ±3 standard deviations.
- For each training sample, 2D slices are extracted at 4 azimuthal directions (0°, 45°, 90°, 135°) through the same spatial location.
- Diagonal slices (45°, 135°) are center-cropped and resized to correct for the √2 elongation.
- Single-channel slices are passed as separate views to a shared patch projection layer.
Training Procedure
- Pretraining method: Masked Autoencoder (MAE) with 85% masking ratio (applied after concatenating patches across views; per-view mask count is not enforced)
- Initialization: ImageNet MAE ViT weights (RGB projection channels averaged to single-channel input)
- Framework: PyTorch with flash-attention kernels
- Hardware: 16 × NVIDIA GH200 GPUs
- Precision: bfloat16 mixed precision
- Global batch size: 2048
- Learning rate: Cosine schedule, base LR = 1.5 × 10⁻⁴, effective LR = base_lr × batch_size / 256, warmup ratio = 0.05
- Epochs: 100 (~1M samples per epoch)
- Sampling: Density-aware sampling from seismic cubes, biased toward regions with sparser spatial coverage
- Decoder: Lightweight 8-layer MAE decoder
Evaluation Protocol
Representations are evaluated with a frozen backbone using a k-nearest-neighbor (kNN, k=5) classifier on patch-level embeddings. Four interpretation benchmarks were used: salt segmentation, package segmentation, injectite mapping, and flatspot mapping, measured by mean Intersection-over-Union (mIoU). Only 100 labeled points per class are used (or a single labeled line for injectites).
Code
The model code and inference pipeline are available at:
https://github.com/NorskRegnesentral/NCS_modelsCitation
If you use this model, please cite:
1@article{ordonez2025ncsmodel,
2 title={The {NCS}-model: A seismic foundation model trained on the Norwegian repository of public seismic data},
3 author={Ordo{\~n}ez, Alba and Forgaard, Theodor Johannes Line and Wade, David and Bugge, Aina Juell and Nese, H{\aa}kon and Waldeland, Anders Ueland},
4 journal={arXiv preprint arXiv:2603.23211},
5 year={2025}
6}
Acknowledgments
This work is funded by The Research Council of Norway through the SFI Visual Intelligence (Centre for Research-based Innovation), grant no. 309439, and the industry partners Equinor ASA and AkerBP ASA. We also thank Equinor and AkerBP for providing access to the seismic data used in the evaluation.
