Approximation Models for Shared Mobility Rebalancing Under Structured Spatial Imbalance

إليك ملخص تقني مفصل للورقة البحثية بعنوان: "نماذج تقريبية لنماذج إعادة توازن التنقل المشترك تحت عدم التوازن المكاني الهيكلي."

1. بيان المشكلة

تعاني أنظمة التنقل المشترك (مثل مشاركة الدراجات، مشاركة السيارات، وطلب المركبات) من عدم توازن مكاني مستمر. حيث تتراكم المركبات في وجهات الإنزال الشائعة، مما يترك نقاط انطلاق الرحلات في حالة نقص في الإمدادات. يتعين على المشغلين إعادة تموضع المركبات الفارغة لاستعادة التوازن، مما يتسبب في تكاليف تشغيلية كبيرة (إجمالي الكيلومترات المقطوعة بالمركبات، VKT).

• التحدي: يتطلب حساب الحد الأدنى لتكلفة إعادة التوازن بدقة حل برنامج خطي للنقل (LP) بناءً على مسافة إيرث موفر (مسافة واسرتاين-1، $W_1$). وبالنسبة لشبكات الخدمة ذات الحجم الضخم التي تضم آلاف المناطق، يصبح هذا الأمر مستحيلاً من الناحية الحسابية (تعقيد $O(N^3 \log N)$)، مما يجعله غير مناسب للتخطيط التكراري أو للمقارنة المعيارية في الوقت الفعلي.
• الفجوة: النماذج التحليلية الحالية (مثل تلك الخاصة بمسائل البائع المتجول) لا تنطبق على إعدادات إعادة التوازن من متعدد إلى متعدد. كما أن النماذج العشوائية السابقة غالباً ما تفترض وجود متوسط صفري لصافي الإمداد، مما يجعلها تفشل في التقاط عدم التوازن المكاني الهيكلي (الاختلافات المنهجية بين كثافة المنشأ وكثافة المقصد).
• الهدف: اشتقاق نماذج تقريبية ذات صيغ مغلقة لأدنى مسافة لإعادة التوازن في مناطق خدمة مستطيلة، بحيث تكون فعالة حسابياً، ومؤسسة نظرياً، ودقيقة بما يكفي للتخطيط التشغيلي.

2. المنهجية

2.1 صياغة المشكلة

يصيغ المؤلفون نظام تنقل مشترك يعمل في منطقة مستطيلة $\Omega = [0, L] \times [0, W]$ بمساحة $R = LW$ ونسبة عرض إلى طول $\kappa = L/W$.

• المدخلات: كثافة المنشأ $f(x)$ وكثافة المقصد $g(x)$.
• مؤشر عدم التوازن ($I$): يُعرف بأنه مسافة التباين الكلي بين $f$ و $g$:
  $$I = \frac{1}{2} \int_{\Omega} |f(x) - g(x)| dx$$
  هذا الرقم الكمي يحدد حجم عدم التطابق ($I=0$ يعني توازناً تاماً؛ $I=1$ يعني انفصالاً كاملاً).
• المقياس: تُصاغ المسألة كمسألة نقل أمثل لتقليل المسافة الإجمالية تحت مقياسين:
  1. مانهاتن ($L_1$): لشبكات الطرق الشبكية.
  2. إقليدس ($L_2$): للحركة غير المقيدة (الحد الأدنى النظري).

2.2 اشتقاق نماذج التقريب

اشتق المؤلفون نماذج موحدة لكلا المقياسين باستخدام إعادة صياغة تدفق بيكمان، عبر تحويل البرنامج الخطي المتقطع إلى مسألة تقليل تدفق متجهي مستمر.

• نموذج مانهاتن ($L_1$):

  • يتم تفكيك مسألة البعد الثاني إلى مسألتين مستقلتين في البعد الواحد (شرق-غرب، وشمال-جنوب) عن طريق إسقاط كثافة عدم التوازن على المحاور.
  • باستخدام حجج القياس والتحليل البعدي، تبين أن مسافة إعادة التوازن لكل مركبة تتناسب مع مؤشر عدم التوازن $I$، والجذر التربيعي للمساحة $\sqrt{R}$، وعامل شكل يعتمد على نسبة العرض إلى الطول.
  • الصيغة: $W_1^M \approx C_M \cdot I \cdot \sqrt{R} \cdot g_M(\kappa)$
    حيث $g_M(\kappa) = \frac{\kappa + 1}{\sqrt{\kappa}}$ هو عامل الشكل.

