Analysis and Design of Spare Strategy for Large-Scale Satellite Constellation Using Direct Insertion under (r,q) Policy

Technische Samenvatting: Analyse en Ontwerp van Reserve-Strategie voor Groot-Schaal Satellietconstellaties met Directe Injectie onder (𝑟, 𝑞)-beleid

Probleemstelling
Groot-schaal satellietconstellaties in de Lage Aarde-Orbit (LEO) vereisen robuust reservebeheerbeleid om prestaties te handhaven ondanks satellietuitval. Hoewel reparatie in de baan vaak onbetaalbaar duur is, wordt het lanceren van vervangende satellieten steeds kosteneffectiever door dalende lanceerkosten en efficiënties van massaproductie. Er bestaan twee primaire strategieën voor aanvulling: indirect (reserves in groepsverband in parkeerbanen brengen voor latere overdracht op basis van drift) en direct (onmiddellijke injectie in constellatievlakken). Dit artikel richt zich op de directe reserve-strategie, waarbij de behoefte wordt aangepakt aan een rigoureus, efficiënt analytisch kader om het voor directe injectie specifieke $(r, q)$-voorraadbeleid (herbestelpunt $r$, bestelhoeveelheid $q$) te modelleren, evalueren en optimaliseren. In tegenstelling tot eerdere studies die mogelijk vertrouwen op computerverzwarende simulaties, streeft dit werk ernaar een snelle, nauwkeurige analytische methode te bieden die geschikt is voor ontwerpoptimalisatie.

Methodologie
De auteurs stellen een kader gebaseerd op een discrete-tijd Markov-keten voor om de stochastische dynamiek van satellietuitval en doorlooptijden van lanceervoertuigen (LV) te modelleren. De methodologie verloopt via de volgende stappen:

1. Systeemmodellering:

   • Staatdefinitie: Het aantal satellieten (operationeel plus reserves) in een enkel orbitaal vlak wordt gemodelleerd als een staat $X_k \in \{0, 1, \dots, N_{sat}\}$.
   • Uitvaldynamiek: Satellietuitval volgt een Poisson-proces met snelheid $\lambda_{sat}$. Een overgangsmatrix $P_f$ wordt geconstrueerd om staatverminderingen te modelleren, waarbij wordt aangenomen dat uitgevallen operationele satellieten onmiddellijk worden vervangen door reserves en dat de reserves zelf niet uitvallen.
   • Doorlooptijddynamiek: De doorlooptijd voor grondbevoorrading wordt gemodelleerd als een verschoven exponentiële verdeling ($T \sim \text{Exp}(\mu_{lv}) + \tau_{lv}$), waarbij $\tau_{lv}$ een constante verwerkingsvertraging is.
   • Beleidlogica: Het $(r, q)$-beleid wordt afgedwongen: als het voorraadniveau daalt tot of onder $r$, wordt een bestelling geplaatst voor $q$ reserves.

2. Analytische Afleiding:
   De kernbijdrage is de afleiding van de stationaire staatsverdeling over een volledige aanvullingscyclus ($\pi_{rc}$) door de cyclus te ontleden in twee distincte perioden:

   • Inter-Bestel (IO)-Periode: De tijd van aankomst van een aanvulling totdat het voorraadniveau het herbestelpunt $r$ bereikt. De auteurs leiden de verwachte staatsverdeling $\pi_{io}$ en de duur $k_{io}$ af met behulp van projectiematrices ($C_r^+$) om staten boven $r$ te isoleren.
   • Doorlooptijd (LT)-Periode: De tijd van het plaatsen van een bestelling totdat de nieuwe reserves aankomen. De verwachte verdeling $\pi_{lt}$ en duur $k_{lt}$ worden afgeleid met behulp van de doorlooptijd-kansdichtheid en de uitvalmatrix $P_f$.
   • Cyclusaggregatie: De algehele stationaire verdeling $\pi_{rc}$ wordt berekend als een gewogen lineaire combinatie van $\pi_{io}$ en $\pi_{lt}$, genormaliseerd naar hun respectievelijke verwachte duur.