• نموذج إقليدس ($L_2$):

  • تم اشتقاقه باستخدام تحليل نمط فوري لمجال عدم التوازن. تهيمن على التكلفة أدنى نمطين من الرتب (شرق-غرب، وشمال-جنوب).
  • يرث النموذج نفس الهيكل وعامل الشكل الخاص بنموذج مانهاتن ولكن مع ثابت معايرة مختلف.
  • الصيغة: $W_1^E \approx C_E \cdot I \cdot \sqrt{R} \cdot g_M(\kappa)$

• إجمالي التكلفة الموحدة:
  إجمالي المسافة المقطوعة (VKT) لفترة زمنية بها $N$ من الرحلات هو:
  $$E[VKT] \approx C \cdot I \cdot (\lambda h) \cdot R^{3/2} \cdot g_M(\kappa)$$
  وهذا يعني أن VKT يتناسب مع $R^{3/2}$ (زيادة خطية في عدد الرحلات $\times$ زيادة في المسافة بمقدار $\sqrt{R}$).

2.3 الحدود النظرية

تثبت الورقة البحثية حداً علوياً صارماً لمسافة واسرتاين-1:
$$W_1(f, g) \leq I \cdot \sqrt{R} \cdot g_M(\kappa)$$
وهذا يؤكد أن نماذج التقريب المستمدة (حيث $C \in (0, 1]$) صالحة فيزيائياً ومحدودة بهندسة النظام وحجم عدم التوازن.

3. المساهمات الرئيسية

1. مؤشر عدم التوازن القياسي: تقديم $I$ كمقياس موحد لقياس عدم التطابق المكاني الهيكلي، وهو متميز عن التباين العشوائي.
2. نماذج التقريب ذات الصيغ المغلقة: اشتقاق نماذج موحدة لمقاييس مانهاتن وإقليدس تدمج صراحةً حجم منطقة الخدمة ($R$)، والشكل ($\kappa$)، وحجم عدم التوازن ($I$).
3. عامل الشكل ($g_M$): تحديد دالة الشكل $g_M(\kappa) = (\kappa + 1)/\sqrt{\kappa}$، مما يثبت أن المناطق المربعة ($\kappa=1$) تقلل من تكاليف إعادة التوازن بغض النظر عن المقياس، بينما تزيد المناطق المستطيلة من التكاليف.
4. التحقق الصارم: تم التحقق من صحة النماذج مقابل 500 حل دقيق للبرنامج الخطي (LP) لكل مقياس ودراسة حالة واقعية، مما أكد قوانين القياس النظرية.

4. النتوات والتحقق

4.1 المعايرة العددية (بيانات اصطناعية)

• مجموعة البيانات: 500 حالة عشوائية لكل مقياس مع مساحات وتوزيعات طلب متفاوتة (منتظم، متكتل، اتجاهي).
• ثوابت المعايرة:
  • مانهاتن: $\hat{C}_M \approx 0.144$
  • إقليدس: $\hat{C}_E \approx 0.119$
  • النسبة $\hat{C}_E / \hat{C}_M \approx 0.825$ تطابق بشكل وثيق التوقع المتماثل نظرياً $\pi/4 \approx 0.785$.
• الدقة:
  • مانهاتن: متوسط الخطأ النسبي المطلق (MAPE) = 18.9%، $R^2_{log} = 0.922$.
  • إقليدس: متوسط الخطأ النسبي المطلق (MAPE) = 25.2%، $R^2_{log} = 0.897$.
• أسات القياس: الأس المستخرج لـ $\sqrt{R}$ هو $\approx 1.02$، مما يؤكد قياس $R^{3/2}$ النظري لـ VKT الإجمالي. أما الأس لـ $I$ فهو دون خطي ($\approx 0.8$)، مما يشير إلى أن حالات عدم التوازن شديدة التركيز لا تتناسب خطياً بسبب تأثيرات الانتشار المكاني.

4.2 المتانة

يظل هيكل النموذج متيناً عبر أنماط الطلب المختلفة (منتظم، متكتل مكانياً، اتجاهي غير متماثل). ومع ذلك، فإن ثابت المعايرة $C$ يختلف بشكل كبير (من ~0.12 للطلب المنتظم إلى ~0.40 للطلب الاتجاهي العالي)، مما يعكس كيف تؤثر بنية الطلب على كفاءة الإلغاءات المحلية.

4.3 دراسة حالة: مدينة نيويورك (يناير 2026)

• البيانات: 31 يوماً من بيانات رحلات المركبات المستأجرة (FHV) عبر 263 منطقة.
• السيناريوهات:
  1. نيويورك (جميع الأحياء): هندسة قريبة من المربعة ($\kappa \approx 1.06$).
  2. جزيرة مانهاتن: هندسة مستطيلة ($\kappa \approx 2.39$).
• النتائج:
  • نجح النموذج في التنبؤ بقيم $W_1$ اليومية باستخدام حلول LP الدقيقة على شبكة طرق OSM.
  • نيويورك: MAPE للمعايرة = 10.9%، MAPE للتحقق = 13.8%.
  • مانهटन: خطأ أعلى (Validation MAPE = 30.2%) بسبب الطبيعة الهيكلية والاتجاهية لتدفقات التنقل التي تبتعد عن الافتراضات المتماثلة للنموذج الأساسي.
  • تأكيد القياس: تتوافق نسبة مسافات إعادة التوازن بين السيناريوهين مع توقع القياس $\sqrt{R}$.
  • تقدير VKT: يقدر النموذج حداً أدنى نظرياً قدره 238,000 كم/يوم من سفر المركبات الفارغة لنيويورك، مما يوفر معياراً للكفاءة التشغيلية.

5. الأهمية والآثار

• المقارنة المعيارية التشغيلية: توفر للمشغلين أداة خالية من الحلول البرمجية (solver-free) لمقارنة أداء إعادة التوازن. إذا لم تستطع خوارزمية التحسين التفوق على الحد الأدنى النظري الذي تتنبأ به هذه الصيغة، فإن التوسع أو الاستراتيجية قد تكون غير مجدية.
• التخطيط الاستراتيجي:
  • تصميم المناطق: تؤكد تقسيم مناطق الخدمة إلى مناطق فرعية مربعة يقلل من تكاليف إعادة التوازن.
  • تحسين التردد: يتيح حساب تردد إعادة التوازن الأمثل ($h$) من خلال الموازنة بين تكاليف إعادة التموضع وجودة الخدمة.
  • تأثير السياسات: يحدد كيف يمكن لبرامج الحوافز التي تقلل من مؤشر عدم التوازن $I$ أن تقلل مباشرة من VKT إعادة التوازن.
• القابلية للتعميم: يمتد المنطق التحليلي لما وراء التنقل ليشمل إعادة تموضع المخزون، والإغاثة الإنسانية، وروبوتات المستودعات، مما يوفر طريقة لربط عدم تجانس الطلب، والهندسة، وتكلفة النقل دون حل مسائل تحسين واسعة النط scale.

6. القيود والعمل المستقبلي

• محدودية المؤشر: يلتقط المؤشر $I$ الحجم ولكنه لا يلتقط الترتيب المكاني المحدد للزيادة أو النقص. يجب معايرة الثابت $C$ تجريبياً لأنماط طلب محددة.
• الهندسة: يفترض الاشتقاق مناطق مستطيلة؛ والعمل المستقبلي سيمتد ليشمل مضلعات غير منتظمة.
• وفورات الحجم: يفترض النموذج الحالي أن التكلفة خطية مع المسافة والكتلة. يمكن للعمل المستقبلي دمج إعادة التموضع الهرمي (مثل الشاحنات للنقل الضخم) لنمذجة التكاليف تحت الجمعية.