3. Optimalisatiekader:
   Met gebruik van de afgeleide $\pi_{rc}$ formuleren de auteurs een optimalisatieprobleem om de totale verwachte operationele kosten ($C_{total}$) te minimaliseren, gedefinieerd als de som van productiekosten, voorraadkosten en lanceerkosten. Het probleem is onderworpen aan:

   • Weerstandsbeperking: Het verwachte tekort ($S$) mag een door de gebruiker gedefinieerde drempel $\epsilon$ niet overschrijden.
   • Ladingbeperking: De totale massa van de bestelling ($q \times m_{sat}$) mag de ladingcapaciteit van het lanceervoertuig niet overschrijden.
   • Beslissingsvariabelen: De gehele waarden voor bestelhoeveelheid $q$ en herbestelpunt $r$.

Belangrijkste Bijdragen

• Analyse van Directe Strategie: Het artikel biedt een gedetailleerd analytisch model specifiek voor de directe reserve-aanvullingsstrategie, waarmee eerder werk dat zich richtte op indirecte strategieën wordt uitgebreid.
• Niet-iteratieve Oplossing: In tegenstelling tot de analyse van indirecte strategieën, die mogelijk iteratieve procedures vereist, is de hier gepresenteerde analyse van de directe strategie niet-iteratief. Het lost de stationaire verdeling op via matrixoperaties, wat de computationele complexiteit aanzienlijk verlaagt.
• Kosten- en Weerstandsmaatstaven: Het kader kwantificeert expliciet de afweging tussen operationele kosten (productie, voorraad, lancering) en systeemweerstand (verwacht tekort), waardoor datagedreven ontwerpbeslissingen mogelijk worden.
• Validatie: Het analytische model is gevalideerd tegen Monte Carlo-simulaties over 100 testcases met variërende parameters (uitvalsnelheden, doorlooptijden, bestelgroottes).

Resultaten

• Nauwkeurigheid en Efficiëntie: Numerieke validatie toont aan dat de analytische methode een maximale relatieve fout van minder dan 1% bereikt in vergelijking met 20-jarige Monte Carlo-simulaties (gemiddeld over 1.000 runs). De analytische methode berekent resultaten in minder dan één milliseconde, terwijl simulaties enkele uren vereisen, wat bevestigt dat deze geschikt is als binnenste lus voor optimalisatie.
• Optimalisatiebevindingen:
  • In een basisscenario (uitvalsnelheid $\lambda_{sat} = 0,05$, tekortdrempel $\epsilon = 0,25$) bleek het optimale beleid te zijn $q=2$ en $r=39$.
  • Rideshare versus Volledig Contract: De analyse onthult dat het reserveren van het volledige lanceervoertuig over het algemeen optimaal is, zelfs wanneer rideshare-mogelijkheden beschikbaar zijn, vanwege de specifieke kostenstructuren die worden aangenomen (een premie van 20% voor prijs per eenheid). Rideshare wordt alleen kosteneffectief in een regime met zeer lage uitval ($\lambda_{sat} \in [0,001, 0,008]$).
  • Gevoeligheid: Naarmate de satellietuitvalsnelheid toeneemt, neemt het optimale herbestelpunt $r$ toe (bijvoorbeeld van 39 naar 42) om het vereiste weerstandsniveau te handhaven, terwijl de totale kosten dienovereenkomstig stijgen.

Betekenis en Claims
Het artikel claimt dat het voorgestelde, op Markov-ketens gebaseerde kader een snelle en nauwkeurige tool biedt voor de analyse en het ontwerp van reservebeheerbeleid voor groot-schaal constellaties. De primaire betekenis ligt in het mogelijk maken van ontwerpoptimalisatie door een computationeel efficiënt alternatief te bieden voor simulatie-intensieve benaderingen. De auteurs stellen dat de geoptimaliseerde resultaten voor de directe strategie dienen als tegenhanger voor vergelijking met de indirecte strategie in een begeleidende studie. Bovendien wordt het kader gepresenteerd als uitbreidbaar naar andere constellationconfiguraties (bijvoorbeeld heterogene agenten) en alternatieve aanvullingsbeleid, hoewel het huidige werk zich strikt richt op de directe strategie onder een $(r, q)$-beleid.